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Abstract: Recently, Christodoulou IChrl proved that trapped spheres can form in evo- 
lution, through the focusing of gravitational waves, from a generic initial state. His 
work is the motivation for the present paper, in which we consider the same physical 
problem, using very different mathematical methods. Our approach is based on a con- 
trolled "far field expansion". By a systematic use of scaling symmetries, we regularize 
Christodoulou's singular "short pulse method", rigorously track vacuum solutions by 
the far field expansion and exhibit trapped spheres that first appear deep inside the far 
field region. Our presentation is self-contained. In the final section, we present a de- 
tained outline of the construction of another, more subtle, expansion that allows us to 
continue the solutions beyond the far field region to within any fixed "finite distance" 
from the (expected) singularity. From a methodological perspective, the underlying aim 
of this paper is the development of a general method for constructing solutions to the 
vacuum Einstein equations by controlled expansions. 



1. Introduction 

Formal and controlled (perturbation) expansions are common tools in mathematics and 
physics. In general relativity, see, for example, IIAnReL HBBML [Cha]. This paper is 
a first step in the development of a hybrid method, combining formal expansions and 
simple tools from the theory of hyperbolic partial differential equations (such as energy 
estimates), to construct generic classical solutions to the vacuum Einstein equations, for 
a wide variety of well posed problems with natural small parameters. The purpose of 
this paper is to illustrate the methodology, by carrying it out in all detail for a concrete, 
physically interesting situation. 

Christodoulou |Chr| showed that strongly focused gravitational waves, coming in 
from past null infinity, generate trapped spheres. One of the most important innovations 
of IChrl is the introduction of a small parameter 8 and a picture in which 8 represents 
the duration of a spherical pulse traveling along a null hypersurface. The amplitude 
of the pulse is scaled so that, roughly speaking, the total incoming energy per unit 
advanced time is proportional to i. Christodoulou refers to this picture as his "short 
pulse method". This physical setup triggered our interest in this problem. 
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To illustrate our approach, we shall construct strongly focused gravitational wave so- 
lutions to the vacuum Einstein equations by means of controlled expansions that exhibit 
trapped spheres. 

In Christodoulou's picture, the limit 5 J, is singular. An important preliminary step 
in our analysis is to identify a one parameter group, indexed by 21 ^ 0, of (anisotropic) 
scaling transformations (see, Section [3} that regularizes the limit 5 = 2l 4 J. 0. We 
work almost exclusively in this regularized picture. Our construction of strongly fo- 
cused gravitational wave solutions in the regularized picture is conceptually and techni- 
cally independent from [Chr|. In fact, our formalism is so different, that any discussion 
of the relation to BChrl is deferred until Section|9](see, Remark l93b . We use: 

• A reformulation of the vacuum Einstein equations in terms of an M 31 valued field <P 
(rather than a metric g), a first order quasilinear, symmetric hyperbolic system (SHS) 
for <P, a R 32 valued field of constraints constructed out of <P, and a "dual" linear, 
homogeneous symmetric hyperbolic system (SHS) for <P$. Namely, for each field <P 
satisfying (SHS) and = 0, one may canonically construct a Ricci flat, Lorentzian 
manifold, and conversely, there is a sufficiently small neighborhood of any point on 
any Ricci flat, Lorentzian manifold that arises in this way. The dual system (SHS) is 
an intermediary tool that is used to show that the constraint field vanishes under 
suitable conditions. 

• A unique formal power series solution [<P], in a natural small parameter, to an ap- 
propriate well posed initial value problem for (SHS) that also satisfies the constraint 
equation [ <P$ ] = 0. Here, [ <P ] is the formal power series constraint field constructed 
out of [0]. 

• A standard, local existence theorem and standard energy estimates for (SHS) that are 
used to show that, under suitable (generic) conditions, the formal solution [>P] is 
an asymptotic expansion for a unique classical solution <P to the same initial value 
problem for (SHS) . Energy estimates for (SHS) are used to show that <P^ = 0. 

We emphasize that: 

• The reformulation of the Einstein equations in terms of <P, <P^, (SHS) , (SHS) has 
very useful analytic (first order symmetric hyperbolic systems) and algebraic (all 
the nonlinearities are quadratic polynomials) properties. 

• In our approach, geometry only appears in the formulation of the problem, and in 
the interpretation of the final results. This is in sharp contrast to IChrl and IChrKll . 
in which the very character of the method is consistently geometrical. 

We have tried to present our construction in a transparent form with as many details 
as possible, so that it is accessible to the general reader without any specific background 
in general relativity or hyperbolic partial differential equations. For example, we include 
the derivation of the single (L 2 ) Sobolev inequality that is used. In fact, the discussion, 
up to the formation of trapped spheres, is entirely self contained, apart from the refer- 



ence to [Tay | for a simple local existence theorem for quasilinear symmetric hyperbolic 
systems defined on the product of a time interval with a torus. For these reasons, this 
paper is longer than it might be. We have, however, omitted lengthy, but straight for- 
ward, direct verifications (typical of general relativity) of many equations and algebraic 
identities. We have included an index of notation (Appendix [A). 

We conclude this introduction with an overview of the contents of this paper. Section 
IDgives completely self-contained but unmotivated statements about the purely algebraic 
reformulation of the vacuum Einstein equations. This reformulation is motivated and 
derived in three appendices: 
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• Appendix [B] introduces the underlying general formalism. It is based on work of 
Friedrich (Fr) and the Newman-Penrose complex tetrad formalism ONPI . 

• In Appendices ICl and iDl we choose a gauge, adapted to the problem we are solving, 
which forces the "angular derivatives" in (SHS) to appear only in very special places, 
with the consequence that the equations are easier to control. 

The reader may be put off by the multi-page equations of Section[2] However, he or 
she should not be discouraged, because later, in Propositions l5 . 21 1531 15 .41 we will derive 
relevant /irrelevant forms of these equations, that exhibit the essential constituents that 
have to be treated carefully, and sweep everything else into "generic terms" that we 
don't need to know much about. We will almost exclusively work with this simple, 
transparent form of the equations. 

Section [3] introduces a number of exact symmetry transformations of the equations 
of Section [2] In particular, the global anisotropic scaling 21 (Definition 13.5b plays a 
central role for everything we do. 

In Section 21 we define a two-parameter family of fields M. a ,^, solving (SHS) and 
the constraint equations. For all parameter values a, 21 ^ 0, the field -M a ,2i corresponds 
to Minkowski space. It will be the background for our far field expansion (see below). 
■M a , a is obtained from M,\ t i by applying the scaling symmetries of Section[3] 

In Section|5] an (asymptotic) characteristic initial value problem for the focusing of 
gravitational waves is informally stated. It is motivated by OChrl . 

A formal solution [ <P ] to the characteristic initial value problem is constructed in 
Section|6] To exhibit trapped spheres, it is sufficient to expand directly around past null 
infinity. That is, we make a far (weak) field (power series) expansion in -, morally 
the distance to past null infinity. The expansion is very simple, it can be written down 
explicitly, and one sees trapped spheres directly in the lowest order term. Somewhat 
surprisingly, trapped spheres appear deep inside the formal radius of convergence of 
the far field expansion, that is very close to past null infinity. 

In Section|7] we reduce local existence for (SHS) to a standard result from |Tay[ . We 
then develop energy estimates that enable us to show, in Sections [8] and [9] that, under 
suitable (generic) conditions, every formal solution [ <P ] is an asymptotic expansion for 
a unique classical solution <P with & = 0, which exists at least up to a small, fixed 
fraction of the formal radius of convergence of [ <P ], and contains trapped spheres. 

The formation of a trapped sphere is important, but not an end in itself. Two more 
pressing questions about strongly focused gravitational wave solutions are: 

• Does a horizon form, under suitable conditions? 

• How do these solutions behave as they are continued beyond the "far field region"? 

To see that a horizon forms in strongly focused gravitational wave solutions, one would 
first have to continue them to include all of "future null infinity". We take a step in this 
direction by showing, in Section|9] that, under suitable conditions, our solutions become 
arbitrarily close to the Schwarzschild solution on a strip that extends from the location 
of the trapped sphere all the way out to past null infinity. As soon as a sufficiently well 
controlled expansion around the Schwarzschild /Kerr family becomes available, one 
could conclude that a horizon has formed. 

In Subsection |9.5| we sketch a method for continuing the solutions out of the "far field 
region" using a more powerful expansion in 21. The expansion in 21 around the regular 
limit 21 I is both more fundamental and more subtle. It is an expansion around a two 
parameter family of decoupled, fully nonlinear two dimensional systems. Fortunately, 
these two dimensional systems can be solved "quasi explicitly". Higher order terms 
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in 21 can be constructed and controlled. In the far field region the expansion in — can 
be rearranged into the 21 expansion, and vice versa. We follow the 21 expansion into 
the strong field regime, as far as it will take us (see, the figure in Subsection |9.5l ). The 
rigorous control of the 21 expansion and its consequences will be reported on in all 
detail in another place. It is a second example of our methodology. 

We appreciate the great effort that Demetrios Christodoulou invested over many 
years to nurture the mathematical study of general relativity at ETH Zurich. 

We thank Lydia Bieri, Joel Feldman, Horst Knorrer and Martin Lohmann for en- 
couragement and helpful conversations. 



2. A Reformulation of the Vacuum Einstein Equations 

Our method for constructing solutions to the vacuum Einstein equations has two parts. 
The first is algebraic, the second analytic. Here, we present the purely algebraic part. 
It is a reformulation of the vacuum Einstein equations that is carefully tailored to the 
constructive analytic tools used in the second, purely analytic part. 

This section compresses the intuition and logic of the three leisurely AppendiceslBl 
[C]and|D]into elementary, but very lengthy, totally unmotivated and, to the contemporary 
eye, unsightly, definitions and statements. 

Let 

(x 1 , x 2 , x 3 , x A ) = (C 1 ,^ 2 , u, u) 
be a coordinate system on the open subset U of M 4 , and 

<%) = ($i(x), $ 2 (x), $ 3 (x)) = (e(x), 7(0;), w(x)) 

any sufficiently differentiable field on U taking values in 

K = {(e, 7 ,»)6C 5 ffiC 8 ffiC 5 | e 3 , e 4 , e 5 , 72, 76 eK}, (2.1) 

a real vector space of dimension 3 1 . Throughout this paper, z is the complex conjugate 
of z G C. 

Remark 2.1. Later on, the complex coordinate £ = ^ 1 + i£ 2 will play the role of an 
"angular" coordinate. See, for instance, Remark |4~T1 

Definition 2.1. To any sufficiently differentiable field <P : U — > 1Z, satisfying the condi- 
tions 

(★) : e 3 > 
(★*) : 3(eie 2 )^0 

at every point oflA, we associate three fields F a ^, r a jk, "W a bjk and a complex frame 
F a = Fj/' on IA. Here and below, small Latin and small Greek indices run from 
one to four. The fields are uniquely determined by: 

• r ajk = -r akj andW abjk = -W abkj andW abjk = -W bajk . 
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The matrix indices (ab), (jk) run over the ordered sequence 
(12) (31) (32) (41) (42) (34) 
The complex frame is written as: 
• (F 1 ,F 2 ,F 3 ,F 4 ) = (D, D, N, L) or, equivalently, 



D 



ei 



'■2 



<>£.- 



N = e 4 ^1 + e 5 



9C 2 



da 



e 31 



The vector fields N and L are always real. 
Proposition 2.1. The field W a bjk has the symmetries 



w 



akij 



Wabjfe = W jkab 

where the matrix g a (, and its inverse g ab are given by 



, gah 



/0 1 
10 
-1 

\0 -1 



g ai w abJk 



/o 1 

10 
-1 

Vo -1 



(2.2) 



(2.3) 



That is, *W abjk has the algebraic properties of a Weyl field. 
Proof. By direct inspection. □ 

Remark 2.2. Later on it will be important to drop (**). That is, to allow the frame to 
collapse. 

The next definition singles out an important class of fields <P(x) and relates them to 
Ricci-flat Lorentzian manifolds, that is, vacuum spacetimes. 

Definition 2.2. Afield <P : U — > 1Z is a vacuum field when: 
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• Conditions (*) and (* *) are satisfied at every point of 14. 

• The Levi-Civita connection for the complex linear metric gonlA is given by 

&(V Fa Fj, F k ) =r ajk 

where g(F a ,F b ) = g ab (see, d2J]lj and 
g(V Fo ^, F k ) d = f i ( - g (F a , [-Fj,-Ffe]) + s{F k , [F a , Fj]) + g(F J: [F k ,F a ])) 

• The Riemann tensor for the Levi-Civita connection is given by 

Rabjk = g([VF J ,Vj ilb ]F 6 - V [FjlPh] F b , F a ) = W oMJk 

Consequently, the Ricci curvature vanishes, since 'Wabjk is traceless. 

• 77ie coordinate functions u and u are both solutions to the eikonal equation. More 
precisely, e%N and L are null geodesic vector fields that are minus the gradients of 
u and u. 

Remark 2.3. For any 7^-valued field <P satisfying the conditions (*) and (* *) the metric 
given by g(F a , Ft,) = g a & is real in the sense that g(X, Y) is real whenever X and Y 
are real vector fields. Over the reals, it has signature (— , +, +, +). Bear in mind that 
g a b are the components with respect to the complex frame F a . If <P is a vacuum field, 
then (U, g) is a Ricci-flat Lorentzian manifold, that is, a solution to the vacuum Einstein 
equations. 

Remark 2.4. It is natural to ask whether all Ricci-flat Lorentzian manifolds arise, at least 
locally, in this way. The answer to this question is given in Appendices ICl and iDl 

Remark 2.5. Assume <P(x) is a vacuum field. Declare L + N to be future directed. 
Let Su^u be the intersection of the level sets of u and u, which by Definition I2.2I are 
null hypersurfaces. The traces of the future-directed second fundamental forms of 
relative to the level sets of u and u are given by g(V dL, D) + g(V^jL, D) = 2j2 
and g(V D N,D) + g(V^N, D) = 2^. By dTTT ), they are real, as they should be. By 
definition, Su,u is a trapped surface when 72 and 76 are strictly negative everywhere 
on S UyU . Equivalently, S vu is trapped if an infinitesimal shift of S v>u along either L 
or N (both future-directed null vector fields, orthogonal to Su fU ) induces a pointwise 
decrease of the area element. 

The basic examples of (closed) trapped surfaces in a vacuum spacetime are the spherical 
SO(3) orbits inside the horizon of a Schwarzschild spacetime. Closed trapped surfaces 
appear in the the formulation of Penrose's incompleteness theorem, see OPenl . 

We need a criterion for a field to be a vacuum field. To this end, we make two more 
definitions. 

Definition 2.3. Suppose, condition (*) is satisfied at every point of the domain IA. Let 
<P(x) = (e(x), j(x), w(x)) : IA — > 1Z be a sufficiently differentiable field, and let the 
weights Ai, A2, A3, A4 be strictly positive functions on IA. The Quasilinear Symmetric 
Hyperbolic System (SHS) for the field $(x) is 



A($)<P = f(<?) 



(2.4) 
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Here, A(<P) = Ai (<P) © A 2 (#) © A 3 (<£) is the first order, matrix differential operator, 
with coefficients that are affine linear functions (overU) oftp, given by 



Ai(#) = diag(L, L, N, L, L) 

A a (*) = diag (L, L, L, L, /V, /V, TV, L) 

/X ± N AiL> 

Ai-D Ai7: + A 2 iV X 2 D 
A 3 {<P) = A 2 ^ A 2 L+_A 3 iV 
X 3 D 

\ 



\ 



X 3 D 
A 3J L+_A 4 iV X 4 D 
X 4 D X 4 L J 



(2.5a) 
(2.5b) 

(2.5c) 



Observe that the "angular" operators D, D only appear in A 3 (<P). Also, 



f(#) = fi(«P)©f 2 («P)©f 3 («P) = 



t53> 



is f/ie quadratically nonlinear vector valued function given by 



fii = S 



r j2 



r j3 



- ei7 2 - ei7i 

- e 2 7 2 - e 2 7i 
e 3 2 5R 78 

2 3? ( - ei74 + ei7 5 
2 3? ( - e 2 74 + e 2 7 5 
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7578 + 7578 - 7673 + 7774 - 7775 + w 4 
2 



- 7676 - 76 2 3? 7s - 1 77 

- 27 6 77 - 3777 8 + 777 8 - w 5 

- 27574 + 27 3 7 5 + 7 3 7 3 + 7474 - 7475 - 7574 + w 3 

- Ai (371 w 3 - 6y 3 W2 + 7ewi — 47 8 wi + 474 w 2 ) 

- Ai (47 2 w 2 + 474 wi + 75 wi) 

- A 2 (27iw 4 - 37 3 w 3 + 27 6 u> 2 - 27 8 u> 2 + 37 4 w 3 ) 

- A 2 (37 2 w 3 + 27 4 u; 2 + 275102 + 77W1) - A3 (71105 + 3lew 3 + 2j 4 W4) 

- A 3 (27 2 u; 4 + 37 5 u; 3 + 277102) - A 4 (37 3 io 5 + 47 6 io 4 + 27 8 w 4 + 7 4 w 5 ) 

- A4(7 2 u'5 — 274 Wi + 475 W4 + 377 w 3 ) 



s 



Definition 2.4. Let <P(x) = (e(x), "f(x), w(x)j :U — > TZ be a sufficiently diffe rentiable 
field, and let Ai, A2, A3, A4 be strictly positive weight functions on IA. The associated 
constraint field 

/V 

5 : 

(2.6) 



&{x) = (*»(*), 4(x), 4(x)) = (t(x), u(x), v(x)) = 




on U taking values in 

TZ = {(t,u,v) G C 5 



is given by 



-23(L>(ei)+e l73 + e 1 7 4 ) 
- 2 3 (D(e 2 ) + §273 + e27 4 ) 
£>(e3) - e 3 7 3 + e 3 7 4 + e 375 
D(e 4 ) - iV(ei) - ei7 7 - ei7 6 + ei7 8 - ei7 8 
D(e 5 ) - N(e 2 ) - e 2 j 7 - e 2 ~/ 6 + e 2 7 8 - e 2 7 8 

£>(7 2 ) - £>(7i) + w 2 - 37174 - 7173 - 7 2 7 3 + 7274 
I>(74) + D(<y 3 ) -W3 + 27374 + 7177 - 7276 + 7373 + 7474 
D( l3 )-D(j 4 ) + N( l2 ) 

- W 3 + 27374 - 27474 + 71 77 + 7276 - 7278 " 7278 

D(7 4 ) - £>( 73 ) - JV( 7l ) 

+ 37i7s + 27373 - 27374 - 7176 - 7i7s ~ 7277 

^(75) - £(77) - 7277 ~ 7475 + 7575 - 7573 ~ 767l 
L(7 6 ) - Z?(7 5 ) - w 3 + 7177 + 7276 - 7375 - 7574 - 7575 
D(j s ) ~ N(7s) - 27376 + 737s - 7 3 7 8 + 7477 + 7674 
■2V (74) + L»(7s) + w A - 27377 + 7476 + 7478 - 747s + 7774 
D(7r) - D(-f 6 ) - w 4 + 37377 + 7476 - 7e73 + 7774 

Ai ( D(wi) + L(w 2 ) + 4-y 2 w 2 + i-j^wi + 75W1) 

A 2 ( D{w 2 ) + L(w 3 ) + 372W3 + 274W2 + 275W2 + 77W1) 

A 3 ( D(w 3 ) + L{wa) + 272W4 + 375W3 + 277W2) 



(2.7a) 



(2.7b) 



(2.7c) 



Proposition 2.2. Suppose, conditions (*) and (* *) are all satisfied at every point oflA. 
Then, the field <P(x) is a vacuum field if and only if there exist strictly positive weight 
functions \\, X 2 , A3, A4 onlA such that <P(x) is a solution to the quasilinear symmetric 
hyperbolic system (SHS) and the constraint field <P^ (x) = everywhere on U. 



Proof. Follows from the Appendices 



See,Propositions|OEI]IE2E3] □ 



Remark 2.6. Proposition ^. 21 together with Definition ^. 21 is a reformulation of the vac- 
uum Einstein equations. When does the solution <P to a well posed problem for (SHS) 
also satisfy <P f = 0? 
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Proposition 2.3. Suppose, that condition (*) is satisfied and there are strictly positive 
weight functions Ai, A2, A3, A4 on U, such that <P(x) is a C 2 solution to (SHS) . Then, the 
constraint field <5" is a classical solution to the "dual" Homogeneous, Linear (over 
K) Symmetric Hyperbolic System (SHS) ; 



= f (<2>, d^) ^ 



In particular, if the data for any well posed problem for the system (SHS) vanishes, then 
the constraint field <2>"(x) vanishes everywhere. Here, A(<?) = Ai © A2 © A3 is the 
first order, matrix differential operator 



A, 
A 2 

A3 = 



diag (L, L, N, L, L) 
diag (L, L, L, L, N, N, L, L, N) 

(t: n +J- 2 l t- 2 d 

A2 A2 A3 A3 



(2.8a) 
(2.8b) 

(2.8c) 



and f(<P, dx 1 !)) is a linear (over R) transformation acting on 3>* — (t, u, v): 



fll\ /fl2 

f (#, a #) <Z>» = (fi©f 2 ©f 3 )^ = I : **ffi : 

?5l/ \% 



ft 



13 



V? 



:i.V 



where <P has been suppressed on the right hand side, and 



- 272^+23(^3^)+ 2 3 (eiui) 

- 2 72 i2 + 23(t 3 ^) +23(e 2 «i) 

(-76 + 2 78 )t 3 - 7 7 t 3 - i 4 f§T - *5§§f + e 3 u 7 + e 3 u 8 
*(73 - 74 + 7s)ii - 72*4 - 7i*4 - £37^7 + eiw 3 - ei¥4 + eiit 5 - ei¥ 6 
. *(73 - 74 + Js)h - 72*5 ~ 7i*5 - h^t + £2^3 - e 2 U4 + e 2 u 5 - e 2 tt 6 
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' -*3^+^-*ysiui + 7iui + ^«i 

- *3fr? ~ - 74M1 + 75"1 - 73^1 - 272^2 - 71 "5 - 72"6 - j;V 2 
~ ^TH + + 73«1 - 74M1 + 75^1 

- 272M3 + 7l"4 + 7l«4 - 71^5 + 72"6 - XI "2 

+ *3^f _ + (73 - 74 + 7 5 ) u l + 7l«3 + 7i"3 - 272W4 + 72«5 - 7l«6 

- ^fff - h - 77U3 + 76 "4 - 27 6 u 5 - 27 8 u 5 + 27 8 u 5 

+ 77"6 + 77"6 + 75"7 - 75""8 - 74W9 + 75 "9 + 73 u 9 
* + *4§^T + *5 |p + 76^3 - 77"4 + 7~7 U 5 + 77^5 

- 276U 6 + 75"7 - 75«8 + 73^9 - 74 M 9 + 75 u 9 + J^V 2 

- h^- ~ J3U2 + 1a u 2 - 75 "2 + 73 U 3 + 74U4 - 75"4 

- 74"5 - 273U6 + 74 u 6 - 72W7 - 71 "8 - J^Vi 
~ *3^f + 73 w 2 - 74"«2 + 75^2 - 74U3 + 75M3 - 73W4 

+ 273U5 - 7 4 u 5 + 7 4 u 6 - 7^7 - 7 2 fi 8 + j^v 3 

+ - +*5§^ - t B g* - 377^7 + 76«7 

+ 76 "8 + 77"8 ~ 376 U 9 ~ 27s Ug + 77U9 + T7U3 



fl3 #» = - ih^ ~ ih^ + t 3 ^ - ~ 

+ 3w 3 Ui + 2W2U2 + 4W2U3 — 2w2Ue + 4u>iU8 + WiUg 

+ [ - £72 - ^76 + ^78 + ^7 8 + (^) 2 N (AO + (^) 2 L (A 2 ) ] Vl 
+ W'D (A 2 ) + ^73 ~ ^74 - ± 2 %]v 2 - ^71^3 

f 23 #« = - »tl^ - + t 3 ^ - t 4 |ff - i 5 fp 

+ 2?l>4Ui + 3W3U3 — 3W3U6 + 2W2^8 + 2u>2M9 

+ [-^~^ + (± 2 ) 2 D(X2)]v 1 

+ [ ~ ^72 - ^76 + ^78 + ^78 + (i) 2 ^ (A2) + (^) 2 i (A 3 ) ] v 2 
+ [i73-i74-^75 + (^) 2 ^(A 3 )],3 

133*^ — — rt iacT — «t 2 -gj2- + t 3-g^ C4 ajr — 15 gjf 

+ W5U1 — 2W4U2 + 2W4U3 — Aw^Uq + 3lf 3 Ug 
-&7V 1+ [-± l5 + (±) 2 D(\ 3 )]v2 

+ [ - £72 - ^76 - ^78 + ^78 + (A 3 ) + (^) 2 L (A 4 ) ] v 3 

Proof. Follows from the Appendices iBl lCllDl See. Proposition lD^4l □ 
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e 4 diag(0, 0, 1, 0, 0)© e 4 diag(0, 0, 0, 0, 1, 1, 1, 0)© 



Remark 2.7. Write A(#) = A"^ = (A^ 1 © A 2 " © A/)^. Explicitly, 

^ie 4 Aiei \ 
Aiei A 2 e 4 A2ei 
A 2 ei A 3 e 4 A 3 ei 
A 3 ei A 4 e 4 A 4 ei 
\ A 4 ei j 

A 3 = e 3 diag(l, 1, 0, 1, l)ffie 3 diag(l, 1, 1, 1, 0, 0, 0, l)ffle 3 diag(0, X u A 2 , A 3 , A 4 ) 

A 4 = diag(0, 0,1,0,0)© diag(0, 0,0, 0,1,1,1,0)© diag(Ai, A 2 , A 3 , A 4 , 0) 

The matrices A p are Hermitian matrices whose entries are linear (over R) functions of 
<Pi = e. The matrix A 3 + A 4 is strictly positive definite by the requirement e 3 > of 
condition (*). We see that, (SHS) is truly a quasilinear, symmetric hyperbolic system. An 
entirely similar discussion applies to (SHS) . See OJohnl . | |Tay| , for linear and quasilinear 
symmetric hyperbolic systems (in the sense of Friedrichs). 



Remark 2.8. By definition, the fields <P and & take values in TZ = R 31 and TZ = 
respectively. The left and right hand sides of (SHS) are in TZ, because f 3 i, f 4 i, fsi, f 22 , 
f6 2 , marked by * , are real. The left and right hand sides of (SHS) are in TZ. In other 
words, (SHS) is equivalent to a real quasilinear symmetric hyperbolic system for an R 31 
valued field, and (SHS) is equivalent to a real linear homogeneous symmetric hyperbolic 
system for an M 32 valued field. 

Remark 2.9. Let *p be the parity transformation (<p • <?) (x) = (- 1) A <P{x) with 

A = diag(l, 1, 0, 0, 0) © diag(0, 0, 1, 1, 1, 0, 0, 0) © diag(0, 1, 0, 1, 0). 

The field <P • ^ solves (SHS) if and only if solves (SHS) . The constraint (<p ■ = 
if and only if = 0. Clearly, *}3 o «p = Identity, and $ splits naturally into 1 1 El- 
even and 7 <}3-odd components. If the ^-odd components of <P vanish at x 6 U, that is 
(CP ■ $)(x) = ${x), then (SHS) implies that L(e 4 ) = L(e 5 ) = at x. 

Proposition 2.4. Suppose (*). Set all the ty-odd components, e\, e 2 , 73, 7 4 , 75, w% iu 4 , 
and the two ty-even components e 4 , es of the field <S> equal to zero, and introduce the 
field & = e 3 © (71,72, 76, 77, 7s) © (wi, W3, W5). In this case, the frame collapses to 
D = 0, N = and L = e 3 and the system (SHS) reduces to 

A(#)# = (subSHS) 
A($) = N © diag(L, L, N, N, L) © diag(Ai2V, A 2 i + A 3 iV, A 4 L) 
f ($) = f 31 © (f 12 , f 22 , f 62 , f 72 , f 82 ) © (f 13 , f 33 , f 53 ) (see, Definition^ 

It is, separately for each £, a quasilinear symmetric hyperbolic system for <S> in the (u, u) 
plane. The components t, U\, u-j, Us, Uq, V\, w 3 o/<2>" vanish. (SHS) reduces to a linear 
(overW), homogeneous symmetric hyperbolic system for <P " = (u 2 , Us, u 4 , U5, Uq)®V2- 



Corollary 2.1. Suppose, <P is a solution to any well posed problem for (SHS) , such that 
all its ?(5-odd components and e 4 , e§ vanish initially. Then, they vanish everywhere. 
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3. Symmetries 

A field transformation S with respect to the open subsets IA, W of R 4 consists of 

• a diffeomorphism from IA to IA' , 

• a map from fields <P = (e, 7, w) : IA — > 72. to fields <£' = (e', 7', «/) : W' — > 72, 

• a map from strictly positive weight functions 71 = (Ai, A2, A3, A4) on IA to strictly 
positive weight functions A' = (X[, \' 2l X' 3 , A4) onW. 

Let x, .t' be Cartesian coordinates on IA and We write 



x' = S-x $'(x') = (S-$){x') A'(x') = (S-A)(x') 

In this section, S will always be linear over real valued functions in its action on <P and 
A. That is, 



(S ■ (/#)) (z') = • z') (5 ■ $){x>) , (5 ■ (M))(x') = fiS- 1 ■ x>) (S ■ A)(x') 



for all / e C(W, M). Therefore, 5 acts pointwise. For this reason, it suffices to make a 
local analysis. For the rest of this section, we make the assumption that x' = S ■ x is a 
local diffeomorphism on E 4 . With this understanding, it is unnecessary to specify the 
domains IA and W . 

Definition 3.1. Afield transformation S is a field symmetry if: 

• (*) and (**) are preserved (see, Definition \2.1\ . 

• <P satisfies (SHS) on IA if and only if S ■ <P satisfies (SHS) on IA' . 

• 4>" vanishes on IA if and only (S ■ <P)^ vanishes on W . 

It is implicit in the last two statements that the weights A appear on IA and the weights 
S ■ A appear on IA' . For afield symmetry S, it follows that <P is a vacuum field on IA if 
and only if S ■ <P is a vacuum field on IA' . 

As in Section 12 let x = (x 1 , x 2 , x 3 , x 4 ) = £ 2 ,u, u). We now define a number of 
field transformations. 

Definition 3.2. Angular coordinate transformation €. Let <£}, € 2 be real functions. 






£ e B {x)\ x 
3=1 
2 A 

£ ffs-0) e B +3(x) 



x' 



x' 



A 



A 



1,2 



1,2 



B = \ 



(e 3 , 7 , w ,/i)(e:- 1 -x') 



We w;7Z also use the notation <£(£) = £ 4 (£) + «C 2 (£), w/zere £ = £ 4 + z£ 2 . 
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Definition 3.3. U(l) transformation 3. Let ( = ((a; 1 , a; 2 ) G U{1). 
x' = 3 • X = X 



( 



#V)= (3-^)(x') = C A ^) 



1 









1^(0 



(x) © 



1 






Vlr^oy 



a:— 3 — 1 



A = diag(l, 1, 0, 0, 0, 2,0,1,-1,-1,0,-2,0, 2,1,0,-1,-2) 
A>(x') = C3-A)(x')=A(x)\ x=3 _ Ux , 

Here D(Q = (ej ^ + e 2 ^)(C) and 7V(C) = (e 4s |r + e 6 ^r)(C). 
Definition 3.4. Global Isotropic Scaling Let 3 > be a constant. 

X -,J ' X {X' , X , -\J^C - ) 

#V) = (3-*)CO = 3 A *(») L^-i.^ 

A= (-1) diag(l,l,0, 1,1, 1,1,1,1,1,1,1,1, 2,2,2,2,2) 
A'(x') = (3.A)(x')=A(x)\ x ^- 1 . x , 
Definition 3.5. Global Anisotropic Scaling 2L Lef 21 7^ 0. 

*V) = (St.*)0O = a^Wl^,., 

A= (-l)diag(2,2,0,3,3, 0,0,1,1,1,2,2,2, 0,1,2,3,4) 
A'(x') = (21 ■ A)^) - diag(l, 2l 2 , 2l 4 , 2l 6 ) A{x) \ x= ^. x , 



The transformation 21 plays a central role in this paper. For a sample calculation, see 
the proof of Proposition l4.ll 

Proposition 3.1. £, 3, 3, 21 are field symmetries. 

Proposition |3.1| is proven in Appendix|E] We will now define an additional transforma- 
tion. It is a composition of two field symmetries, and therefore itself a field symmetry. 



Definition 3.6. Pole-Flip transformation Flip Q . Let a 7^ be a constant. The Pole- 
Flip transformation Flip Q is the composition of a U(l) transformation and an angular 
coordinate transformation. Precisely, 

Flip Q = 3 o € where C(0 = ~ |- <C(0 = f 



and£, = i 1 +ii 2 and€(0 = C 1 ^) + i£ 2 (£). 
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Remark 3.1. With this choice of ( and <£, we have C" 1 o 3 ° £ = 3 _1 and £ o £ = 
Identity. Therefore, Flip Q o Flip Q = Identity. The field symmetry Flip Q will be used 
to match constructions between two "angular" coordinate patches. For Minkowski, £ 
will be stereographic coordinates (scaled by a) based on the north and south poles. 



4. The Doubly Scaled Minkowski Field M a ,% 

Fix the coordinates (x 1 , x 2 , x 3 , x A ) = (£*, £ 2 , u, u) as in the preceding sections. For all 
pairs fx, A > 0, set 



Strip(/i,A) = K 2 x (0,/i) x (-oo, -A 1 ) 
Strip^ = Strip(oo, oo) = R 2 x (0, oo) x (— oo, 0) 



(4.1) 



Definition 4.1. For all a, 21 ^ 0, let M a ,% ■ Strip^ — > 1Z (see (12. U ) be the field 



M a : 



I Pali e aM 
1 






/ \ 

Pa,m ^o,,% 

Pa)n X »M 


-PaM 



\ 



f \ 







/ 



(4.2) 



//ere £ = £ x + i£ 2 anof 

ft,,a(u,u)=a 2 «-u e a , a (e) = f (1 + f^l^l 2 ) A o , a (0 = -g£ (4.3) 
We will often consciously suppress the subscripts a, 21 on the functions p, e and A. Set 



The decomposition - 



S(u,u) = ^+u. 
■Mr will be used over and over again. 



(4.4) 



Proposition 4.1. For all a, 21 ^ 0: 

(a) M a ,<n = (Co 21) ■ Mi,i on Strip^, where £(£) — a£. 

W .Ma,a = Flip f ■ A^ a , a on stri Poo n {£ 7^ 0}. 

(c) M a ,% = Z ■ M a .% on Stri Poo , for all Z > 0. 

fii) A^a.a is a vacuum field on Strip^ (see, Definition ^. 2h 

(e) The Lorentzian manifold associated to M a .% is isometric to the open subset of 
Minkowski space given by (14.51 l below. For this reason, we refer to M a ,% o,s the 
doubly scaled Minkowski field. 

Here, C, 21, Z and Flip Q are field transformations defined in Section\3\ 
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Proof. For (a), set C a = diag(a, a, I, a, a) © Is © I5 and let A be the matrix in 
Definition |3.5l Then, 

((£ o 21) • Mi A ) {x") = (C • (21 • Mi,i))(/) 

= c a (a-Aii,i)( a /)| x , =e _ I . ie „ 

= c a A A4 1)1 ( a! )| )e=a _ I . a , =a _ I . (e _ 1 . !e „ ) 

= C 2i A ^l ljl (f^ )M ",^ U ") 
= Af a , a (a/') 

Similarly for (b) and (c). Part (d) is also verified by direct calculation. It suffices to check 
(d) for A4i.ii because the general case follows from (a) and Proposition 13. II Recall 
that the definition of a vacuum field is independent of the choice of weight functions 
Ai , . . . , A4. For (e), see Remark PTTl below. □ 



Remark 4.1. The Riemann curvature tensor of the Lorentzian manifold associated to 
the vacuum field Ai a .% = (e, 7, 10) on Strip,^ vanishes, because w = (see, Definition 
It is isometric to the open subset of Minkowski space given by 



{(*°, 



< IXI, X 



^{0}x{0}x [0,oo)} 



(4.5) 



where (X°, X) are the standard Minkowski coordinates, and 



X 1 



V2\QL\ 



(21 V 




2FV 



/ 22l f l \ 

n. ^ 



V2I21I 1 + f |£| 2 



221^2 



The level sets of u = 2"3|2l| (X° - |X|) < and u = 2-i\±\(X° + |X|) > are 
null hypersurfaces. They intersect in a standard sphere of radius |X| = 2~ 2 with 
the north pole removed, on which — £ is the standard stereographic coordinate system. 
The southern hemisphere corresponds to |£| < 1 4 1. 

Remark 4.2. The limit A4o,o = li m aj.o A4 21,21 exists on Strip^. By taking the limit 
of (d) in Proposition 14.11 it is a solution to (SHS) with M = 0. Observe that the 
associated frame is degenerate, because D = here. (See, Proposition ^. 41 ) 

Remark 4.3. For each a / 0, the limit Ai a .o = lim^o M- a ,% exists on Strip^. By 
taking the limit of (d) in Proposition 14. 11 it is a vacuum field. The Lorentzian mani- 
fold associated to M. a is isometric to the open subset of Minkowski space given by 
{|X°| < |X|, X° + X 3 < 0}: 



^(X° + X*) = u ^(X°-X*) = u 



V2 X 



X 1 
X 2 



= -V2 




Observe that — (X°) 2 + |X| 2 = —2uu. The field A4 a ,o is independent of £ and u, that is, 
translation invariant in these directions. The level sets of u and u are null hypersurfaces. 
They intersect in standard Euclidean planes. 
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5. The Far (Weak) Field Ansatz 

We now formulate a physically interesting initial value problem, motivated by llChrl . 
for (SHS) (see, Section[2]). It is set up so that <P$ always vanishes, that is, <P is a vacuum 
field. Under appropriate conditions, the corresponding vacuum spacetime (see, Remark 
I2.31 i contains trapped spheres. 

Fix the coordinate system (x 1 , x 2 , x 3 , x 4 ) = (l- 1 ,u,u) on the infinitely wide 
strip Strip^ C M 4 , see ( 14.11 ). and make the far ("past null infinity") field Ansatz 

*(a:) = M aM (x) + u- M V(x) (5.1) 

on U. See, (14.2b for the definition of A4 a ,<2i{x). Here, 

M = diag(2, 2, 2, 3, 3) © diag(l, 2, 2, 2, 2, 2, 2, 3) © diag(l, 2, 3, 4, 4) 
ff(x) = (^i(x), ^ 2 (x), #3(1)) = (/(x), z(x)) G ft 

Our basic Ansatz ( 15. Il l, Minkowski plus asymptotically small corrections (assuming, 
IP = 0(1) as 11 — ► —00), is completely naive. The only subtlety, lies in the choice of 
the diagonal matrix M that prescribes the far field asymptotics of the system, and is 
ultimately a statement about the physics of propagating gravitational waves. Anyone 
who has made formal or rigorous perturbative calculations or constructions in classical 
or quantum physics knows from experience that one must "play with the expansion" 
until, "one sees what is going on". We have followed this traditional route to the matrix 
M, However, the only real justification is that it works. 

Definition 5.1. Let DATA(£ , u) ■ K 2 X (0, 00) — > C be smooth and vanish when u < u , 
for some fixed u > 0. Set, step by step down the left column and then down the right 
column, #(0)(£,u) = (/(0)(f,u),u>(0)(£,n),2(0)(£,u)) eTZequalto 



wi(0) = 


DATA 


*s(0) = 





w 7 (0) = 




wa(0) = 




zi(0) = 


-|z^(o) 


w 4 (0) = 




^2(0) = 


2(e^ + 2A)^ 1 z 1 (0) 


w 6 (0) = 





MO) - 


2(e^ + A) d^z 2 (0) - d?(wr{0)*l(P)) 


/i(0) = 


ea-V(o) 


24 (0) - 


2 e£ 9-^3(0) -2d„> 7 (0)z 2 (0)) 


/ 2 (0) = 


-ze^V(O) 


w 3 (0) = 


-0-^(0) -A0"V(O) 


/3(0) = 


-5Rw 8 (0) 


w 6 (0) = 


-9-^2(0) 


,A(o) = 


-4ea~ 1 5Rcj 5 (0) 


wa(0) = 


9-^3(0) -4i5- 1 9(Acj 5 (0)) 


/ 5 (0) = 


4ed- 1 9w 5 (0) 



w/iere e = e a ,a A = A a ,2i are defined in ( 14.31 ) a«c/ ^ = \ (wt ~ ^W 1 )' an d 



{d^g)(u) = / dtt'fldi'). (5.2) 
Jo 
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(Asymptotic) Characteristic Initial Conditions: Informally, the initial conditions for 
the field <P on Strip^ are: 

lim &(£,u,u) = &(0)(£,u) (5.3a) 

u — > — oo 

u, u) = when u < u (5.3b) 

Remark 5.1. The function DATA is arbitrary, there are no hidden symmetry assumptions. 
For the construction of classical solutions, we will assume a transparent smallness con- 
dition (see, Theorem l8.U . 

Definition 5.2. Let a, 21 ^ 0. Let 

[M,, a ]W = £r=o(^) fe M aM (k)(t,u) 
be the formal expansion in — for the Minkowski vacuum field -M a ,a, see ( 14.2b , in which 

\\\ = + [S] = -EZ Q (z) k % 2{k+1) u k+1 - (5.4) 

In Section|6]we construct a unique formal series [!^](x) = X^o(«) fc ^(&)(£>2£) on 
Strip^ C M 4 such that [<P] = [M a ,% } + u~ M [<P] is a formal power series solution 
to (SHS) satisfying d5.3ab and (I5.3bb . In particular, &(k)(£,u) = when u < w for 
all k > 0. By construction, the associated formal constraint field [ ) vanishes "as 
u — > — oo", and therefore identically because it is a formal solution to (SHS) . 

Remark 5.2. Equation (15 . 3bb stipulates that <P coincides with the Minkowski vacuum 
field M a ,<& when u < u Q . On the other hand, d5.3ab is an asymptotic initial condition 
at "past null infinity" u — > — oo. Where does the prescription for •P'(0)(^, u) in Defini- 
tion [5TT] come from? The system (SHS) (see, Definition 12. 31 l for [<P] is equivalent to a 
recursive hierarchy of inhomogeneous linear equations for •P'(fc), k > 1, and equations 
for the initial term ^(0). Additional equations for 5^(0) have to be satisfied so that the 
formal constraint field ] (see, Definition 12.41 ) vanishes "as u — > — oo". Altogether, 
there are just enough equations to completely determine all but one of the components 
of ^(0), namely wi(0) = DATA. Our prescription amounts to an "upper triangular" 
arrangement of these equations that generate tf'(O) out of DATA. See, Section[6] 

The goal for the rest of this section is to bring the equations of Section 2 into a 
relevant /irrelevant form that exhibits the essential constituents that have to be treated 
carefully, and sweeps everything else into "generic terms" that we don't need to know 
much about. 

Proposition 5.1. In this proposition, ignore Definition \2.4\ and regard <P(x) and <P^ (x) 
as independent, sufficiently differentiable fields on Strip^ with values in 1Z and 1Z, re- 
spectively. Set 

M = diag(2, 2, 2, 3, 3) © diag(l, 2, 2, 2, 2, 2, 2, 3) ffi diag(l, 2, 3, 4, 4) (5.5a) 
E = diag(4, 4, 4, 6, 6) © diag(2, 4, 4, 4, 4, 4, 4, 6) diag(0, 0, 0, 0, 0) (5.5b) 
M* = diag(2, 2, 2, 3, 3) diag(2, 2, 2, 2, 2, 2, 3, 3, 3) diag(0, -1, -2) (5.5c) 
= diag(4, 4, 4, 6, 6) © diag(4, 4, 4, 4, 4, 4, 6, 6, 6) diag(2, 2, 2) (5.5d) 
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and 

$(x) = M a ,%(x) + ir M W(x) W(x) G K (5.6a) 

^(x) = u- M> ^(x) G ft (5.6b) 

\ J ( x ) = u 23 j = 1, 2, 3, 4 (see, Definition\Z3l (5.6c) 

The systems (see, Section^ A(<P)<P = f (#) and A(#) ^' = f c^) /or ^ and 
<p" are equivalent to the following systems for W and W^: 

A a ^{x,<J?)<P = ia,*(x,V) (5.7a) 
Aa.afx, = f a , a (x,f,^)f« (5.7b) 

w/zere A a , a (x, <f) = A£ a (x, ) gfjr one/ A 0) a(ar, ^) = A^x, and 

A a,siM = u E {u- M A»($)u- M ) (5.8a) 

£,,*(*, «?) = - A^)(g§^- M ) # + f(*) - A^(<P)^X Q , a ) 

(5.8b) 

A£ a (x, = u E * (u- Mt A»(<l>) u- M> ) (5.8c) 
f a<Ql (x,W,d^)=u Et - Mf [-A^)(^u- Mt ) +l^,d x ^)u- Mi ) (5.8d) 

7/7 d5.8t , <P, d x <S> have to he expressed in terms ofW, d m \P using d5.6a| l. We will sometimes 
drop the a, 2t and write A(x, W), f(x, X P), A(x, f (x, d x ^). They are notationally 
distinguished from A(<2>), f (<P), A(<P), f(#, by the number of arguments. 

Remark 5.3. The matrices A M (x, <?), A M (a;, are Hermitian, so that (I5.7ab and (I5.7bb 
are also symmetric hyperbolic. They are affine linear (over W) functions of the field \P. 
The linear (over R) transformation f (a;, d x ^) depends affine linearly (over R) on 
\P(Bd x )P. On the other hand, f (x, \P) is a quadratic polynomial in the components of 
\P without constant term. There is no constant term, because M. a .% is a vacuum field. 
By direct inspection, neither derivatives of e a a nor derivatives of \ a ,% appear in the 
term A^($)-^M a ,<&- See, (23) and Q . 

Definition 5.3. Lef foe defined as in equation ( 14.4b . 

• V is a generic symbol for a quadratic polynomial in the components of the fields W 
and \P without constant term, whose coefficients are (complex) polynomials in —, 21, 

S, e ai 21, X a ,%, A a ,2l- 

• is a generic symbol for a polynomial in the components of the fields W and \P and 
all their first order coordinate derivatives, whose coefficients are (complex) polyno- 
mials in —, 21, 5*, e a ,%, ^a,a> A a ,2i, and all their first order coordinate derivatives. 

We use the same symbols V and for a vector or matrix all of whose entries are 
polynomials of this kind. 
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Remark 5.4. The vector fields D, N and L corresponding to <P = M. a ,% - 



-M 



D = -fe a , a | 



N 



2 Q Q d 



L = — + X f 

du ii 2 " " fin 



W are: 



(5.9) 



Here -2- = M-^- 4- i-£_) 
Below, e = e 0) a and A = A a a- 

Proposition 5.2. Lef = (f,LJ,z). The system d5.7al > to£es f/ze relevant / irrelevant 
form: 

\ 



L 



N 



(h\ 




1 


BU>i 


h 






—i euji 


h 
h 






2 e $l(oj4 — uJ 3 — 
—2 e 9(^4 — CJ3 — 


CJl 
















UJ 3 






— z 2 — Awi 


UJ4 






— AoJi 






\z% + 2i 9 (A (uj4 — cj,5 


ffs\ 






/2(/ 3 + Bw 8 )\ 


^5 


-I 




^6 




S 


\ujrj 






V o / 



(5.10a) 



(5.10b) 



±£> iV 



\ o o 



u 





0\ 

rL ±D 

±£> 

_U Z U 

±D N+±rL D 

U p, 

D Lj 



fzi\ 
W 



\ 




4Az 5 

\2l 2 z,5 — 2A;?4 — 30J7Z3/ 



(5.10c) 



Proof. By direct (machine) calculation. □ 

Proposition 5.3. Lef 'f' 5 ' = (s,p, y), see ( I5.6bb . ant/ reca/7 Definition ^. 4\ Then, 



'«1 

P2 

P6 
\P9j 



■pi 



= 7? tt 



_ <pt — 




-uN 



= -uN 



. pi 




—L(lu 7 ) — UJi 
L(u> 6 ) 

u D{ujj) — u D(loq) + 0J4 — uj 3 — AXui 



i-ptt 



(5.11a) 



(5.11b) 



(5.11c) 



21) 









= = ( 







uD-AX L 
lo-j u D — 2A 
2w 7 




(5. lid) 



Proof. By direct (machine) calculation. □ 

Remark 5.5. Every generic symbol that appears in ( 15.1 lb . has no constant term as a 
polynomial in the components of IP, & and their first coordinate derivatives. There is no 
constant term, because M a ,w. is a vacuum field. 

Proposition 5.4. Let W$ = (s,p, y), see ( |5.6bt . The dual system d5.7b| > takes the rele- 
vant/ irrelevant form: 



/8!\ 




S2 




Si 




S5 




Pi 




P2 




P3 




P4 




P7 




W 



( 



e (pi 



M.P6 

-MPs 






- Ps) + e (p 6 
p 5 ) - ie(p 6 

Vi 




- Pz) ~ MP4 ■ 

- Pi) + Kpa ■ 



\ 



Pz) 



P 5 
P 5 )J 



U 



(5.12a) 



u 2 




u 2 



(5.12b) 



(5.12c) 



Above, the symbols T>* are linear over R generic transformations, in the sense of Defi- 
nition\5~3\ 



Proof. By direct (machine) calculation. □ 

Remark 5.6. The overall factors u E and u Et appear in ( I5.7al i and (I5.7bb . so that these 
systems are line by line (up to a permutation of the lines) equivalent to their rele- 
vant/irrelevant counterparts in Propositions 15.21 and 15 .41 



6. Formal Solutions 

In this section we consider formal power series 

oo 

on Strip^ C M 4 , see ( 14. Il l, where for each k > 0, the coefficient function \P(k) = 
&(k){£,u) is a smooth field on R 2 x (0, oo) taking values in 1Z. By Proposition [53J 
the associated formal constraint field [&'] is itself a formal power series in —, that is, 
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Remark 6.1. It also follows from Proposition 15 . 3 1 that, for each k > 0, the coefficient 
function g^(fc) depends only on \P(£), < £ < k. 

We shall construct all formal power series solution [ ] on Strip,^ C M 4 to (SHS) , more 
precisely to the system ( I5.7al ). with initial conditions ( I5.3ab and ( 15 .3bb . We will use 
d5.7bb to show that the associated formal constraint field ] vanishes. 



Definition 6.1. Regard the components of ^(k) and^P(k), k > 0, and their formal first 
coordinate derivatives, as an infinite family of independent abstract variables. 
Set Vq = 0. The generic symbol V k , k > 1, is an arbitrary polynomial in the com- 
ponents of &(£) and < £ < k — 1, and all their first coordinate derivatives 
( g| jr !f r (£) and ■^ : 'W{£), p = 1,2, 3), whose coefficients are (complex) polynomials in 
21, u, e a ,%, \a t w, A ,2i, and all their first coordinate derivatives. It is further required 
that the polynomial Vk have no constant term, that is, Vk vanishes when &(£) and 
g|jr^(^) vanish for all < £ < k — 1 and p, = 1,2,3. We use the same symbol Vk for 
a vector or matrix all whose entries are polynomials of this kind. 

Proposition 6.1. Substitute [ M. a ,% ) (see, Definition ^. 2D for M.„ m and [W\ (see, d6.ll )) 
for & in (I5.7ab . Then is a formal power series solution to ( I5.7ab if and only if its 
coefficients i'(fc), k > 0, satisfy a system of the form 

zi(k)=V k k>0 (6.2a) 

z 2 (k)=V k k>0 (6.2b) 

z S (k)=V k k>0 (6.2c) 

£z 5 (k)=V k k>0 (6.2d) 

(l-S k0 )z 4 (k)=- 1 ^(e-§ i + 2X)z 5 (k)+V k k>0 (6.2e) 

&wi{k)=-zi(k)+V k k>0 (6.2f) 



i» 2 (k)=-(2-6 k0 )$t(tJ 1 (0)iJ 1 (k))+Vk k>0 (6.2g) 

^u> s (k) = -z 2 {k)-\u>x(k)+V k k>0 (6.2h) 



- u>i(k) — -Awi(fc) + Vk k>0 (6.2i) 



u 5 (k) = - 1 ± I (u i (k)-uj 3 (k))+Vk fc>0 (6.2j) 

u> 6 {k)=Vk k>0 (6.2k) 

u 7 {k) = V k k > (6.21) 



£ W8(ft) = z 3 (k) + 2i% (A(w 4 (A) - w 3 (At) - w 6 (*))) + ^ fc > (6.2m) 

^/i(fc)= ewi(A;)+n fc>0 (6.2n) 

^/a(fc) = -»ewi(*)+n fc>0 (6.2o) 

fs(k) =-j^$tu> s (k)+V k k>0 (6.2p) 



-§gh{k)= 2e%t(uj A (k)~Lo z (k)-Lo b (k))+Vk fc>0 (6.2q) 

^/ 5 (fc) = -2e9(w 4 (*)-^^(fc)-W6(A!)) +n- fc>0 (6.2r) 

//ere, e = e 0) a, A = A a ,2i. /« < |6.2gb , ( I6.2kl i, ( |6.2p i, (6.2qi, ( I6.2rb . f/ze generic symbol 
Vk is real valued when <P(£)(t;, u) £ IZfor all < £ < k. 
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Proof. Substitute the formal series ( 16. U into the relevant / irrelevant form of system 
(15.7al > given in Proposition l5.2l Collect all coefficients of common powers of i. □ 

Lemma 6.1. For aZZ a, 21 ^ ami a/Z DATA(£, u) as in Definition \5.1\ there is a unique 
formal power series solution d6.U to J5.7at satisfying the initial conditions J5.3at and 
( I5.3bb in the sense of formal power series. 

Proof. The initial condition ( I5.3ab forces the zeroth coefficient function i'(O) of ( 16. U 
to be given exactly as in Definition l5.ll Observe that *?(0) satisfies the fc = equations 
in ( 16.21 ). The coefficient functions <P(k), fc > 1, are constructed by induction. For each 
step fc, equations ( I6.2al ) to ( I6.2rl ) are solved exactly in this order to obtain •P'(fc). The 
right hand side is explicitly known by induction and the "upper triangular" structure 
of ( I6.2al ) to ( I6.2rl ). Whenever jj^ appears on the left hand side, it is inverted using 

the operator d^ 1 of Definition 15. II because the constant of integration is zero by the 
initial condition d5.3bt . By induction, one also verifies that W{k), fc > 0, vanishes 
when u < u , so that d5 .3bb is satisfied at all orders. It is essential at precisely this 
point that the generic polynomial Vk in Definition ^, ll has no constant term. Finally, by 
Proposition ^. II there exists a formal power series solutions satisfying the hypothesis of 
the lemma. The construction given here is forced at every step, and therefore generates 
a unique formal power series. □ 

Proposition 6.2. For all a, 21 ^ 0, all u > 0, and all DATA(£,w) as in Definition 
15.71 there is a unique formal power series [ W ] on Strip^, which satisfies d5.7ab and 
[ ] =0, and which in addition satisfies (I5.3bt and u)\ (0) = DATA. 
Moreover, for all k > 0, the value o/!f r (fc) at (£, u) e R 2 x (0, oo ) depends only on the 
restriction o/DATA(£, u) and its derivatives of all orders to the half-open line segment 
{£} x (0,u] (formal finite speed of propagation). 

Proof. Suppose such a formal power series (16.11 ) exists. The condition wi(0) = DATA, 
the fc = equations in (16.2b and the initial condition ( I5.3bb taken together give the 
formulas in Definition 15. II for all the components of ^(0), apart from 077(0), 2:2(0), 
23(0), 24(0) and U!e(0). The remaining five formulas are exactly what is required for 
the vanishing of ps(0), j/i(0), 1/2(0), 2/3(0) and pe(0), see ( 15. 1 lcb and (15. 1 ldb . Here, 
!?»(0) = (s(0),p(0), 2/(0)). Therefore, this solution must coincide with the unique so- 
lution of Lemma 16.11 and gives the uniqueness statement of the proposition. Formal 
finite speed of propagation follows from an examination of the construction of [ \P ] in 
the proof of Lemma loTTI 

To prove existence, we only have to show that [ ] = 0. Note that 

• [ ] is a formal power series solution to the linear homogeneous system ( I5.7bl ). 

• [ipt] = when u < u . 

• ^"(0) = 0oni 2 x (0,oo). 

The first bullet follows from Proposition 12.31 because [&] is a formal power series 
solution to ( I5.7ab . The second follows from d5.3bb . which implies [ <P } = [ M. a ,% ] when 
u < Uq, and from [ M.^ a ] = 0. For the third, note that ^5(0), j/i(0), 1/2(0), 2/3(0) and 
Pq{0) all vanish on R 2 x (0, 00) by the prescription of ^(0) in Definition 15. II By the 
first two bullets and by equation ( I5.12al ). we conclude, step by step, that si(0), S2(0), 
Pi(0), P2(0), ^3(0), P4(0), S4(0), S5(0), ^7(0), p&{Q) also vanish. The first equation in 
( 15. 12bb gives s 3 (0) = 0. It remains to show thatpg(O) = OonR 2 x (0, 00). By ( 15. 1 lcb , 

pg(0) = -2(e ■§. + 2 A) oj 7 {0) + 2e 4o; 6 (0) + w 4 (0) - Zirfti). 
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DefinitionOnow implies (^) 2 p 9 (0) = 0. Therefore, p 9 (0) = on R 2 x (0, oo). 

The three bullet statements imply, by induction on k > 1, that = on the set 

R 2 x (0, oo). In fact, at each step k, one verifies, in the given order, that yi{k), y2(k), 
ys(k) all vanish by ( 15. 12cb . P\{k), pz(k), Pz{k), Pi{k) all vanish by ( 15. 12ab . p§{k), 
Ps{k) both vanish by (|5.12bt , pr(k), ps(k) both vanish by J5.12ct . pg(k), S3(fc) both 
vanish by d5.12bt , and Si(fc), sa(fc), S4(k), ss(fc) all vanish by (15 . 12al . □ 

Proposition 6.3. For all k,R>0, all < |2l| < \a\ < 1, and all DATA, 

||#(fc)||c*(Q) < Pfc,fl(l|DATA|| cB+2 * + 3 (c) ) Q = £) 4 | ft |(0) x (0,2) 

where [W] is the corresponding formal solution in Proposition ^. 21 and pk t R : R — * R, 
is an infinite family, indexed by k, R > 0, of universal polynomials without constant 
term. Here, D r (0) is the open disk of radius r > in f/ze (£ , £ 2 )-plane. 

Remark 6.2. The uniformity of the estimate in a, 21, when < |2l| < |a| < 1, will be 
exploited later. In particular, it is compatible with taking the limit a = 21 J, 0. 

Proof. Observe that: 

• ||e a ,2i||c«(Q) < ¥ and ll^o,allc«(s) ^ ¥ fora11 ^ > 0- 

• ll^u 1 5llc R (Q) < 2 |j.g|jc R (Q) for all i? > and all functions g(£,u) on Q. 

The existence of polynomials po,ii> -R > 0, follows by direct inspection of & a gi(0), see 
Definition 15. II The existence of polynomials Pk.n., i? > 0, is shown by induction over 
k > 0. At each step k > 1, we use d6.2| i. By the inductive hypothesis and Definition 
16. II there is a polynomial p' k R (depending only on k and R) so that each generic term 

Vk on the right hand sides of ( 16.2b satisfies H'PfeHc^Q) < Pfc.fl(ll DATA llc J? + 2,! + 2 (Q))- 
We can assume that p' k R has no constant term, because Vk does not have one (see, 
Definition ^. U . Now, the existence of Pk,R, R > follows directly from estimating 
the non generic terms on the right hand sides of ( I6.2ab to ( I6.2rl ). exploiting the upper 
triangular structure. Only in one equation, (I6.2el i, a coordinate derivative appears. □ 

Remark 6.3. Fix DATA and let [ W a .% ] be the formal power series solution in Proposition 
16.21 The indices have been added to make the dependence on a, 21 7^ explicit. One 
can show, by induction, that *^a,a(fc)(£, u), k > 0, are polynomials in 21. Just follow 
the construction of [#a a ] given in the proof of Lemma loTTl and use the observation 
that ea.a an d -^a,a are polynomials in 21. 

Let *p be the parity field symmetry, see Remark |2T9l Then • [ !?a,2l ] = [ S'-a.-a ]• 
This is a direct consequence of <}} ■ .Ma, a = A^-a,-a> the uniqueness statement in 
Proposition ^. 2l and the fact that cji(0) is ^-even. Therefore, the ^-even CP-odd) com- 
ponents of !^a,a(fc)' ^ ^ 0' are even (odd) polynomials in 21. 

Let [i'o.o] = lima j. [^a,a] (the limit is taken coefficient by coefficient). We have 
[^0.0] = [A^o.o] + u - [<Pb,o]- The ^-odd components all vanish. By inspection, 
the e4, es components also vanish (see, Remark |2~9| ). Therefore, [#0.0] satisfies the 
hypothesis of Proposition ^. 41 in the sense of formal power series. The field [<?o,o ] is a 
formal solution to (subSHS) , and [$o ]=0. 

We now match constructions between two stereographic charts. 

Proposition 6.4. Choose a, 21 7^ 0. Pick DATA CT (£, u) as in Definition ^. 1\ for a = — , + 
and let [ W 7 ] be the associated solution in Proposition \6.2\ The following statements 
are equivalent: 
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• ^DATA^f e, «) = -|^r DATA-^f |, u) w/ien £ £ 0. 

• Flip^ • [<2> CT ] = [<?- a ] wheni^ O.Here, [V ] = [M aM ] + u- M [& cr }. 

• Flip^ • [$ a ] = [!f- cr ] w/zen ^0. 

//ere, — a = + vvrtew a = — , and conversely, —a = — when a = +. 

Proof. The equivalence of the last two bullets follows from Proposition 14.11 (b), and 
the fact that Flip o_ commutes with multiplication by u~ M . Each of the last two bullets 
implies the first. Just look at how Flip .a. acts on the component wi. The first bullet im- 
plies the last two, because Flip a. is a field symmetry, and by uniqueness in Proposition 
l6.2| (more precisely, by formal finite speed of propagation). □ 

It is convenient (see Subsection l8.2b to make the 
Definition 6.2. For all (£, u) £ R 2 x (0, oo) with £ ^ 0, set 

(Fli Pf -DATa)(£,u) = pDATA(^|,u) 

7. Energy Estimates 

In this section, we prove an abstract local existence theorem for a general class of quasi- 
linear symmetric hyperbolic systems, with a concrete breakdown criterion. Then, we 
develop appropriate energy estimates. These tools are applied in Section|8] 

Convention 7.1. In this section, 

x = (x\ x 2 , x 3 , x 4 ) = (£\ e, u, u) 

q=(q°,q\q 2 , q 3 ) = (t, £\ £ 2 , 
t = u + u 

Let D r (£) C R 2 be the open disk of radius r > around £ and, generally, B r (p) C R N 
be the open ball of radius r > around the point p. 

For any parameter vector a = (ai, . . . , afc) £ (R+) fe , the notation X < a Y signifies 
that X < CY for a constant C = C(a) > that depends only on a. Dropping the 
subscript a, the notation X < Y means Y < CY for a universal constant C > 0. 

7.1. Sobolev inequality. 

Lemma 7.1. Let 9 q = (g|r, g^j, ^f). /f 6 € (0, 2], then for all C 2 -functions /(q) = 
/0?\ 9 2 , 9 3 ) °« f « e cylinder CYL = Di/ 4 (0) x (0, 6) C M 3 which vanish for q 3 < 1/4, 

sup |/(q)| < ( £ ll^/lli=(CYL)) 1/2 , «e«0- 
qSCYL |a|=2 
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Proof. Let B = L>i/ 4 (0) x {0} be the base of CYL and S 2 the unit sphere in M 3 . For 
each q £ CYL, let T q C S 2 be the set of all quotients £ = jE^Sy where p £ B. Set 
KC) = IP-ql- We have 



|/(q)| < jj^j- f dA s2 (C)F(C) 



-qls 2 Jr a cs 2 

where 



F«)=f drr \{C H(/)(q + rC)C 



since, Taylor's theorem and the support properties of / imply |/(q)| < F (£) for all 
C 6 r q . Here, if (/) is the Hessian of /. Let C q C M 3 be the convex hull of B U {q}. 
By the Schwarz inequality, 

i / r r l (® \ i/2 / r \ 1/2 

< irTT ( / J As2{c) ldr ) (/ d 3 y| J ff(/)(y)l 2 ) 

where, r 2 has disappeared into the measure d 3 y and \M\ = (tr M T M) 1 ' 2 is the Eu- 
clidean matrix norm. Also, observe that 1(C) < 3 and |iq|g2 is bounded below by a 
universal constant, for instance 7r/100. By construction, (7 q C CYL, and the proof is 
finished. □ 

Lemma 7.2. Let b £ [1,2]. Then for all C 2 -functions /(q) = /(g 1 , q 2 , q 3 ) on the 
cylinder CYL = £>i/ 4 (0) x (0,6) C M 3 , 



1 /2 

sup, l/(q)l < ( E K/H Wyd) , «e4 (7.D 

\a\<2 



qGCYL 



Proof. By reflection symmetry, it suffices to show ( 17.1b when q £ CYL satisfies q 3 > 
I > i. Fix a smooth transition function ?/> = ip(q 3 ) : M — ^ [0, 1] equal to on (—00, j] 
and equal to 1 on [i 00). Then, 

Lemma im 

l/(q)| 2 = l(V^/)(q)| 2 < E IK(^/)lli»(CYL) 

|o|=2 

£ !l^llc 2 ([o,i]) E IK/IIWyd □ 

|o|<2 



7.2. Finite speed of propagation for a general class of symmetric hyperbolic systems. 
We show finite speed of propagation in the context of the following hypotheses: 

(FSO) U C M 4 is open and A = U x B r (Q) where B r (0) cl p ,0<r< +00, P £ N. 
(FS1) M M (g, 0), = 0, 1, 2, 3, is a symmetric PxP matrix on A. Moreover, h(q, 0) 

is an M. p valued function on A. Both M M and h are C 1 on ^4 and extend, with their 

derivatives, continuously to A, and M° > on A. 
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(FS2) 0i, 02 are ^-functions on U with values in B r (0) C M p , that are solutions to 
the symmetric hyperbolic system 

M(q, 0)0 = h(q, 0), M = M^. (7.2) 

Both 0i, 02 extend, with their derivatives, continuously to U. 
(FS3) U is either the set Cone or the intersection Cone fl(lx Half). 

Above, all the quantities are real. Cone is any set of the form 

{(i,q)G]R 4 : |q-q | H 3 < v\h - t\, t e (i 0) *i)} (7.3) 

where v > 0, qo € M 3 and icb *i £K with <o < ti are arbitrary, and Half is any open 
half-space in R 3 . We refer to v as the velocity of the set Cone. 

Lemma 7.3. Suppose (FSO), (FS1), (FS2), then the difference T = 2 — 0i satisfies the 
linear homogeneous symmetric hyperbolic system 

M(q)T = H(q)T (7.4) 

where M(g) = M(g, 0\(q)) and H{q) is a square matrix. Moreover, M(g>), its first 
derivatives and H(q) are continuous on U and extend continuously to U. 

Proof. Adding and subtracting, 

M(g, 0i)T = - (M(g, 2 ) - M(q, 0i))0 2 + h(q, 2 ) - h(q, 9i). 
Set S = (1 - s)0j + .S02, and 

<M«> - i:(jf ^ 2 ^ fi («.*w))0w + i^^c. *<«»■ 

The proposition follows from the fundamental theorem of calculus. □ 

Suppose (FS3). In this case, let 

S = (dU)n ((to,*i) x M 3 ) 
B = (SW) n ({t } x M 3 ) 

be the "lateral boundary" and "base" of U. Note that S is a piecewise smooth hypersur- 
face in M 4 . Let = 0^ dq^ be a smooth 1-form on the smooth components of S, such 
that 9(X) > for every vector X pointing out of U. 

Proposition 7.1. Given (FSO), (FS1), (FS2) a«c/ (FS3), the difference T = 2 -0i vanishes 
identically on U when 

T\ B = 0, (7.5a) 
# M ]VP(g, 0i(g)) > 0, along the smooth components of S . (7.5b) 
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Proof. We abbreviate M = M(q, 9i(q)). 
Define the C 1 "energy current" vector field 
= r T M^T on U. By hypothesis, p 
and its derivatives extend continuously to 
U. For t G [t , ti], let £(£) be the integral 
of the component ?'° overt/ n ({t} x R 3 ). 
The Euclidean divergence theorem gives 

m= [ d,r- f o» 

where C = (to,t) x R 3 , because, by ( |7.5ab , £7(to) = 0. For U = Cone, see the 
nearby figure. By d7.5bl ). the outward flux j SnC (j, v) of j through S n C is positive. By 
construction, = T T KT, where K = 9 f ,M' 1 + iJ T + H, since M' 1 is symmetric, 
and T is a solution to j7Aj . By Proposition 17.31 K is continuous on W and therefore 
bounded. Also, there is a constant k > such that |T T i-£TT| < kj°, because M° is 
strictly uniformly positive definite on the compact set U. It follows that 

E(t) < [ \T T KT\ < k j dsE{s), 
June Jt 

for all t € [tojti]- Consequently, E(t) = on the interval [to,ti]- In other words, 
T = 0. □ 

Proposition 7.2. Let U = (to, h) x X, -where X C M 3 is open, and to,h € R w/f/i 
£ < ti. Suppose (FSO), (FS1), (FS2) and fef T = 6> 2 - @i. Further suppose 

(i) T\ {to}xX = 0. 

f//J M° > a ana? —6 < M 1 < b, i = 1,2, 3, o« A, for constants a,b > 0. 
Part 1. Sef v* = ^3 (6/a) > 0. T/zen T vanishes at (t, q) eU if 

dist R 3(R 3 \A-, q) > u* |ti — 1 |- (7-6) 

Part 2. If, in addition, M 3 > on A then Y vanishes at (t, q) €lA if 

dist R 3 (Halfq n (M 3 \ X), q) > v* \ti - t \ (7.7) 

where Half q = {y £ M 3 | y 3 < q 3 }. 

Proof (of Part 1). Let q £ X satisfy ( |7.6t . Observe that, for the set Cone C U with 
velocity v* , base at time to and vertex at (ti, q), we have 

M M"| SxBp(O )>O (7.8) 

by the choice of v* and (ii). Here, the lateral boundary S and the 1-form 9 are just as in 
the discussion above Proposition 17. 11 We can now apply Proposition ^. l| to Cone, and 
conclude T\^^ = 0. □ 

Proof (of Part 2). Suppose that M 3 > on A. Let q e X satisfy d7J] i. The boundary 
of the set Cone n Half q c U, where Cone has velocity v* , base at time to and vertex 
at (ti, q), has two smooth components. On the "round" one, the inequality ( 17.81 ) holds 
again by the choice of v*, and on the "fiat" one by M M M = M 3 > for 8 proportional 
to (0, 0, 1, 0). We have r| ConenHa|f = 0, by Proposition|7j] □ 



(ti, qo) 




28 



7.3. Existence / breakdown theorem. 

Assumptions for the Existence /breakdown theorem: All the quantities here are real. 

(EBO) Let U = (-00, T) x R 3 , where Tel and letA = Ux B 2 (0) C R 4 x R p . 
(EB1) ]VP(g, 0) is a symmetric P x P matrix on A. Particularly, M° > |. 
(EB2) /i(g, 0) is an R p valued function on A. 

(EB3) Both h and M M are smooth on A and their derivatives of all orders extend con- 
tinuously to A. 

(EB4) JC C Q C R 3 with JC compact, Q open, such that on (— oo, T)x (R 3 \/C) x _B 2 (0), 
the matrix M' 1 is constant and denoted by P M and h(q, 0) = fl(t) 0, where fl(t) 
is a matrix depending only on t = q°. It is assumed, ]X1 > \. We can naturally 
extend JVP and /i, by W and ^(i)6>, to (-oo, T) x (R 3 \ /C) x R p . 

We now formulate and prove an existence theorem for the quasilinear symmetric hy- 
perbolic system 

M(q,0)0 = h(q,0) M = M^. (7.9) 



Proposition 7.3. Suppose (EBO), (EB1), (EB2), (EB3), (EB4). 

Part 1. For each to < T, there is a t± G (to, T] and a smooth solution : [to, ii) x 
R 3 — ► R p of (17. 91 > w/f/z trivial initial data, 0(to, ■ ) = 0, smc/i f/zaf supp (9 C [io, ^i) x 
B r (Q) for some finite r > 0, a«c/ 

%ii)xS)cBi(0)cl P (7.10) 

anrf 5«c/i f/iaf t\ ^ T implies either one or both of: 

(Break) x : ©([t ,*i) X Q) ^ 1 (0)cM P . 

(Break)2-' 77ie vector field d q is unbounded on [to, t±) x Q. 

Part 2. Suppose in addition that M 3 > on A ant/ ft.(g, 0) = when q 3 < \. Then the 
solution of Part 1 vanishes identically for q s < i. 

Proof (of Part 1). Fix a smooth transition function ip = ip(\0\) : K — » [0, 1] which is 
equal to 1 on (— oo, |) and equal to on (|, oo). It is for this reason that B 2 (0) appears 
in (EBO). Set 

N = V'M+(1-V)^, g = 4>h+(l-ip)]I9. 

By construction, g and the symmetric matrix N p are smooth on B = U x R p , and their 
derivatives of all orders extend continuously to B. Note that N° > | onB and there is 
a constant b > such that -b < N' 1 < b, i = 1, 2, 3 on [t , T] x R 3 x R p . The latter 
statement follows from the fact that N p is constant (= ty/V 1 ) on the complement, in 
[to, T] x R 3 x R p , of the compact set [t , T] x JC x B 2 (0). Fix the velocity u* = 2V3 6 
(see Proposition ^. 2b . 

We want to reduce our existence /breakdown theorem to ]Tay| . To do this, fix L > 
big enough so that 

JC C Cube = 3 
dist R 3(c»Cube, JC) > l + v*\T-t \ 
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and smoothly extend N and g from (— oo, T) x Cube x R p to spatially periodic matrix 
and vector valued functions on (— oo,T) x R 3 x R p . With these preliminaries, the 
hypotheses of Proposition 2.1 on page 370, Proposition 1.5 on page 365 and Corollary 
1.6 on page 366 in |Tay) are all satisfied, and there is a r e [to,T] and a spatially 
periodic smooth solution : [to, t) x R 3 — > R p with trivial initial data at t = to to 
the symmetric hyperbolic system corresponding to the spatially periodic extension of 
N and g, such that, if r ^ T, then the vector 

(0, d q 0) G R P 8 (R 4 ® K p ) 

is unbounded on [to, r) x R 3 . (There is one caveat: ]Tay| , for convenience, considers 
systems defined for all time. By direct inspection, his argument applies to any open 
subinterval of R.) 




K.+ (2L,0,0) 



OCube + (2LZ) 



Let J = K, + (2LZ) 3 . By construction, the spatially periodic system introduced 
in the last paragraph reduces to p<9 = ft{t)0, on (-oo, T) x (R 3 \J)x R p , and 
admits the trivial solution. Intuitively, "signals can travel at most a distance v*\T — to |", 
which is less than the distance between JC and d Cube. This intuition is formalized 
by applying Proposition 17.21 to the open set (to,T — e) x (R 3 \ J) for arbitrarily 
small e > 0. Consequently, vanishes at every point (i, q) € [to,i~) x R 3 with 
distK3(^, q) > | + v*\T — i |, because 0\t=t o = 0- It follows from our choice 
of L that the periodic solution vanishes in a neighborhood of [to, r) x (d Cube). For 
this reason, and because of (EB4), the modified field 



[*o,t) 



9 q 



0, ifqe (R 3 \Cube) 
0(q), ifqe Cube 



(7.11) 



which we continue to call 0, is a smooth solution to the non-periodic system N6> = g 
with trivial initial data. Moreover, if r ^ T, then the vector (0, d q 0) is unbounded. 

Suppose r^T.We show that (0, d q 0) is bounded on V = [t , r) x (R 3 \ Q), and 
consequently, unbounded on [to, t) x Q. For any time t 2 G (to, t), decompose 

[to,T) xR 3 = V1UV2UV3 

Vi = [to,r) x (R 3 \Cube) 
V 2 = [t ,t 2 ] x Cube 
V 3 = (t 2 ,r) x Cube 



By ( 17.1 U . the vector (0, d q 0) is bounded on V n Vi and, by compactness, also on 

vnv 2 . 
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Let 

d = dist R3 (/C, K 3 \ Q) > 0. 

To verify that (0, d q 0) is bounded on VDV3, we choose £2 £ (£o, r) with v* |t — £2 1 < 
d/2. Intuitively, this choice means that "signals can travel at most a distance d/2 in the 
time interval (£2,1")", which is less than the distance between JC and R 3 \ Q. There 
exists (see, for instance, ]Tay2[ ) a unique smooth solution 2 : [£2,X) x R 3 — > R p 
to the linear system tyl0 2 = fl(t)0 2 with initial condition (0 — 02)\t=t 2 = 0. In 
particular, (0 2 , d q 02) is bounded on the compact set VflV 3 . By (EB4), is a solution 
to the same system on [£o,t) x (R 3 \ JC). As in the last paragraph, Proposition I7.2I 
implies that = 2 on [t 2 , r) x (R 3 \ Q) D V H V3. Consequently, (6>, <9 g 6>) is 
bounded on V n V3, and we are done. 

The final step is to remove the transition function ip. Let 



X = 



he[t ,T) 0([to,t] x Q) c Si(0) c R p }. 



We show that X = [£o,£i), where £1 £ (£o,t]. First, to £ I since the initial data 
vanishes. Second, if t' e X, then [to, f] G X- Let = supX. If = t, then 7^ X. 
If £1 < t, assume by contradiction £1 £ X. Then, the compact set <9([£o,£i] x Q) 
is contained in a ball £? r (0) C R p of radius r < 1. However, 9t<9 is bounded on 
[£o, £1] x R 3 , since supp R3 (dt0)(t, ■) C Cube is compact for all £ G [£o, £i]- Therefore, 
ti + e E X for all sufficiently small e > 0. 

The smooth solution of ( 17.9b that we are looking for is 0\[t o .t 1 )xR 3 - Indeed, it has 
trivial initial data, support contained in [£o,£i) x Cube and satisfies (17.101 i. If q £ 
[£o, £1) x Q, then \0(q)\ < 1 and ip(\0(q)\) = 1, by the definition of t\. In this case, 
the system N6> = g reduces to ( PT9l i. On [£ , £1) x (R 3 \ JC), (EB4) directly implies that 
the system also reduces to (|7.91 l. 

If £1 7^ T, there are two alternatives: £1 < r < T and t\ =t <T. For the first, we 
use the continuity of on [to, £1] x Q to conclude that 



0(Mi) X C) = &([to,h] x Q) £ Bi(0) 

since £ x ^ X. That is, we have (Break)!. For the second, r 7^ X, and (0, 0) is 
unbounded on [£o,£i) x Q. Since is bounded, (Break)2 applies. The proof of Part 1 is 
complete. □ 

Proof (of Part 2). Let Half = {q £ R 3 | q 3 < |}. The assumption h(q, 0) = 0, when 
q 3 < i, implies that Q\ = is a solution to N0 = g on (£ , £1) x Half. Also, 2 = 
is a smooth solution, and (0 2 — <9i)|{t }xHaif = 0. The assumption M 3 > on A 
implies ])d 3 > and consequently, N 3 > on U x R p . At last, Part 2 of Proposition 
17.21 applied to the open set (£o,£i — e) x Half with arbitrarily small e > 0, forces 
2 — 0i = to vanish on [£ , £1) x Half. □ 



7.4. Energy Estimate. 

Assumptions for the energy estimate: All the quantities here are real. 

(E0) U = X x 0(b) where 0(b) = R 2 x (0,6) C R 3 , b £ [1,2] and X = (£ ,£*), 
-00 < £ < £* < -1. 
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(El) M M (g) is a symmetric PxP matrix on U. Particularly, for all g G U, \ < M° < 2 
and M 3 > 0. We assume the integer P is less than some big absolute constant, say 
P < 10 9 . 

(E2) H (g) is a P x P matrix on 

(E3) Src(g) is an R p valued function on U. 

(E4) (9(g) is an R p valued function on U, which is a solution to the linear, inhomoge- 
neous, symmetric hyperbolic system 

M(g)0 = PT(g)0 + Src(g), M = M^. (7.12) 

(E5) Fix a non-negative integer R. Then, M M , (resp. H, Src) are C R+1 (resp. C p ) 
functions on U, all of whose derivatives of order < R + 1 (resp. < P) extend 
continuously to 

(E6) supp 0(t, ■ ) and supp Src(i, • ) are contained in a ball in R 3 with radius indepen- 
dent of t. Moreover, and Src vanish identically when g 3 < i. 

Let M p = M Pl © R Pa © R Pa . We decompose 

= (0i, 2 , 6> 3 ), Src = (Srci, Src 2 , Src 3 ). 

Each PxP matrix is decomposed into nine blocks of size P m x P n , where m, n = 
1, 2, 3. Especially, 

M A ' = (M^ n ) mi „ =li2 ,3, H = (P m ,i)m,7i=l,2,3- (7-13) 

(E7) The matrix ]VP is block-diagonal, ]VP = diag(Mf, M$, M%), and (M 2 )ij = 

/•*(?) foj (afr + ^r)> h 3 = 1> • • • > -P2. for some function 
(E8) P M are constant symmetric PxP matrices, with | < P° < 2, p 1 = 0, = 

and P 3 > 0. 

(E9) fl(t) is a P x P matrix depending only on i with R Pl © R p ' 2 © R P3 block-form 

/ \ 

V o |t|-i# 2 ir 1 ^/ 

where |?2, ^3 are constant matrices and flj, < 0. 
Definition 7.1. For every open X C R 3 , f/ie energy of f contained in X at time t is 




(7.14) 



anaf f/ze supremum norm 

Sup^{/}(<) sup S up|d Q /(i,q)| (7.15) 

|a|<fc qSA" 
a6N* 

for any scalar, vector or matrix valued C -function f. As usual, we denote d a = 
n^=o(^) aM where = , for any a = (0^)^=0,1,2,3 € Nq\ The pointwise norm 
I • \ is always the Euclidean norm (for matrices, \A\ 2 = Iy{A t A)). 
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(E10) There are constants ci > and J > such that for all 

\t\ 2J+2 Eg {b) {Sr Cl }(t)) 



\t\ 2J E R {b) {Src 2 }(t) 
\t\ 2J+2 E R {Src 3 }(t) 



(Ella) Assume R > 4. There is a constant c 2 > such that for all t £ X: 

\t\ 2 E*{M» - WW 



I 4 ! 2 E 0(b){ H ln - flln}(t) 



> < ci 



J 0(b)1 

l*| 2 E§(b){ H 3n - fi3n}(t) 



(Ellb) Assume i? > 0. There is a constant C2 > such that for all iel: 



I* I Sup, 



(max{l,_R» 
0(6) 



> < C 2 . 



\t\s« 9 ^ b) {H ln -p ln }(jt) 

Su P < ^ b) {H 2n -W2n}(t) 
\t\Sup^ b) {H 3n -W3n}(t)} 

Proposition 7.4 (Energy Estimate). Suppose the hypotheses (E0) through (E10) hold, 
and, also, either (Ella) or (Ellb) holds. Let Jq > and assume J > Jo, see (E10). 
Then, there are constants C3(X) € (0,1), c^{X) > depending only on X = 
(R, Jo, \fli\, |#2|, l-^l), s«cA f/zflf c 2 < c 3 (X) a«J < c 3 (X) imply that 



Eg ib) {9}(r) < Ci (X) 



l*o| J y/Eg (b) {0}(t o ) + c a 



(7.16) 



/or all t € X (see, (EO)for the definition ofX). 
Proof. In the proof, we denote S fl = -E^Vm an d Sup( fl ) = Sup 



0(6) 

Preliminaries 1: For a function / with values in 



(R) 



1,2,3, we define the 



energy naturally associated to the linear symmetric hyperbolic system jl .\2\ 

E?{/}(*)=/ d 3 q(/ T A/°/)(t,q) , Ef {/}(<)= £ E?{fl a /}(t) 



|q|<_R 



(7.17) 



See, ( 17.13b . This energy is comparable, by (El), to the one defined in (17.14b . Namely, 



E R {f}(t) < 2 Ef {/}(*), Ef {/}(*) < 2 £*{/}(*). 



(7.18) 



If i? > 2 and / is a vector or matrix valued C R function, Lemma |J2l impries: 

Sup (fl - 2 ) {/}(*) < fl y / £TO). (7.19) 
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If a e No, \a\ < R and i? > 4 (to be used in the case (Ella)) implies 

E R {fxh}{t) < R E R {h}{t) E R {f 2 }(t) (7.20) 
E°{ [d a ,h}f 2 }(t) < R £ E R - l {d^h}{t)E R - l {f 2 }{t) (7.21) 

101=1 

whereas R > (to be used in the case (Ellb)) implies 

E R {fif 2 }(t) <r (SapW{/i}(t)) 2 ^{/ 2 }(t) (7.22) 
£°{ [9 Q ,/i]/ 2 }(t) < fl ^ (Sup^-^/iKi)) 2 ^- 1 ^}^) (7-23) 

101=1 

Inequalities ( 17.221 ) and ( 17.23b are direct consequences of the product rule. The in- 
equalities ( 17.20b and ( 17.211 ). require, in addition, R > 4 and the Sobolev inequality 
(17.19b . In fact, by the product rule, 

E R {fiM(t) <r Yl I d 3 q | ^"/x (*, <l) | 2 | q) | 2 

,_, , JO(b) 



\a\+\0\<R 



For each pair of multiindices (a, j3) with \a\ + \/3\ < R, at least one of | or| or |/3| is less 
than or equal to R — 2, say a. Then, by the Sobolev inequality, 

sup \d a h(t, q )\ 2 <E R {h}(t) 

q<EO(6) 

Inequality ( 17.20b follows at once. This argument works for R > 3. An entirely similar 
argument gives (17.21b . but with i? > 4. 

Preliminaries 2: In this subsection, t € X and a € Nq, |a| < i?, are arbitrary. We 
apply <9 Q to ( 17.12b and obtain (all the derivatives make sense classically) 

M(<9 Q <9) = fld a + (Sf , S£) + <9 a Src (7.24) 
(5?, 5?, ^) = a Q ((iJ-^)6>) + [a Q ,^]6> + [M^-^,9 Q ]a M 0. 

If > 4, d720t . (ITITT i and (E9) imply that 

E°{sn < R { e -m+ m " L m * +E e^-w* 1 }} e r {o}. 

J=l 1 ' M=0 

If i? > 0, dT22l . (177231 and (E9) imply that 

,2 , |^2| 2 + |^ 3 | 2 



E°{sn<R{Y,( s ^ {R) i H »- : P^)'' 



i=i 



E (Snp< fl >{M" - P M }) 2 } £ R {<9}. 



M=0 



34 



If (Ella) or (Ellb), it follows from the inequalities just above that 

\t\ 2 E°{snm 

E°{Sn(t) \ <r clE R {0}(t) (7.25) 
\t\ 2 E°{Si}(t)\ 

where c, = max{c 2 , \t*\-\\jl 2 \ 2 + |^ 3 | 2 ) 1/2 }- 

Estimates: We derive the energy inequalities 17.331 stated below. For i = 1,2,3, 
define "energy currents" associated to = (<9i, <9 2 , 6*3) (see (E4)) 



tf[9 i ]{q) = (eTM?e i )(q). 



(7.26) 



(Warning: we never sum over repeated "lower" indices.) The important current identity 



d^[9i] = Of (d^M?) t + 2 0f(M? d^Oi), 



(7.27) 



follows from (Mf ) T = Mf , see (El). For t el, let r_ = max{t , t - 2}. Let 

T>i(t) =2? 3 (t) = {(t,q) €U\ te (t_,t)} (7.28) 



D 2 (r) =D 1 (r)n{g|g 3 -g°<6-r}. 



For the case i = 2, see the nearby figure, 
where to < r-x < to + b < t\ < t* . E nergy 
estimates are obtained by integrating ( |7.27t over 
T>i(r) C U = T x O and applying the Euclidean 
divergence theorem. The divergence theorem gen- 
erates integrals over the boundary dT>i (r), which 
we now discuss. Recall (E6). The q° = r bound- 
ary contributes E°{6>i}(r). There is no contri- 
bution from the q 3 = 0, by (E6). For i — 1,3, 
the q° = t_ boundary contributes — E°{6>i}(r_), 
and the contribution from q 3 = b is non-negative, 
by (El). If i = 2, there is always a boundary con- 



(7.29) 



Tl -b 



tribution from q 6 — q° = b — t and it vanishes by 
(E7). If i = 2 and r < to + b, there is an addi- 
tional boundary contribution at q° = to which is 

> -E°{0 2 }(t o ). 




The discussion of the last paragraph literally transposes from (9; and \6>i] to d a 0i 
and j!f[d a &i\, for |a| < R. The current [d a 0i] is C 1 and extends, with its deriva- 
tives, continuously to U. The preceding analysis of the boundary terms gives the general 
inequalities 

^{d a e i }(r)-k i (r)EP i {ff'e i }(T-)< / d 4 q d^[d a 0,Kq). (7.30) 

fori = 1, 2, 3. Here, by definition, fci(r) = kj,{r) = 1 for all t, whereas /c 2 (r) vanishes 
when r_ > to an d is equal to 1 when r_ = to. Summing over |a| < R, 

Ef{0 t }(r)-fc t (r)Ef{e ?; }(r_)< / d 4 q ^ ^fFe.K?)- (7.31) 



\a\<R 
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The current identity (17.27b . with d a 0i in the role of Oi, and ( I7.24l i imply 

d^[d a 0i](q) (7.32) 

= 2(d«9,) T l (V^- + \ (S fl Mf)(a a i ) + ((PSrc,) + s?\. 

< 3=1 J 

For i = 1, 2, we directly estimate the right hand side of ( 17.311 ). by using Schwarz's 
inequality for the spatial part of the integral, and (E6), (E10), ( 17.251 ) and ( 17.18b . For i = 3, 
we first exploit fl 3 < (see (E9)) to drop the term 2(d a 3 ) T fl 33 (d a 3 ), and then go 
on as before. We also use the estimate I^M^I = |5 M (1VP - jwP)! < R c 2 < 
c* that holds when either (Ella) or (Ellb) is assumed (in the first case, we use 
d77T9>). Abbreviating E z = Ef {6>J and E = Ei + E 2 + E 3 , we have for all rel: 

Ei(r) - Ei(r_) < x J^^ V^W) (c* VW) + ^j) 

E 2 (r)-E 2 (t ) <x J dt v/E^^VeTW+c* VW) + j^j) (7.33) 

E 3 (r) - E 3 (r_) < x ^ ^ (Te^I) + v^(*) + ||f 

where X is defined as in the proposition. 

For each A = (Ai,A 2 ,A 3 ) € (0, oo) 3 , define 

J(A) = \tel sup |t| 2J E,(t) < Al i= 1,2,3} 

L re[t ,t] } 

Assume A satisfies (recall that J > Jq > 0, by assumption) 



-4i 


> 


|t |VEi(io) 


Ai 


>f c i 


Ai 


> § c,|A| 




A 2 


> 


2|i |VE 2 (to) 


A 2 


> 8Cci 


A 2 


> 8C(Ai 4 


-c|A|) 


A 3 


> 


|i |VE 3 (to) 


-4, 




-4, 


+ 


c\A\) 



(7.34) 



where \A\ 2 = A\ + A\ + A\ and where C — C(X ) > is the maximum of the three 
constants of proportionality in the inequalities d7.331 l. It is a direct consequence of the 
inequalities (17.331 ), ( 17.341 ) and the continuity of I 3 t n E^(t) that J {A) is an open 
and closed sub-interval of I which contains to. Therefore, J (A) = X. To see that J {A) 
is open in X, first observe that for every r e J {A), the inequalities (17.33b , ( 17.341 ) imply 
the strict inequalities E;(t) < (A;|t|~"') 2 , and then use continuity. 
For each A > 0, set 



A{\) = A (l , 1 + 8C, ! + £(! + 8C)) . 



The three rightmost inequalities in ( 17.34b are homogeneous (degree 1) in A, and hold 
for A(X), A > 0, if and only if they hold for A(l), which is the case if c* > is 
sufficiently small depending only on X, because it is true for = 0. The definition 
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of c» right after ( PT25l l implies c* < (1 + \ fi 2 \ 2 + |ftt| 2 ) 1/2 c 3 (X). Consequently, the 
condition on holds if Cs(X) is suitably small. If 

A > A = 2|to|- / v /E(t n ) + max{8, J^Cci > 

then the remaining inequalities in ( 17.34b hold for A(X), that is J(A(X)) = I. By the 
definition of J {A), we have J(A(X )) = T. By dTT8l l, inequality ( T77T6l > follows if 
C4(X) is sufficiently big. □ 

Remark 7.1. Once the system of inequalities ( 17.33b has been established, the rest of the 
argument is abstract, in the sense that it holds for any three functions E,, i = 1,2,3, 
satisfying ( 17.33b . 



7.5. Refined energy estimate. We proved a finite speed of propagation theorem for for- 
mal power series vacuum fields [&], see Proposition ^. 21 The refined energy estimate 
obtained in this subsection plays a similar role for a classical vacuum field 
To make the last statement more precise, recall that the energy estimate for the symmet- 
ric hyperbolic system ( 17.12b was obtained by integrating the divergence current identity 
( 17. 27b over appropriate open subsets C C M 4 . We now construct more refined sets C 
which allow us to estimate the energy "localized in the (S, 1 ,!; 2 ) plane". We will be 
guided by the basic requirement that the boundary integrals in the divergence theorem 
have definite signs. That is, the boundary of C must be non-timelike (with respect to the 
symmetric hyperbolic system). 

Convention 7.2. Until further notice, we use the coordinates 

x = (x 1 , x 2 , x 3 , x 4 ) = £ 2 , u, u). 

Recall the matrix differential operators A(<P) and A(<P) associated to <P = (e,j,w). 
See, ( 12.5b , (12.8b . Suppose 8 is a one-form and suppose (see, (12.2b ) 

fN DY 

M^>°> M JV">O, °Ad L) -° (7 ' 35) 
The last inequality is in the sense of Hermitian matrices. Then 

(9 (tl A"(#) > 0, M A /i (#) > 0. (7.36) 

For each xq = (^o,u ,uq) £ M 2 x (0, oo) x (— oo, 0) and choice of constants fco, ki, f > 
0, where d < u and t) < |uo| _1 , set 



where 



C= (J (£> r(3il „)(£o) x {(«,«)}) 

(u,u) £ B 

T= (J (dD r ^, u ){£o) x {(u,u)} 
(«,«) e B 



B=(0,« o -i))x(-(x),- R i I ^)cR 2 , 
r(u,u) — k + ki \uq - u\ ■ \\uq\~ 1 - \u\~ \. (7.37) 



37 



More geometrically, C is a disk bundle over the base B, and T the corresponding circle 
bundle. The set C is an open subset of K 4 . Note that, r : B — * (ko, fco + ^1 l^ol/l^ol), a 
bounded set. The set C has piecewise smooth boundary. We concentrate on the smooth 
piece T C dC here. Let 8 be a 1-form along T whose kernel coincides with the tangent 
space to T and for which 6(X) > if X is a vector pointing out of C. We choose 

2 

i=l 



1 

Kl 



du + ki 



■ du, 



e-e 

— Col 



(7.38) 



Proposition 7.5. Let x G T. If e%(x) > and the inequality 

ki > 2 max j^=VN 2 + M 2 , M 2 VM 2 + M 2 } (7.39) 
holds at x, then (17.361 holds at x with 9 given by (17.381) . 

Proo/ If M L*< > and > and det M ( g f ) M > 0, th en ( 17351 ) and therefore 

( 17361 ) hold. The condition e 3 > implies A1 I/ M > 0. By d739l . 



fV + ^es + fciM 2 \u -u\ > ^ > 



which implies 0^7^ > 0. Finally, e 3 > and ( 17.391 ) imply 

e 3 k 



^e 4 + £ 2 e 5 + fci H" a |mo -Ml)- l^d + ?e 2 | 2 > 



and therefore det M ( g £ ) M > 0. □ 



Remark 7.2. Proposition 17 . 5 1 will be applied as follows. Fix <P, and consider symmetric 
hyperbolic systems with differential operators given by A(^) or A(<P). Then, if the 
assumptions of Proposition [73] are satisfied for all points x 6 T' C T, the boundary 
integral J^, (j, v), where j is the energy current vector field, is non-negative. 

Convention 7.3. Observe that the definitions of C and T depend only on the parame- 
ters ko, ki, £), £o> Uq> u o- For the rest of this paper, we make the specific choice of 
parameters 

= 1CT 3 , k = i, ki = yg0~\ uo = -g^" 1 , u =b + d (7.40) 

For each fee [1,2] and £o G K 2 , we denote the corresponding sets by C(£o,b) and 
J-(£o, b). The base B is given by B = (0, fe) x (— oo, — and the radius function 
r(u, u) takes values in (j, h) on B. 

Recall the far field ansatz (see, Section[5|) <P = M. a ,% + u~ M 'P, where 'P = (/, ui, z). 
The ansatz depends on the scaling parameters a and 21. 

Convention 7.4. For the rest of this paper, the parameters a, 21 G R are restricted by 

< |2t| < |o| < 1CT 3 (7.41) 
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Proposition 7.6. Let d be fixed as in ( 17.40b . Let £ S R 2 anrf 6 € [1, 2], Assume a, 21 € 
K satisfy ( 17.4 U . //J ;n addition, 

W{x)\ < 5 



at x = (£,U, u) <G .F(£o) 6) ant/ |£| < 4 |§|, f/ie« f/ze assumptions of Proposition 17.51 

Proof. The first five components of ^ satisfy < |tf^| < 5. Consequently, the first 
five components of <P satisfy 

MM < \u\\p-^e aM \ + \u\\^fi\ < |e 0i2[ | + ^l/il < f |a| + 50 i = 1, 2 

e 3 = 1 + £/s > 1 - T^l/sl > 1 - 50 2 > (|) 2 
H 2 | e4 | = ^|/,|<^<5c. ^ = 4,5 

In the first line, we use |£| < 4 1 = The proposition follows by direct inspection. □ 
Convention 7.5. For the rest of this subsection, we use coordinates 

q={q°, q\ q 2 , q 3 ) = (t = u + u, £ 2 , u) 
Assumptions for the refined energy estimate, i) is defined in (17.40l >. 
(REO) I = (t ,f) where -oo < t < t* < £ G K 2 , fee [1,2], 

W =U ({*}>< 0(£o,M)) C M 4 

o(Zo,b,t)= |J (^(( )a )Ko)xM), 

«e(o,6) 

r / (t J M) = | + T^|& + 5-^|-|23-^|. 
(RE1) - (RE9) are formulated identically to (El) - (E9). 

(RE10), (RElla), (REllb) are formulated identically to (E10), (Ella), (Ellb) with the un- 
derstanding that and Sup^L are replaced by E§^ o b ^ and Sup^ o b t y see 
(T77T4b and (177131 . 

(RE12) Let the 1-form 6 be as in (17381 . Then, 6», t rvF 1 > on 

(SW) fl(lxl 2 x (0, b)) 

Remark 7.3. U is a bundle over X with fiber 0(£o, &j i) C R 3 at t e X. The fiber is an 
open disk bundle over the w-interval (0, b). An equivalent description of the fiber is 

0(£o,MH{q = (£\£ 2 ,«)elR 3 I t-u)eC(£o,b)} 
for each t 6 X. It is important that 

(SW) n (I x I 2 x (0, 6)) c .F(fo, b). 
For each t E X, the map r'(i, • ) : (0, b) — ► (4, |) is decreasing. 
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Proposition 7.7 (Refined Energy Estimate). Suppose that the refined hypotheses (REO) 
through (RE10) and (RE12) hold, and, also, either (RElla) or (REllb) holds. Let J > 
and assume J > Jo, see (RE10). Then, there are constants C3(JT) £ (0, 1), c<i(X) > 
depending only on X = (R, Jo, \f^2\, |-^3|)> such that C2 < ca(X) anrf < 
Cs(X) imply that 



for all t £ X {see, (E0)for the definition ofX). 

Remark 7.4. We use the same names for the constants in the assumptions and statements 
of both energy estimates, Propositions l7.4l and l7.7l This was done for convenience, and 
does not imply that there is any relationship between them. 

Proof. This proof completely mimics the proof of Proposition 17.41 with a few modifi- 
cations. First of all, our previous conventions that E R = E^,^ and Sup( fl ) = SupS^ 

are replaced by the conventions E R = Eq,^ q b ~ and Sup'- 72 ' 1 = Sup^^ o b t y Also, the 
definitions ( 17.171 i are replaced by 



The inequalities (|7TT8T >. (17191 . (l7T20t . (I7T21I ( 177221 d7T23l still hold wit h thes e modifi- 
cations. The only one that requires discussion is the Sobolev inequality ( 17.19) . For this 
purpose, let CYL = Di (0) x (0, b) and : CYL — > 0(£ o , &, *) be the diffeomorphism 

<K£,Ji) = (fo +4r'(i) £, Then, 



The second inequality follows from Lemma [7721 The first and third inequalities are 
direct consequences of the chain rule, because all derivatives of order up to R — 1 of the 
Jacobians of <fi and have finite sup-norms on their domains of definition depending 
only on R, especially, independent of £o, b and t. 

Observe that in (T7.281 l. (I7.29K the set U is now given as in (REO). Estimate ( 17.30b still 
holds. By construction, D\ (r) is a disk bundle over the (t, u)-rectangle (r, r_ ) X (0, b). 




(7.42) 




<* sfE^Jf^m) <R ^Eg (iobt) {f}(t). (7.43) 
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The boundary dD\{r) has five components, four of them arising as disk bundles over 
the boundary of the rectangle, the fifth is a circle bundle over the interior. The treatment 
of the first four components is unchanged. The fifth is accounted for by (RE12). The 
domains ^(r) and D^ir) are handled in the same way. 
The rest of the proof is completely unchanged. □ 



8. Classical Vacuum Fields 

In Section|6]we have constructed formal power series vacuum fields [ ^ ] by solving an 
initial value problem for d5.7ab . The goal of this section is to prove the existence of an 
actual, classical, vacuum field <F, for which [ >P } is rigorously an asymptotic expansion. 

Convention 8.1. We adopt the conventions of Section [7] We use the coordinates q = 
(t, q) (see, Convention 17. lb . Keep in mind that u = u(q) = q° — q 3 . To conveniently 
translate between the x coordinate system (Sections [2] through [6]l and the q coordinate 
system (Sections|7]and|8]l, we abuse notation and write f(q) instead of f(x(q)), for any 
function /. It is also implicit that partial derivatives are adapted to the new coordinate 
system. For example, the matrix differential operator A M (x(g), ^F) g|jr is abbreviated 
as A»(q,V) » . 



8.1. Preparatory Definitions and Estimates. The goal of this subsection is to make the 
necessary definitions and estimates so that the Existence/Breakdown Theorem and the 
Refined Energy Estimate can be applied to ( I5.7ab and ( I5.7bb - 

Convention 8.2. (I5.7al > and (I5.7bb are equivalent to real symmetric hyperbolic systems 
for 1Z = M 31 and 1Z = M 32 valued fields, respectively. See, Remark l2~8l and Proposition 
15.11 This equivalence will be implicit each time the Existence/Breakdown theorem and 
the Refined Energy Estimate are applied to (I5.7ab and ( I5.7bl ). or to equivalent systems. 

Convention 8.3. In this section, C m is a vector space over R with dimension 2m. A 
linear map from C m to C™ is, by convention, linear over E. It can be represented either 
as a 2n x 2m real matrix, or as an n x m complex matrix which may have the complex 
conjugation operator C as matrix elements. We adopt similar conventions for the real 
subspaces K c C 5 C 8 C 5 and K c C 5 C 9 C 3 . 

Convention 8.4. The notation F(q, f, d q f, . . .) displays the explicit pointwise depen- 
dence of F on q, f(q), d q f(q), . . . 

To put ( I5.7ab in the form required by Propositions 17.31 and 17 . 71 we use 

(51) a, a G K satisfy Convention|74] 

(52) [!?■] = X)^=o(m ) fe is the formal power series solution in Proposition 16.21 
corresponding to DATA(£, u) = DATA(q) which vanishes for q 3 < u Q = i. There- 
fore, [ty] vanishes when q 3 < i by Proposition ^. 21 Fix an integer K > 0, and set 
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(53) The field \P is expressed in terms of (h, a, £) by 

V = (f,u,z) = V K + (h,a,£). (8.1) 
Let S — (Ex, S 2 , S3) be the field given by 

(Si, S 2 , S 3 ) 

= (4,4,4,4,4) © (hi,h 2 ,h4,, h s ,ai, a 2 , 03, 04,08) © (/i3, 05, 06, 07) 

There is a permutation matrix 7r so that 

(/i,a,4 = 7r(Si,S 2 ,S 3 ) (8.2) 

The field S = (S x , S 2 , S3) takes values in tt" 1 ^ C C 5 ffiC 9 ffiC 4 . The permutation 
is required for Proposition l7.7l see (S6). 

(54) System (I5.7at is abbreviated as A(q, & = f(q, (see, Convention WW. Some 
of its properties are discussed in Remark l531 System ( I5.7ab is equivalent to 

B(g, S)S = Q(q,E)E + Src(g) , B = B"^ (8.3a) 
B(g,S) =7r _1 A(g,aO{:+7rS)7r (8.3b) 

Q(g,S)i7 = 7r- 1 -^- / ds'f- A(q,snII)<P K + i(q, V K + s'nE + snll) 
ds s=oJ \ 

(8.3c) 

with the source term 

Src(g) = 7T- 1 (ffoflfc)- A(q,9 K )¥ K ). 

The transformation Q(q,E) acting on ir~ 1 TZ is linear over R. Note that the bracketed 
expression in ( I8.3cb is a quadratic polynomial in s and s'. The operator ^ | s _ Q J Q ds' 
selects certain combinations of its coefficients. 

(55) The matrices B^ and Q are affine linear (over R) in S. Let B M (g) and Q(g) be the 
R linear maps given by 

B»(q)II= A| s=o B^(g, S 77) , Q( g )i7 = £| s=0 Q(g, si7) 

We have, B^(g, S) = B*%, 0) + B»(q) E. Similarly for Q. 

(56) The three by three C 5 © C 9 © C 4 block-decomposition of B is 

B = diag(Bi,S 2 ,B 3 ), B 2 = lgL, B 3 = 1 4 N, 

and B\ is the 5x5 Hermitian matrix operator on the left hand side of (15.1 Pel ). The 
block-decomposition of Q is denoted Q = (Q mn )m,n=i,2,3- 

(57) Qi, Q2, Q3 are constant 9x5, 4x9, 4x4 matrices. Their nonzero entries are (C 
is the complex conjugation operator): 

($1)51 = -1 (Qi)n = -1 «3i)93 = 1 

(#a)i9 = -l-C ($2)27 = C (^2)28 = -1 

«?3)ll=-2 (03)22 = -1 
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Observe that $ 3 < 0. Let$(t) = ($ mn )m,n=l,2,3 be the C 5 © C 9 © C 4 block matrix 



Q{t) = ' 

(58) B M are the constant, diagonal, C 5 © C 9 © C 4 block matrices 

B = diag (17, U, U, U, V) © 1 9 V © 1 4 U , & 

where {7 = and V = ^ + Note that B° = 1 18 , B 1 = B 2 = 0, B 3 > 0. 

(59) Let s(x) be the smooth function that vanishes when x < and is equal to e~ x l x 
when x > 0. Let -0 = ip(q) : K 3 — ► [0, 1] be the smooth cutoff function 

*(3-||^) s(l-\q*-l\) 



*p._d_ 
1 dqv 



V>(q) 



«(3 ~ 1**1*0 + s 0!^ " I) <I - l<? 3 - !|) + s 0<? 3 - II 
where C= (£\£ 2 ) = (?\g 2 ).Let 



/C = 7J 3 | f |(0)x (il) c Q = 7J 4 | f |(0)x (0,2) 

By construction, supp R 3 ip C K, and %p is equal to 1 on Z) 5 1 a. | (0) x (§,§)■ For each 
integer R> 0, the bound ||V'||c H (R 3 ) 1 is independent of a and 21, see (SI). 

(510) Define 

M»(q,E) = ipW(q,E) + (1 - 
H(q,B)=1>Q(q,S) + (l-1>)Q(t) 
h(q,E) = H(q,E)E + ^Src(q) 

where ip = ip(q) is given in (S9). 

(511) If and S" (2) are both smooth solutions to B E = Q E + Src, see (S4), then 
their difference T = E^ — E^ 1 ' isa solution to 

B(q,E^)T = GT 

GU = £\ s=Q (Q(q,^)( S n)-B^q, S n)^ + Q(q, S n)E^) (8.4) 

where G(q, E^\ E^ 2 \ d q E^ 2 '^ acts on 7r _1 7?. linearly over R. The bracketed ex- 
pression in ( 18. 4t is affine linear in s. The operator ^ | selects the coefficient of s. 

Definition 8.1. Each entry to the left of the vertical bar is a generic symbol for a poly- 
nomial (with complex coefficients) in the (components of the) quantities to the right and 
their complex conjugates. 
,-1 



J 
J 

n 
n 

Qk 
Qk 



u 

1 ( That is, a generic symbol for a complex number) 
linear over R in 3^(0) 
linear over R in 21, e, A 

u~ l , 21, S, e, A, &(k)k=o...K+l } and their first derivatives 

u , 21, u, Sk, e, A, \P(k)k=o...K+i, and their first derivatives. 
It has no constant term as a polynomial in W(k) and its derivatives. 
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and 



There is no subscript K on the generic symbols J , J ' , 7i, TC because they represent 
polynomials that are required, by definition, to be independent of K. Precisely, neither 
their coefficients nor their degrees depend on K. By contrast, the presence of the sub- 
script K on the generic symbols Qk, Qk indicates that they represent polynomials that 
are allowed, by definition, to depend on K in an arbitrary manner. Precisely, their co- 
efficients and degrees may be functions of K. 

Above, S K is defined by S = - Ef=o(^) fe %? ( - k+1) u k+1 + ^tSk, where as before 
- = -i + Jr, see (l44t and ( t5~4l >. 

Proposition 8.1. 

B»(q,Q) = Jl" + u-'H + vT 2 Q K (8.5a) 

Qm(?,0) = Qm{q) + u^H+u^U + u~ 2 G K (8.5b) 

Q2n(q,0) = Oan(g) + n+ H + u- x Q K (8.5c) 

Qan(g,0) = 03n(<z) + Ht^+u-^J + u- 2 Q K (8.5d) 

and 

Srci(g) = u~ (K+2) gK (8.6a) 

Src 2 (g) = u- (if+2) g,f (8.6b) 

Src 3 (g) = u-( K+ Vg K (8.6c) 

B"(?) = u" 2 ^ (8.7a) 

Qin(«) = (8.7b) 

2 «(g) = J (8.7c) 

Q 3n (g) = ^ 2 i7 (8.7d) 

Remark 8.1. This proposition is a detailed examination of large |«| behavior of the con- 
stituents of the symmetric hyperbolic system d8.3ab . To convey its significance, it is 
helpful to suppress all but the derivatives in d8.3at and analyze the caricature scalar 
ordinary differential equation 

H u J)£f = q(«,/)/ + s(«) (8.8) 

In this remark, u plays the role of q°. Suppose f(u) is a solution to this equation on 
(— oo, T), T < 0, with asymptotic data lirriu^-oo f(u) = and h(u, f(u)) > 0. How 
can we estimate f(u)l For all u\ < T, 

/w-ri^fd.*^)^ (8.9, 

J-oo \Ju b(s,/(s))/b(ii,/(ii)) 

If there were constants ^4 > and 6 > c > with 

b(«,/(u))>6 q(u,/(«))<c|u|- 1 \s(u)\ < A\u\- 2 (8.10) 
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for all u € (—00, T), then ( 18.9b would imply 

sup \.f(u)\<-^-^- (8.11) 

«£(-oo,T) — c \i 1 

The apparent difficulty is that the functions in ( 18.10b depend on the solution f(u). How- 
ever, if it can be shown that a strictly weaker bound than ( 18.111 ). say ( 18.111 ) with A 
replaced by 2A, implies ( 18.101 ). then an open-closed argument justifies ( 18.1 11 ). More 
precisely, one would first cutoff —00 by a finite value, argue by continuity, and then 
remove the cutoff. 

To apply this reasoning, assume, in analogy with ( I8.3al ). that b, q are affine linear in /: 

b(«, /) = b(«, 0) + b(«) / b(tt) = ($b)(u, 0) 

q(u, /) = q(u, 0) + q(«) / q(u) - ($q)(u, 0) 

Also, in analogy with d8.5al ). d8.5bb . d8.7ab , d8.7bb . assume that there is a constant e > 0, 
so that 

|b(u,0) -M\ < e|u| _1 |q(«,0)| < e|u| _1 \h{u)\ < l^l" 1 |q(u)| < M -1 

For convenience, suppose b 7 = 1. The last inequality in (18.101 ) is an analog of (18.6al ). If 

e, A, ITI" 1 are sufficiently small, (8.12) 

then ( 18.111 . with A replaced by 2 A, implies ( 18.101 ), with b = \ and c = \. It follows 

from an open-closed argument that ( 18.111 ) is a genuine estimate for /(it). 

To interpret ( 18.121 in the light of our analogy, observe that the generic symbols H, Ti. in 

the second column in ( 18.5al ). ( 18.5bb can be made small by making i'(O) (equivalently, 

DATA) and the angular scaling parameter a small. 

We conclude the present discussion with the following remarks: 

• The analog of the step from ( 18.101 ) to ( 18.1 lb for the system ( 18.3ab is provided by the 
energy estimate. 

• Neglecting Q for the moment, d8.5db . d8.6cl ). d8.7db are similar to d8.5bl ). d8.6al ). 
d8.7bl ). since |£| _1 + u^ 1 is 0(u~ 2 ) as u — > —00 uniformly for u in a compact 
set. The interpretation of d8.5cb . d8.6bb , (18.7cb is different, because (18.8b is not the 
appropriate toy model problem for the equation satisfied by 3%. In fact, 32 satisfies 
an ordinary differential equation along the short integral curves of L, so that less u 
decay is required. 

• The inequality for q in d8.10b remains true if we add any non-positive constant to 
q. Analogously, the matrix Q{t), see (S7), appearing in (18.5cb . (18.5db . has only non- 
positive eigenvalues. This is implicitly exploited in the proof of the energy estimate. 

Remark 8.2. In d8.6ab all but the last component of Srci are actually u~^ k+3 'Qk- It is 
to accommodate the last component that we truncated the formal power series [ <P ] at 
K + 1 rather than K, see (S2). 

Proof (of Proposition I&71 ). The proof is by direct verification, using (S4), (S6), (S7), 
Proposition 15.21 Remark 15.41 and Definition 15.31 See the Supplement to Proposition 
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18.11 (Appendix IOT>. To give the flavor, we schematize the calculations for a few repre- 
sentative cases. Let C be the complex conjugation operator. Now 



matrix 


component 




B i (g,0)-P 1 

Q22M)-022(<?) 

Q33M) - $33(9) 


(4,5) 
(6,5) 
(1,1) 


-wi(O)<7-wi(O) + i0* 
(f + &9k) - (— 



in agreement with (I8.5ab . (I8.5cb and ( I8.5db . 

To verify d8.6ab . note that Srci = Qk has no constant term as a polynomial in 
<P(k) and its first derivatives. This follows directly from the definition of Src and the 
properties of / given in (S4). If S is replaced by Sk , see definition ( 18.11 ), then Srci = 
Qk- There is an overall u~( K+2 >, by construction of the formal power series solution 
[ & ] . This implies ( f8T6ab . □ 

To put d5.7bb in the form required by Proposition l7.7l we use 

(51) LetS - ' = (Eh, Eh, Eh) be the field given by 

(Eh, Eh, Eh) 

= (2/1,2/2,2/3) © (si,S2,Si,S5,Pl,p2,p3,P4,P7,Ps) © (s3,P5,P6,P9) 

where = (s,p, y) is the constraint field. There is a permutation matrix tt so that 

(s,p,y) = w(Eh, Eh, Eh)- 

The field E% = (Eh, Eh, Eh) takes values in n^U C C 3 © C 10 © C 4 . 

(52) System ( I5.7bb is abbreviated as A(q, V)^ = f (q, d q ^)^ (see, Convention 
18.11 ). Some of its properties are discussed in Remark |531 System d5.7b| i is equivalent 
to the linear, homogeneous symmetric hyperbolic system 

B(q,*)E* = Q(q,*,d q *)S* , B = B"^ 
B M (g,<F) = K- l A»(q,&)K 
Q(q^,d q V)=n- 1 f(q,V,d q V)n 

The transformation Q acting on 7? ~ 1 1Z is linear over R. Moreover, B M depends affine 
linearly over R on <P, and Q depends affine linearly over R on W © d q <P ' . 

(53) The three by three C 3 © C 10 © C 4 block-decomposition of B is 

B = diag( J B 1 ,B 2 ,S 3 ), B 2 = 1 1Q L, B 3 = 1 4 N 7 

and B\ is the 3x3 Hermitian matrix operator on the left hand side of ( I5.12cl ). 

(54) Q x , Q 2 , Q3 are constant 10 x 3, 4 x 10, 4 x 4 matrices. Their nonzero entries are 
(C is the complex conjugation operator): 

($1)5,1 = 1 ($2)1,9 = ($2)4,10 = -1 

($3)1,1 = -1 ($2)1,10 = ($2)2,8 = ($2)3,7 = ($2)4,9 = -C 
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Note that Q 3 < 0. Let Q(t) = {Q 



mn Im.n- 



= i 2 3 be the C 3 



block matrix 




(S5) jB^ are the constant, diagonal, C 3 © C 



10 



block matrices 



P = 1 3 U®1 W V®UU , $ = $"-2 



with £/, V" as in (S8). Note that $° = l ir , J3 1 = $ 2 = 0, J3 3 > 0. 

Definition 8.2. Each entry to the left of the vertical bar is a generic symbol for a poly- 
nomial (with complex coefficients) in the (components of the) quantities to the right and 
their complex conjugates. 

£ B | u~\ 21, S, e, A, i'(O), ^-^(0), and first derivatives 
Ql \ u~\ 21, S, e, A, «p-(0), <f - S>-(0) 
J | /fte (?", but it has no constant term as a polynomial in <P — \P(0), d q (fP — ^(0)) 
21, u, Sq, e, A, ^(0), and first derivatives 

where So is defined by S = — 2t 2 u + it" 1 So, see j4.4\ and ( 15. 41 ). 
Proposition 8.2. Suppose (SI) to (S5). Then 



V 

hn(q) 
hn(q) 



u-*Ql 

u- 2 g* 



3n(g) +W" 2 ^ 



(8.13a) 
(8.13b) 
(8.13c) 
(8.13d) 



Moreover, 



-3, Si, S 2 , Pi, p 2 , P3 = U £ 



s 4 , s 5 , p 4 , P7: ps = " 1 Q i + uQ\ 

Finally, it is a consequence of d5.7b| ) f/iaf 



(8.14a) 
(8.14b) 





( 


S5 




Pi 




P7 




W 


V 



h~Pn)+ e(p 6 -p 3 )\ 
*e(p 4 -p 5 ) -ie(ft -p 3 ) 


MPS -P 3 ) ~ A (P4 -P 5 ) 
<Mp 6 ~ P3) + A (P4 - P5)/ 



(8.15) 
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Proof. The first part, ( 18.13b . and the last part, ( 18.15b . follow directly from Proposition 
15.41 Remark l5~4l and Definition l5.3l It is entirely similar to the proof of Proposition l8. II 
To prove ( 18.141 , write \P = •P'(O) + (9— ^(0)), and consider each object on the left 
hand side of ( 18.14b as a polynomial in \P — ^(0) and — ^(0)), with coefficients 
possibly depending on •f'(O) and ^^(O) (see, Proposition 15.31 ). The idea is that the 
constant term of this polynomial is of the generic form u~ 1 Q\ Everything else is of 
the form Q\ or u Q\, respectively. The fact that the constant term is of the generic form 
u~ x Q^ is an essential part of the construction. It is the fact that i'(O) is built so that the 
first term in the formal power series of the constraint field, (0), vanishes (this follows 
from the vanishing of the formal constraint field and Remark l6TTT l. □ 

Everything we have done in this section so far was to prepare for 

Proposition 8.3 (Main Technical Proposition). Fix K > as in (S2). Suppose (Si) 
through (Sll) and (SI) through (S5) all hold. Let R > 4 be an integer. Set 

Y=(r, K, ||data|| c « + 2/c +6(q) ) (8.16) 

Let d be as in d740t . Fix c' 2 G (0, 1) and T e (-oo, -fl -1 ). There are constants 
Cq(R) € (0, 1) and Cj(Y) G (0, 1), non-increasing in all their arguments, such that 
Parts 1, 2 and 3 below hold whenever 

\a\<c 6 (R)c' 2 , ||DATA|| CR+ 4 (Q) <c 6 (J?)c' 2 , iTp 1 < c 7 (Y) c' 2 (8.17) 

Part 1. The system M(q, S) S = h{q, S) in (S10) satisfies (EB0) through (EB4) and the 
assumptions of Part 2 of Proposition \73\ with: 



(EB0) - (EB4) 


T 


P 


M>(g,6>) 


h(q,0) 






Q 


K 


(SI) - (S10) 


T 


31 




h(q,E) 




Q(t) 


Q 


K 



The table indicates that the symbols in the first row, appearing in the general (EB0) 
through (EB4), are given by the specific objects in the second row, appearing in (SI) 
through (S10). 

Part 2. If t < T and S : [t ,t + e) x R 3 -> n^n (e > 0) is a C°° solution to 
M(g, S) S = h{q 1 S) which vanishes identically at to, then 

M K+1 sup Je* Q) {~}(to) < (cj(Y)r 1 (8.18) 

Here, the energy Eq^ bt \ is defined as in the Refined Energy Estimate, Proposition \7.7\ 

Part 3. We distinguish three alternative systems, denoted by (Sysl), (Sys2) and (Sys3), 
that are given in the columns of the table below. This part of the Proposition applies to 
each of these systems individually. To evaluate the entries in the table, we require the 
following information. First, 

(Sysl): 6 = 2 Co € K 2 
(Sys2): b = l £ e-D 2 |£|(0) 
(Sys3): b = l Co€Z> 2 |f|(0) 
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Second, fix to < T and define the open set 

V = |J {(}xOK„,M) C M 4 . 

t&(t ,T) 

For (Sysl) and (Sys3) there is a single field S defined on V taking values in 7r _1 7?.. For 
system (Sys2), there are two fields, S' 1 ' and S^ 2 \ of this kind. The various fields satisfy 
the conditions: 

( i) They are C p and their derivatives of order < p extend continuously to V. Here 
p = oo for (Sysl), (Sys3) and p — 1 for (Sys2). 

( ii) They are solutions to 



M(q,S)S = h(q,S) 

B(q, = Q(q, 5«) + Src(g) 

B(q,~)~ =Q(q,~)+Src(q) 



See (S4) and (S10). 
(Hi) They vanish when q 3 < i. 
f/vj For a/Z t G (to, T), f/zey satisfy 



sopgj fo)1>t) {sW}(t) < c 6 (ij)4 



/or (Sysl) 
/or (Sys2) 
/or (Sys3) 



/or (Sysl) 
/or (Sys2) 
for (Sys3) 



To sfafe f/ie conclusions, recall the notation T = a ( 2 > — B ^' and the usage ^ for the 
constraint field associated to & = Sfn + tt E, see (S2). Finally, <Z™ = 7r i", as /« (SI). 



Conclusion 1: E(V), E^>(V), E^'(V) C £i/ 2 (0) C n^TZ M 31 . 

Conclusion 2: 77ie assumptions (REO) through (RE12) /zoM, vv/f/z (RElla) /or (Sysl) and 
(REllb) for (Sys2) and (Sys3), provided that the symbols in the first column of the table 
below (appearing in the general (REO) through (RE12)j are given by the specific objects 
in the other three columns. 



(REO) - (RE12) 


(Sysl) 


(Sys2) 


(Sys3) 


1 = (to,**) 


(*o,r) 




(to,T) 


b 


2 


i 


1 


(Pl,P2,Ps) 


(10,15,6) 


(10,15,6) 


(6,18,8) 


M"(g) 


M"( 9 , ~(q)) 




B"(?,#(g)) 






G{q,EM,~W,d q S(*)) 


Q(<7,<M^) 


Src(g) 


■0Src(g) 








9{q) 


S{q) 


T(q) 


**(«) 

r 


w 






W{t) 


Q(t) 


Q(t) 


M) 


R 


R 








Cl 


c' 2 








c 2 


c' 2 


c 2 


J 


K+l 


> 


> 


Co 


Co 


Co 


Co 
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Conclusion 3: For (Sys3), if in addition to + 1 < T, then 

sup |t|Snpg> < Y (l + sup |*|SapgL lt) {S}(*)) (8.19) 

Proof. We begin with a warning. 

First Warning. In the course of this proof, we produce a finite chain of smallness as- 
sumptions on cq(R) and Cf(Y). It is essential, for the purpose of showing that the far 
field expansion is truly an asymptotic expansion to a classical solution of d5.7al i. that 
these smallness assumptions depend only on R and Y, respectively. To give a repre- 
sentative example, suppose quantity <r cq(R). Then there is a legitimate smallness 
assumption on Cq(R) making, say, quantity < 1. By contrast, there is no legitimate 
smallness assumption associated to quantity <y Cq(R). We can take a more relaxed 
attitude to the system i5Jbi . because it is only necessary to demonstrate uniqueness. 

Convention 8.5. In this proof, the constants of proportionality in <^ and <y are always 
non-decreasing in R and the components of Y, respectively. 

Overall Preliminaries. For all n > and < k < K + 1 and j3 E N$ with \/3\ < 1 + R, 
the following estimates hold on (— oo, T) x Q: 

\d^\ = |2t| 6/30 < \a\ \d^u\ < 2 

\dPe\ < % \a\ \d X\ < f \a\ 

\d fi S\< R l \d^3 K \< Y l 

Only the estimates on ^(0) and \P(k) require discussion. They follow from Proposition 
[Ql and (l8T6t and (fQ7t . 

It follows from the estimates just above, the product rule and ( 18.171 ), that for all 
a G No with |a| < R, 



n = 


0,1 


\d a (u- {n+1 


] Qk)\ 


<Y 


IT]- 1 < Y c 7 (Y)c' 2 


(8.20a) 


n = 


1,2 


\t\ K+n \d a {u- {K+n] 


Qk)\ 


<Y 


1 


(8.20b) 


n = 


0,1 


\t\ n \d a {u 


~ n H)\ 


<R 


c 6 (R) c'. 2 


(8.20c) 


n = 


0,1 


\t\ n \d a (u- 


n H)\ 


<R 


\a\ <r c 6 (R) c' 2 


(8.20d) 


n = 


0,1,2 


\t\ n \d a {u 


' n J)\ 


<R 


1 


(8.20e) 






\t\\d a (\tr+ u 


' l )J\ 


<R 


i^r 1 < R c 7 (Y)d 2 


(8.20f) 



at every point of (— oo,T) x Q. In this instance, the constants also depend on the 
particular polynomial represented by the generic symbols. Observe that in the second 
inequality, one does not use the property that Qk has no constant term as a polynomial 
in •P'(fc) and its derivatives. 

Second Warning. It is crucially important that whenever < r appears in an estimate (for 
example, d8.20ct , d8.20d| ), d8.20e| i, d8.20fl )) that the generic symbol on the left hand side 
has no subindex K, see Definition 18. II On the other hand, whenever <y appears (for 
example, d8.20al ). d8.20bl )). the generic symbol on the left hand side is allowed to carry 
a subindex K. 



\d p u~ n \ < {R . n) \t\~ n 
Wfl)|<« c 6 (R)c' 2 
\d^(k)\ < Y 1 
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Preliminaries for Part 3. For Part 3, it is necessary to supplement the Overall Prelimi- 
naries. Let (Vi, V2) be the open cover of V given by 

Vi =Vn(Oo,T) x Q), V 2 = Vn ((to,T) x (K 3 \/C)). 

The sets Q, /C are defined in (S9). 

• For (Sysl), observe that the Overall Preliminaries apply to Vi . On V2, we have ip = 0, 
and the equations simplify, see (S10). The estimate ||V- , llc H (K 3 ) ISr 1' see ( S9 )' W1 U be 
used on the transition region for -0. 

• For (Sys2) we have V = V\. In this case, the Overall Preliminaries will suffice. 

• For (Sys3) we also have V = Vi. However, in addition to the Overall Preliminaries, 
we require the estimates 

|!F-#(0)| , |0,(ff-!P(O))| < Y 1 (8.21) 

\t\\u- 2 g\\ , \t\\d q {u^Gl)\ < Y \T\- l <yc 7 {Y)c' 2 

n = 0,l: |t|"| U -(" +1 ^»| <y iri" 1 ^ c 7 (F)c 2 , 

|W| < M<c 6 (i?)c 2 

|t||(|i|- 1 + «- 1 )^l <y irr^yc^y)^. 

on V. Estimate ( 18.211 ) follows from 

<P-V(0) = K j:\±) k V(k)+irS (8.22) 
fc=i 

and condition (iv) in the Proposition. The rest are consequences of ( 18.211 ) and the 
Overall Preliminaries estimates. 

Proof of Part 1. One verifies (EBO) and (EB4) by direct inspection, apart from the in- 
equality in (EB1). However, 

±<M°( g ,~)<2, M 3 (q,S)>0, (8.23) 

for all (q,E) £ (— 00, T) x R 3 x B 2 (0) when c e (R) and Cr(Y) are small enough. 
Here B 2 (0) C M p as in (EBO), with P = 31. Th e first inequality in ( 18.23b completes the 
verification of (EB1). The second inequality in ( 18.231 is used for the first supplemental 
hypothesis in Part 2 of Proposition 17. 3 1 
To check ( 18.231 ). recall from (S10) that 

M»(q, S) = i> B»(q, S) + (1 - ^) J8" (8.24) 

is a convex combination of B M and B M on (— 00, T) x R 3 x i?2(0). It suffices to verify 
( 18.23b for and B M separately. For B M , see (S8). For B'' it suffices to verify 

^ < 1 + ^/3 < 2, \ < 1 + ^(1 + ^/3) < 2 (8.25) 

for q 6 supp R3 ip C Q, see (S6) and Remark |5~4l Here, fa is one of the components of 
& = + 7T S 1 = (/, w, z), see (S3). To check ( 18.25b , note that 

A"+l 

If I < |^(o)| + ^ i; ^ |#(*)| + 
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The three terms are respectively <r cq(R) and <y pL- <y C7(y) (see, Overall Pre- 
liminaries) and < 2. By the choice of Ce(R) and cj(Y) (see, the First Warning), we can 
make 

\&\ < 3 , when (<?, S) e (-oo, T) x Q x S 2 (0) (8.26) 

Consequently, I/3I < 3, and therefore, (18.25b holds because j^y < rry < 5, see ( 17.401 ). 

To validate the second supplemental hypothesis in Part 2 of Proposition 17.31 it is 
necessary to show that h(q, 0) = ipSrc(q) = for all q e (— 00, T) x R 3 with q 3 < i. 
To do this, observe that <Pk = there, see (S2). 

Remark 8.3. Later on, in the proof of Part 3, we need the analogous inequalities 

±<B°(q,&)<2, B 3 (q^)>0 

for all (q, S) <G (—00, T) x R 3 x #2(0), with the same smallness assumptions. These 
inequalities can again be reduced to ( 18 .25b . 

Proof of Part 2. To prove d8.181 l. rewrite M(g, E) S = h(q, S) as 

a t 5(q) = (M?(q,S))- 1 (- M i (q,S)d i S + H(q,S)S + i>Src) (8.27) 

t=l,2,3 

Here, M°( • , S) is invertible on an open neighborhood of {to} x K 3 in the set [to, to + 
e) x R 3 . This is a consequence of ( I8.23l l, and the assumption S(to, • ) = 0. 

By repeated differentiation of ( |8.27t with respect to t, we obtain an expression for 
dY l E{q), for any m > 1. Restrict the result to {t } x R 3 and simplify it using the 
assumption d^S(to, • ) = for all (3 G Nq with fio = 0. In every surviving multi- 
derivative d@S(to, ■ ) we must have 1 < /?o < m — 1, and each one is recursively 
expressed using d^S(to, ■ ) with 1 < n < m — 1. This procedure generates an explicit 
expression for d™S(to, ■ ) in terms of the quantities in (18.281 just below, and their 
derivatives. 

By differentiation with respect to the remaining coordinates q, we find (inductively) 
that for every a £ Nq with \a\ < R, the function d a S(to, • ) is a polynomial in 
(M°((i , • ), O))^ 1 as well as derivatives of order < R - 1 of 

M"((*o, -)>0) , H((t , -),0) , M"(to, •) , F(t , •) , V'Src(t , ■) (8.28) 

This polynomial has no constant term as a polynomial in ?/>Src(to, ■ ) and its deriva- 
tives. This has two consequences. Consequence A, d a S(to, ■ ) vanishes on R 3 \ JC, by 
the support of ip. Consequence B, 

\t \ K+1 \d a S(t ,q)\ < Y 1 

for all q G Q. To verify this inequality, check (using ||^||c J? (R 3 ) ^-R 1 m ( S9 )> me 
first inequality in (18.23b . Proposition 18.11 and the Overall Preliminaries) that the matrix 
(M°((i , ■ ), 0)) _1 and the derivatives of order < R — 1 of all the terms in ( 18.28b are 
bounded in absolute value on Q by <y 1. At this point, we have \d a S(to, q)| <y 1. 
To get the stated decay, use ( 18. 6b (even though one can get a better result). It is here 
that one exploits the fact that the expression for d a S(to, ■ ) has no constant term as a 
polynomial in ijj Src(£o, • ) and its derivatives. 
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The proof of Part 2 is completed by combining Consequences A and B with 

E 0(i ,2,t ){ S K t o) ~-R ( Su Pe?(f ,2,to){ S K*o)) 

and making a suitable choice of C7(Y) (see, the First Warning). 

Proof of Part 3, Conclusion 1. Follows from condition (iv) in the Proposition, by suit- 
able choice of Cq(R). For (Sysl), we also use R > 2 and the Sobolev inequality (17.43b . 

Proof of Part 3, Conclusion 2. To start with, we check that (RE10) and (RElla) or (REllb) 
hold, when Cq(R), c?(Y) are made sufficiently small. We will freely use the Overall 
Preliminaries, the Preliminaries for Part 3, Propositions 18.11 [iOl and the inequality 

E^ oM {f}(t) <r (Su P ^ o . M) {/}(t)) 2 . (8.29) 

• (Sysl): Let t G (t a ,T). For (RE10), we have 

|t| 2if+4 ^ (£oj2 , 4) {VSr Cl }(t) = \t\ 2K+A E^ t) {^u-^g K }{t) < Y 1 

Therefore, the left hand side is < (c7(y))~ 2 , when Cr(Y) > is small enough (see, 
the First Warning). Similarly, for ?/>Src2 and i/jStcs. In these two cases, one could 
get a better decay estimate, but we don't need it. 

For (RElla), we verify the first and second inequalities, the other two are similar. The 
second goes 

\t\ 2 Eg {M) {H ln (q,E)-Q ln }(t) 

= \t\ a ^ o ,2,t)MQin(q,0)-Qin)+i>Qi n (q)S}(t) 

= \t\ 2 Egu ,2,t)Mh^ + + + *& s } (*) 

<R\tf ^( & , a ,*){V'(^ + iW+^S)}(i) + |t| 2 ^ (eo>ait) {V»^}W 

The first term is (cg(i?)c' 2 ) 2 , the second is <y (c7(F)c 2 ) 2 . Here, we use condi- 
tion (iv) in the Proposition. By the First Warning, the result is < (c' 2 ) 2 if Ce(R) and 
07(F) are small enough. The first goes 

\t\ 2 E* {M) {M»(q,E)-W}(t) 

= \t\ 2 E§ (io , 2 , t) {^(in + ^Qk) + 1> S}(t) 

The estimate is completed just as above. 

• (Sys2): There is nothing to check for (RE10). For (REllb), 

l^uP^ , M) {B^( 9 ,- (1) )-P-}(t) 

= 1*1 s »PoLm) { i B " (9- °) - + (?) } (*) 
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which is < c' 2 when cq(R) and Cf(Y) are made sufficiently small. It is important 
here that V = Vi, see, the Preliminaries for Part 3. For the second, third and fourth 
parts of (REllb), observe that, by (18.4b . 

(G-Q)n= (Q(q, 0) - Q(t))n + (Q(q) S«)/7 

+ &U( - B"(g, »n) ^ + Q(g, sn) . 

Therefore, for the second, 

\t\ Sup^ 0il;t) {G 1 „(g, ~ W ~< 2 \ d q ~W) Q ln (t)}(t) 

= + in + ^Qk + + ^ Jdq ~^ + 1 jew}® 

which is < c' 2 when c$(R) and ct(Y) are made sufficiently small. The third an fourth 
inequalities in (REllb) are checked in the same way. 
• (Sys3): There is nothing to check for (RE10). For (REllb), 

l*|Sup^ , ljt) {B' 1 ( g ^)-p' i }W = |t|Supg^ ;li4) {^^}W 
1*1 SnP^ , llt ){Qln(?, d q V) - Q ln }(t) = \t\ Sup^ ol t) { ^ } (t) 

snpS eo ,i, t ){4n(?,^W - Q 2n }(t) = s» P ( ^ olt) {n + iG*}(t) 

\t\Su P { °l o l t) {Q 3n (q,V 7 d q V) - Q 3n }(t) 

= l*l*pSLi,*){(iir + =^ + ^'K t ) 

which are all < c' 2 , when cq(R) and 07(F) are made sufficiently small. It is important 
here that V = Vi, see, the Preliminaries for Part 3. 

We are now finished checking (RE10) and (RElla) or (REllb). 
Next, we check (RE1). In order, (RE1) follows from 

• (Sysl): inequalities (18.23b . since S(q) € Bi (0) C M 31 for all q e V by Conclusion 1. 

• (Sys2): the discussion of ( 18.25b . Conclusion 1 and V = Vi. 

• (Sys3): Remark 1831 Conclusion 1 and V = Vi. 

To check (RE12), we can assume |<?| < 3 on Vi, by (18.26b and Conclusion 1. By Propo- 
sition [7T6] and (SI), this implies (RE12) for (Sys2) and (Sys3), because V = Vi. For (Sysl), 
we use convexity, see d8.24b . The term is again handled by Proposition |T6j using 
the support properties of ip. It therefore suffices to verify that 

(9 M P^>0 on (dV) n ((t ,T) x R 2 x (0,2)) 

with 8 defined by ( 17.38b . This is a consequence of 9(U) = ki\u — u\/\u\ 2 > and 
6{V) = fciHul" 1 - Kl _1 | > 0, see (S8) and ( PT40I 

(RE5) holds, by condition (i) in the Proposition, for (Sysl), (Sys2) and (Sys3). In partic- 
ular, for (Sys2), G(q, S^\S^\d q E^) is C°. 

(RE4) holds, by condition (ii) in the Proposition, for (Sysl), (Sys2) and (Sys3). For (Sys2), 
see (Sll). For (Sys3), recall that solves d5.7bb because W solves d5.7al ). 

Note that condition (iii) in the Proposition and (S2) imply that \P = + tt — = 
and ^ = when q 3 < |. In particular, Src = there. These facts imply (RE6) for 
(Sysl), (Sys2) and (Sys3). 
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The remaining assumptions, (REO), (RE2), (RE3), (RE7), (RE8), (RE9) are verified by 
direct inspection. 

Proof of Part 3, Conclusion 3. We have to prove 

sup |t|Supgl lt) {S*}(t) <y (l+ sup |t|Supg> M) {S}(f)) = « 

It follows from the Overall Preliminaries, the Preliminaries for Part 3, as well as (18.14b 
and ( 18.151 ) (see, PropositionlOi that on V = Vi, 



(A) 


|0* 


<R 


1 


by DefinitionliHl 


(B) 


|#-#(0)|, |0,(#-#(O))| 


<Y 




by iS.22) 


(Q 


is! 


<Y 




by (B), DefinitionEIl 


(D) 


l^iUS'sl.W.H.bil.bal.lPs 


<Y 




by (A), (C), (|8.14a|) 


(E) 


|»« 


<Y 


K 


bv (A), (O. d8.14aD, d8.14bD 


(F) 




<Y 




by (E), DefinitionO 



For each point (ti, qi) G V with ti > to + 1, consider the line segment 

Seg = ((ti, qi) - 0, 0, 1)) n {g e M 4 | g 3 > 0}. 

We have Seg C V. To see this, view V as an open disk bundle over the (t, u) -rectangle 
(to, T) x (0, 1). The projection of Seg to the (t, u) plane is injective and contained in the 
base, because t\ G (to + 1, T). At each point in the image of the projection of Seg, the 
corresponding point on Seg is contained in the fiber, because the radius function r(u, u) 
(see (|7.37t ) is a decreasing function of u on the base for fixed u = t — u, and because 
the endpoint (ti, qi) is contained in the fiber, by assumption. 

By Conclusion 2, is a C°° solution to RS^ = which vanishes when q 3 < 
i. In particular (18.151 ) holds. Recall L = e?,(-^> + where <P = (e,j,w), and 

h < e% < 2 (see (RE1)), and that Seg is an integral curve of L. The last three sentences, 
Seg C V and (F) imply, by integrating the equation for p^ in (18.15b along Seg, that 

\Pi\SYK\t\- 1 (8.30) 

on Seg for each endpoint (ti, qi) G V with t\ > to + 1. 

Finally, integrating the remaining equations in d8.15t along Seg, and using (|8.30t 
and (D), we obtain | S4 j , |ss|, |f>7|, \ps\ <r /s: | ^ | 1 on all admissible segments Seg. 
Therefore, \ ~$\ < Y k^I" 1 on V when t G (t + 1,T). □ 

8.2. Construction of classical vacuum fields. 

Theorem 8.1. Let (£,u, u) be the usual coordinates on the truncated strip 

Strip(l, A) = R 2 x (0, 1) x ( - 00 , -A" 1 ) 
of width 1, for each A > 0. Suppose < |2l| < \a\. Assume the functions 
DATA f7 (^, u) : R 2 X (0,oo) -> C 
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a G { — , +}, are smooth, vanish when u<\, and are Pole-Flip compatible 

DATA CT = Flip a. • DATA ^ (8.31) 

for all € (K 2 \ {0}) x (0, oo), see Definition \6.2\ Let [ r P cr } be the formal power 

series solution corresponding to DATA' 7 . Fix integers R> 4, K > and an e € (0, ^). 
Set 

B = (R,e) , C =(r, e, K, max l|DATA ff || cR+2K+6 ( c(a ^ 2) )) 

w/iereC(a, 21, 6) = D 4 | f |(0) x (0,6) for each b > 0. 
Let 

b = b(B) c = c(C) 

be constants in (0, 1). Ifh and c are f7zac/e sufficiently small depending only on B and 
C, respectively, then the Existence and Uniqueness statements below hold whenever 

< |2l| < |a| < b, max ||DATA ,T || C H + 4 (c(a a 2)) < b (8.32) 

<re{-,+} 

Existence: 

Part 1: There exists a pair ('f' - , !f r+ ) of Pole-Flip compatible ^-fields 

V a : Strip(l,c) -> ft, 

which are both solutions to (15. 7 al l, vanish when u<\, extend with their first derivatives 
continuously to Strip(l, c) ana? satisfy 



lira |w| c sup 

>->— °° agMj: |a|<l 



a a (^-^(0))(-,u)|| O0(c(O)8lil)) =0 (8.33) 



Part 2: 77ie constraint fields )", (^H")" associated to the fields in Part 1 vanish, and 

= (M a ,a+w- A/ ^-, X a , a + U - M f+) 
are a pa/r of Pole-Flip compatible vacuum fields (see, Definition \2.2\l with initial data 

Part 3: The fields in Part 1 are actually of class C R ~ 3 , and extend with their derivatives 
of order < R — 3 continuously to Strip(l, c). Moreover 



sup |u| K+1 sup 



K 



^(fc)(.) 



i 

C°(c(a,8t,l)) ~~ C 



|a|<ft-3 

(8.34) 

Uniqueness: Asrame (<f r_ , | f r+ ) a«(i (^ — , have all the properties listed in Part 1. 
Then they coincide on Strip(l, jt^) C Strip(l, c). 

Remark 8.4. Part 1 asserts the existence of a solution on I U II U III. However, Unique- 
ness in Theorem 18.11 refers only to I. We actually prove uniqueness on I U II. By a 
standard finite speed of propagation argument, which we do not carry out, the domain 
of uniqueness can be extended to III. 
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Remark 8.5. The DATA CT are given for u E (0, oo), just for convenience. By construc- 
tion, the restriction of ^(O) to u G (0, 1) depends only on the restriction of DATA" 7 to 
u £ (0, 1), see Definition l5.ll It now follows from Uniqueness in Theorem l8.1l that , 
ty + are determined on Strip(l, by the restriction of DATA CT to u e (0, 1). 

Proof. Theorem l8.1l is formulated in the coordinate system (£,u, u). Almost the entire 
proof, however, is given in the coordinate system q = (t, £, u), where t = u + u. 

We assume (SI) through (Sll) and (SI) through (S5), where a, 21, K in (SI) through 
(Sll) are identified with their occurrences in the statement of Theorem 18. II and where 
DATA in (S2) is identified with either one of DATA CT , a = — , +. By direct inspection, our 
assumptions and identifications are consistent, when we make the legitimate smallness 
assumption b < 1CP 3 , that ensures that a, 21 satisfy (SI). This condition is subsumed in 
d837l i below. 

Convention 8.6. For the entire proof, c 3 ( • ) and C4( • ), as well as c 6 ( • ) and c 7 ( • ) are 
defined as in Proposition [T7] (Refined Energy Estimate) and Proposition [83] Further- 
more, cg(i?) > 1 will always denote a constant, such that the Sobolev inequality 

holds for all (£ , b, i)el 2 x[l,2]x (-co, — and all vector valued C R functions 
/. See (17.431 for the Sobolev inequality and ( 17.401 for the definition of 5. 

The smallness condition on b. Set J = e, with e as in Theorem l8.ll and 

x = faJoMMM) 
x* = (o,Jo,|0i|,K? 2 |,|03|) 
x = (o, j 0) I^U^U^I) 

where Qi and Q i are fixed as in (S7) and (S4), respectively. Also set 

c' 2 (B) = i min {c 3 (X), c 3 (X*), c 3 (X)} e (0, 1) (8.36) 

The right hand side of d8.361 l only depends on B, which justifies the notation c' 2 (B). 
We impose the legitimate smallness condition 

b < min{l(T 3 , c 6 (i?) C 2(B)} (8.37) 

It is the only smallness condition on b in the entire proof. 
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The first smallness condition on c. Set 

Y=[R,K, max ||DATA <T || c <R+2K+6(g)) 

Observe that C(a, 21, 2) = Q where Q is defined in (S9). Therefore, Y depends only on 
C, and so does the right hand side of ( 18.381 ). justifying the notation T(C). We impose 
the legitimate smallness condition 

c * wmr-2 < 839 > 

There will be one more smallness condition on c, later in the proof. 

Convention 8.7. The system J8.3at corresponding to DATA " will be denoted ([83aji cr . We 
sometimes suppress the superscript a and simply write DATA and (18.3al i. in which case 
the discussion applies equally to DATA ff and J8.3a| > gr for a = — , +. 

Convention 8.8. In every application of Proposition 18. 31 the c 2 of Proposition 18 . 31 will 
be the c' 2 {B) of dQ6b . 

Remark 8.6. We have Flip » • (M a ,% + u- M W) = M a .% + u~ M (Flip a. -W) since 
the transformation Flip o_ is linear, maps M. a ,<A to itself and commutes with the matrix 
u~ AI . In particular, if ^ solves (I5.7at . then Flip^. • <P solves d5.7at . 
Observe that Flip^. • tf^- = ^^ a , because DATA CT are Pole-Flip compatible, by assump- 
tion. Therefore, Flip^ • + 7r E) = ^ CT + Flip «_ ■ EE). In particular, if EE solves 
diOal r. then tt" 1 Flip f • (ttH) solves d8~3al T' T . 

Convention 8.9. For the rest of this proof, we consciously abuse notation and write 
Flip a. • E for 7T _1 Flip o_ • inEE). See, (S3) for the definition of n. 

21 21 v ' 

Guide. The proof now proceeds through a sequence of 10 steps. Each step begins with 
a statement (in italics), that is then proven. 

Step 1. The assumptions of Proposition \8~3\ up to and including ( 18.171 ) are satisfied for 
all T < T(C), //DATA = DATA CT for a = -, +. 

By direct inspection. Recall that c-?( ■ ) is non-increasing in all its arguments. 

Step 2. For each to < T(C) there is a ti(fo) G (to, T(C)] and a smooth solution 

B to : [io,*i(*o)) x M 3 ^ tt^TZ S W 31 

to the system M(g, 5)5 = h(q, 5), with trivial initial data, EE to (to, ■ ) = 0, vanishing 
identically for q 3 < i, so that t\ (to) ^ T(C) implies either one or both of (Break) 1 or 
(Break) 2 (see Proposition 1 7. 3D . 

Apply Parts 1 and 2 of Proposition |7.31 in the context of Part 1 of Proposition ^. 31 with 
T = T(C). 

Remark 8.7. EE t() is a 1 -parameter family of solutions, parametrized by to < T(C). 
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Step 3. ti(t ) = T(C),forallt < T(C). 
For each £o G R 2 , we introduce the set 



Jfa,to) = Ue Mi(*o)) sup 2 ){S t0 }(r) < (c 6 (i?)c 2 (B)) 2 } 

rG [to,* 1 



It is an interval and closed as a subset of [to, ti(to))- By Part 2 of Proposition 18 . 3 1 and 
by (EH), 



l/^0(? ,2,to){- t o}( < °) - C7 (y) |T(C)| K+1 



< C6(gK(g) 



Therefore, to G J(£o, to). By continuity of the energy, i/(£o, to) contains at least one 
point different from to- For every t* G i7(£oi to), t* > to, the assumptions of Proposi- 
tion !8.3l Part 3, (Sysl) are satisfied with T = t*. It follows from Conclusion 1 that 



S f .(9)eBi/ 2 (0)cR 31 for all q G U te(to ,t*) W X 0(6>> 2 > *) (8-40) 

By Conclusion 2 and 

/v + 1 > Jo, 4(B) < c 3 (X), t* < T(C) < -l/c' 2 (B) < -l/c 3 (X) 

we can apply the Refined Energy Estimate (|7.42t in Proposition [777] m the context of 
(Sysl). Combining ( |7.42t with Part 2 of Proposition |8.3l one obtains 



2c 4 (X) 2c 4 (X) 



cr(y)M^ - c 7 (y)|T(c)|^+i 

c 6 (i?)c 2 (i?) c 6 (i?)c 2 (i?) 1 ■ 
- 2c 8 (i?) ~ 2 



for all t G (to , t* ). The second inequality is self-evident. For the third, use ( 18.38b again. 

The continuity of the energy Ef^^ 2 T \{St }(T) for r G [to,ti(to)) implies that 
( 18.41| > holds for r = t*, and, consequently, for all r G ^7(f , *o)- 11 follows that 
t7(£c>)to) is also open as a subset of [to,ti(to))- The set i7(£o:to) is nonempty, open 
and closed as a subset of [to,ti(to)), and we conclude »7(£o,to) = [^0,^1 (to))- The 
upshot is that ( 18.40b holds with t* replaced by ti(to), for all £o G K 2 , and therefore 
S to ([t ,h(t )) x Q) c Bi /a (0) C K 31 because 

[to,*i(*o)) x Q C [t ,ti(t )) xl 2 x (0,2) = |J |J {(}x%2,i). 

£ eR 2 te[t ,ii(*o)) 

(8.42) 

Now, (Break)]^ (see Step 2) is excluded. The inclusion ( 18.421 ). the fact that i/(£o, to) = 
[to, ti (to)) and the Sobolev inequality, ( 18.351 ), exclude (Break) 2 . By Step 2, we conclude 
that ti (t ) = T(C), for all t < T(C). 

Remark 8.8. A byproduct of the proof of Step 3 is: 

The inequality ( ED holds for all r G [t , T(C)) and £ G R 2 . (8.43) 



59 



Convention 8.10 (for Steps 4 and 5). We introduce a new field that is used in the next 
two steps. It is the restriction of S to to [t , T(C)) x W, where 

W = D S|f| (0)x(0,l) c R 3 

Consciously abusing notation, we will denote this new field by the same symbol S to . 
The new field S to is smooth and extends, with all its derivatives, continuously to 

[to,T(C)) x Wc M 4 . 

Step 4. 3^ is a solution to (18 . 3al >. for each to < T(C). 
Define the open cover (Wi , W 2 ) of W, 

Wi = w n (m 2 x (o, D) w 2 = wn(R 2 x(i,i)) 

If q G Wi, then S" to (g) = (see, Step 2) and Src(g) = (see, (S2) and (S4)). In this 
case, St is self-evidently a solution to d8.3at . On the other hand, if q G W2, then 
7/1(9) = 1 (see, (S9)). By direct inspection of (S10), the equation M(g, 5)S = S 1 ) 
collapses to J8.3al >. 

Remark 8.9. Recall from Remark 1731 that, for all (£ , b, i)el 2 x[l,2]x (-00, -±), 
the set 0(£o, &j t) is a bundle over the u-interval (0, b) whose fibers are disks centered 
at £o> with radii < |. 

Step 5. For each t < T(C), let S% be the solution to ( 18.3al > '. Let 

y(to) = [i ,T(C))x {|||| <|e|<||||}x (0,1) 

= (J U {r}xO(6,l,r) C y(to) 

re[to,T(C)) |f|<|«o|<2|f| 

Note that both S% and Flip .a • S^ a (see, Conventiori \8.9\) are defined on y(to). We 
claim that S? o = Flip a. • S^° on Z(t ) C y{t ). 

The argument is by finite speed of propagation. For any | |j- | < |£ 1 < 2| § |, let 

2"(&,*o) = 

{t G [*o,T(C)) J 5£ = Flip. • S t 7 on |J {r} x 0(£ O) l,r) C Z(t )} 

T€[t ,t] 

OurgoalisZ(£ ,to) = [to, T(C)), for each |f | < | Co I < 2|f |. First of all, T(£ 0j i„) is 
an interval that contains to, because S° a and S^* 7 have trivial initial data. By continuity, 
it is closed as a subset of [t , T(C)). We now show thatZ(£ , to) is also open as a subset 
of [t ,T(C)). 

By the second to last inequality in ( 18.431 ), and by ( 18.351 ), 

sup supg j T) {S£}(r) < § c 6 (R)c' 2 (B). 

re[t ,T(C)) 

Let t' G X(^o, to)- By the definition of Z(£o,to), and by continuity 

sup Supg> iT) {Flip ft -S-°}{t) < lc 6 {R) C ' 2 (B) 

T"G[t ,**) 
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for some t* e (£', T(C)). The assumptions of Proposition l8.3l Part 3, (Sys2) are satisfied 
with T = t*, = Ef o , 3^ = Flip f • 3^° ' . Conclusion 2 of Proposition [83 

enables us to apply Proposition [TT] (Refined Energy Estimate) for (Sys2), with J = \. 
The assumptions of Proposition l7.7l are satisfied because 



J = \ > Jo, c' 2 [B) < c 3 (X*), t* < T(C) < -l/c' 2 (B) < -l/c 3 (X*). 

(8.44) 

In the present case, the Refined Energy Estimate ( 17.42b becomes 



N 1/2 v / ^^n^) + 







(£o,i,t)1. j IV ) - -4v- ; | r | 1/2 

for all r S (to,t*), where T = Hf — Flip a. • Furthermore, the energy on 

the right hand side is zero. The vanishing of the energy on the left hand side implies 
[to, t*) C I(£o> to). Hence, T(£,o, to) is an open subset of [to, T(C)). 

Convention 8.11 (for the remaining steps). We introduce a new field for the remaining 
steps. For each to < T(C) and a £ { — ,+}, it is the map (see, Convention l8.9b 



[to, T(C)) x R 2 x (0, 1) -► ir^K 



fs&(<z) if|Cl<2||| 
\ Flip f --7(g) if|C|>i|f| 



It is well defined on the flip-invariant [i ,T(C)) x < |£| < 2|f |} x (0, 1), 

which is contained in £(to) U (Flip a • -Z(io)) by Step 5. It coincides with 3^ Q on the 
set Z(to) of Step 5. Consciously abusing notation, we will denote this new field by the 
same symbol 3° Q . 

Step 6. For each t < T(C) and a G { — , +}, 



"to 



Flip . • 37° on [to, T(C)) x (R 2 \ {0}) x (0, 1) 



The field 3Z ■ [to, T(C)) x R 2 x (0, 1) — > 7T 1 1Z is smooth, vanishes when q 3 < |, 
and extends, with its derivatives of all orders, continuously to [to, T(C)) x R 2 x [0, 1]. 
Moreover, it is a solution to J8.3al > CT on its entire domain of definition and satisfies 

sup sup Supg- 2) lT) {» t CT }(r) < ±c 6 (R)c' 2 (B) (8.46a) 

|fo|<2|f | re(to,T(C)) 

su P sup [r\^Su P ^\ r) {SZKr) < 2c *W^ (8.46b) 
|&|<a|f | re(u,.T(c)) kS0 ' ' ; c 7 (y) 

The main point before (I8.46at is that Flip a. is a field configuration symmetry. The in- 
equalities ( 18.46ab . ( I8.46bb are consequences of ( 18.43b and the Sobolev inequality (18.351 ). 
Be aware that that the 3 to in (18.43b is related to the present 5g by ( 18.45b . 

Step 7. For each to < T(C) — 1 and a £ { — , +}, 

t 1 

sup sup E o{( } {(3Q }(r) < ( y, Jo) --. 

ISol<2|f| re(t +l,T(C)) 1*0 1 
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where ^ s * ne constraint field associated to the field ( 18.451 l. 

For any |£o| < 2 ||M and any t* s (to + 1, T(C)), the assumptions of Proposition [83] 

Part 3, (Sys3) are satisfied with T = t*. By Conclusion 3, R > 4, A" > 0, and (18.46bb . 

sup |r|Sop^ 0jliT) {(SS)*}(r) < (y , Jo) 1. 

re(to + l,t*) 

By continuity, this holds for t = to + 1 as well. Therefore, the energy satisfies 

l*o + l| 2 i£ (eo , Mo+1) {(^)'}(*o + 1) < {YM 1. (8.47) 



By Conclusion 2, we can apply Proposition 17.71 (Refined Energy Estimate) for (Sys3), 
with J = |, X = (tp + 1, **). The assumptions of Proposition lTTl are satisfied, because 

J = \ > Jo, c' 2 (B) < c 3 (X), t* < T(C) < -l/c' 2 (B) < -l/c 3 (X). 
The Refined Energy Estimate ( 17.42b and ( 18.47b imply 

SU P £ o«o,i,r){( s S) K r ) ~PVo) r: t 71 £(y;j ) in- 

re(to+l,t*) 1*0 + 1-1 |*o I 



Step 8. For all a e {- , +} and all t x <t 2 < T(C), 



1 



SUp Sup £ «o.1,t){"£ -S£}(r) <(y,j ) t—t 

l?o|<2|f I re(t 2 ,T(C)) ' l*2| 

Foranyt* G (t%, T(C)) and |£o| < 2| ^ |. the assumptions in Proposition l8.31 Part 3, for 
(Sys2), are satisfied with t = h, T = t* and = S£, S (2) = S£. By Conclusion 
2, we can apply Proposition ^. 7l (Refined Energy Estimate) with J = |, see ( 18.44b . The 
Refined Energy Estimate ( 17.42b implies 

pO r-a gg -,/ v < / / 2 NjW^ if^K^) < 1 

for all re (<2, **)■ For the second inequality, see (I8.46bb . 

Step 9. There exist 

E a : (-oo, T(C)) x M 2 x (0, 1) tt^K, cr G {-, +} 

with S a = Flip a. • Ez~ a (see, Convention \8.9i , that vanish when q 3 < i, are C R ~ 3 , 
and extend, with their derivatives of all orders < R— 3, continuously to (— oo, T(C)) 



x 



[0, 1]. Moreover, they are solutions to both ( 18. Sab ' 7 a«c/ (^ cr )" = 0, anc/ 



sup sup Supg 3) 1t) {~ ct }(t) < ic 6 (i?)c' 2 (B) (8.48a) 

l?o|<2|f | re(-oo,T(C)) 

sup sup |r| K+1 Supg- 3) lT) {E CT }(T) < ( y )Jo) 1 (8.48b) 

|Co|<2|f| re(-oo,T(C)) 

sup |rr +1 Supg- 3 ; (0)x(0il) {» CT }(r) < (y , Jo) 1 (8.48c) 

re(-oo,T(C)) *' 'Br 1 v ' 
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For each (3 G (0, 1), introduce the compact set 

Xp = [T(C) - (3-\T(C) - p] X l> 2 |a|(0) x [0, 1] 

For every sequence t n — > — oo, with t„ < T(C) — the sequence of fields E^ is, 
by Step 8, a Cauchy sequence in L 2 {Xp). Set S"\ Xf} = L 2 - lim^-oo Sf . By ( |8.46at , 
the 1-parameter family E% with i < T(C) - is a bounded subset of C R_2 (Af(3), 
the space of 7r _1 7£ = R 31 valued functions of class C R ~~ 2 on the interior of Xp, that 
extend continuously, with their derivatives of all orders < R — 2, to the boundary. By 
Arzela-Ascoli, there is a subsequence that converges in C R ~ 3 (X@). Therefore, E a \x l3 
is in C R ~ 3 (Xp). It follows that S" is C R ~ 3 on the interior of U/3e(o l) ^Z 3 ' anc ^ extends 
with its derivatives of all orders < R— 3 continuously to U/3e(o i) ^ = (^°°;T(C))x 
D 2 |f |(0)x [0, 1]. By construction, = Flip a •»'- <T on ,fyn(Flip » -Afg). Hence, the 

pair of fields E a have unique C R ~ 3 Pole-Flip compatible extensions to (— oo, T(C)) x 
R 2 x (0, 1), which extend with their derivatives of all orders < R — 3 continuously to 
(-oo, T(C)) x R 2 x [0, 1], as required by Step 9. 

It follows directly from Step 6, that the pair E a has all the desired properties, including 
the bou nds ( I8.48ab ( I8.48bl i (recall that R - 3 > 1), with the exception of (E a f = 0, 
(18.48cb and (I8.48db . It is implicit in our construction that for each (3 G (0, 1), there 
is a sequence t n — > — oo so that E° — > S' ' in ^iXp), and therefore (S^ )" — > 
(S 17 ) 8 in C°(A^). Now, by step 7, (S^)*!^ = 0, for all /3 G (0, 1). By Pole-Flip 
compatibility, (E rT )^ = everywhere. The estimate ( I8.48ct follows from (I8.48bb when 
|£| < 2|f |. For < |£| < 4|||, it also follows from jgggp , by using Pole- 

Flip compatibility and Lemma IfTI in Appendix IF) To verify d8.48dl i. observe that the 
assumptions of Proposition 18.31 Part 3, for (Sys2), are satisfied for any |£ | < 2||j|, 
t < T < T(C), = E ^ = E° ' . By Conclusion 2, (RE1) holds in the context of 
(Sys2), which implies ( T8.48dt for a G (— oo, T(C)) X D 2 |f |(0) x (0, 1), and for general 
q by Pole-Flip compatibility. 

Step 10. The fields E a in Step 9 are unique in the following sense: Suppose, for some 
h < T(C), the ^-fields E a : (-oo,i x ) x R 2 x (0,1) -> -k~ x TI are Pole-Flip 
compatible, extend with their first derivatives continuously to (— oo,ti] x R 2 x [0, 1], 
are solutions to i&.3ab ' J . vanish when q 3 < i, and satisfy 

lim sup |r| Jo Supgl lT) {E CT }(r) = (8.49) 

— 00 |Co|<2|f | U(t °. ■ ' 

Then, E a = E a on (-oo, t x ) x R 2 x (0, 1). 

For every sufficiently negative r < ti, the assumptions of Proposition 18. 31 Part 3, for 
(Sys2), are satisfied for any |£ | < 2|f |, t < T = t, E^ = E° , E^ = E" . The 
condition that t is sufficiently negative is used to verify hypothesis (iv) of Part 3 of 
Proposition 18.31 By Conclusion 2, we can apply (checking, the additional hypothesis 
similarly to ( 18.44b ) Propositionl7T7lwith J = Jo- It follows from ( 17.421 that 

E o( io ,i,r){- j( T ) < ( c 4(^ )) rnsjo • 

We take the limit to — > — oo, keeping r fixed. By d8.48b| ) and ( 18.491 l, and the fact 
that 2 Jq < 1, we conclude that the energy Eq(£ \ T ){_^ a ~ ^ a }( T ) = 0. Exploiting 
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the Pole-Flip compatibility, £ , °'(t, • ) = S a (r, • ) for all sufficiently negative r. To 
demonstrate that S a = S a on (— oo, t\) xR 2 x (0, 1), we make a closed-open argument 
almost identical to the one in the proof of Step 5. 

We finally return from the q = (t, £, u) to the x = (£,u,u) coordinate system, and 
complete the proof of Theorem 18.11 The x-set Strip (l,c) is contained in the q-set 
(-oo, T(C)) x R 2 x (0, 1), by the smallness condition (839). 

Existence in Theorem\8J] follows from Step 9, with <P a = Hsfe + ir E° ' . We only have 
to check ( 18.331 . ( 18.341 and conditions (*) and (★*) (see, Definitions 12. 1 l and 12.2b . and 
apply Proposition ^. 21 Write 

!^(. > «)-!^(0)(-,«) = i £(i) fe <^(fe+l)(-)+7T »*(■,**) 

k=0 

The coefficient functions ^(k + 1), appearing on the right hand sides, are estimated 
using ||!f r0 '(fc + l)||c H + 1 (C(a,a,2)) 1 (see, the Overall Preliminaries in the proof of 
Proposition IOI and S a is estimated using ( I8.48ct . Now, ( I8.331 l follows. Also (18.34b 
follows, with an additional legitimate smallness condition on c depending only on 
(Y, Jo). Condition (*) is a consequence of the inequality \ < e% < 2 on Strip (l,c), 
which follows from (I8.48db . Here, is a component of 3 a = (e, 7, w) = M a M + 
u~ M ^ a . Finally, the equation L(e\~e2 — eie.2) = —2j2 — &1&2) (a consequence 
of the first two lines of ( |2.4t ) implies that 3(eie2) cannot change sign along the inte- 
gral curves of L, and therefore S(eiea) < on Strip(l, c) because it is negative when 
u<\. This implies (**). 

Uniqueness in Theorem \8. 1 I follows from Step 10. □ 
9. Conclusions 

Theorem 18. II is the central mathematical result of Sections [7] and [8] It establishes the 
existence and uniqueness of classical vacuum fields, under appropriate smallness con- 
ditions. In this section we extract a number of corollaries of Theorem l8.ll 

First, we show (Proposition 19. It that the far field formal power series is an asymp- 
totic expansion for the classical solution of Theorem l8.ll In this sense, the solution is 
quasi-explicit. Using this result, we can systematically exploit Theorem 18. II by trans- 
ferring properties of the formal solution to corresponding properties of the classical 
solution. 

The high point of llChrl is the demonstration that strongly focused gravitational 
waves generate trapped spheres. We recover this result (Proposition 19. 21 i. One would 
like to show that the trapped sphere is just a shadow of a black hole, that is, a hori- 
zon forms. We take a first step by showing (Proposition 19.51 that for special initial 
conditions, the solutions of Theorem l8. 1 I become arbitrarily close to the Schwarzschild 
solution on the upper edge of the strip, between the trapped surface and past null infin- 
ity. This would be used in conjunction with a controlled perturbation expansion around 
the Schwarzschild /Kerr family to construct the global exterior of a black hole. 

An important feature of our approach is that the scaling limit a = 21 J, is reg- 
ular, see Remark 1531 Informally, "(SHS) plus = 0" (see, d5.7ab ) degenerates to a 
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2-dimensional system, separately for each £. Its solutions break down in finite time. We 
speculate in Subsection [93] that an appropriate expansion in 21 around this two dimen- 
sional system can be used to explore the spacetime singularity. 

We also obtain, in Subsection 19.61 a class of solutions distinct from that of II CM . 
corresponding to more general initial data at past null infinity. This class arises by from 
the limit 21 J, with a > fixed. 

We assume, without further comment, the definitions and conventions of Section [2] 
through Section|6]and Theorem l8.ll However, see the Index of Notation, AppendixlAl 

Proposition 9.1 (Asymptotic expansion). Let <P a = M a ,<& + u~ M *I/' y be the pair of 
Pole-Flip compatible vacuum fields of Theorem \8. l\f or K = 0. For each L > 0, 



sup 

\a\<R-3 



E 



k=0 



C°(c(a,8t,l)) 



o 



\L+1 



as u — > — oo. In other words, the far field formal power series [W 7 ] is an asymptotic 
expansion for 

Proof. Observe that the conditions imposed on the data a, 21, DATA CT , R, e in Theorem 
18. H are independent of K. Therefore, Theorem [8T| can be applied with the same data 
for all K > 0. For each K > 0, we obtain a pair of Pole-Flip compatible vacuum 
fields on Strip(l, \k) C M 4 , where \k > depends on K through the vector C in 
Theorem 18. II The vacuum fields corresponding to any pair K, K' > coincide for 
sufficiently negative u, by the uniqueness statement of Theorem l8.ll In particular, this 
is true for K = 0, K' = L. The bound ( 18.34b for the K' = L vacuum field implies the 
proposition. □ 



9.1. Three Points of View. It is helpful to consider the focusing of gravitational waves 
from three perspectives, that yield three different pictures. 

Regularized Picture (R) 
t 

Field transformation £ o 21 
where = a£ 

High Amplitude Picture (H) 

Isotropic scaling transformation -3 
^ with scaling constant 3 = 2l 4 

Finite Mass Picture (F) 



First, recall from Section[3]that the Isotropic Scaling J, the Anisotropic Scaling 21 and 
the Angular Coordinate Transformation £ are field symmetries (see, Definition 13.1b . 
Their isotropic respectively anisotropic character refers to their action on the frame and 
the coordinate system. Both scalings are, at the level of the Lorentzian metric, global 
conformal transformations. 

The pictures are fixed by the table 
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Regularized 


High Amplitude 


Finite Mass 


Background 


M a ,% 






Data 




a _ V(f£> u) 


2TV(f£, 2l- 4 u) 


Domain 


Strip (l,c) 


Strip(l,c2l 2 ) 


Stri P (2l 4 ,c2l- 2 ) 


Hemisphere 


Kl<lfl 


l£l<i 


iei<i 



The Regularized Picture is the arena of Theorem 18. II The other two are obtained, as 
indicated above, by scaling. Row 1 displays the Minkowski background fields. Row 2 
gives the functional form of the initial data at past null infinity. More precisely, 

data£(£,u) = if(^u) = M a , m + u- M &% (9.1a) 

data£(£,w) = 2TV(f£, u) $ a H = M 1A + u- M ^ (9.1b) 

data£(£,u) = STV(££, 2T 4 u) $f = M hl +u- M ^ (9.1c) 

The first equation gives two new names for the DATA CT (^, u) of Theorem 18. II Row 3 
displays the functional dependence of the domains on the scaling parameter 21, with the 
notation Strip(/i, A) = M 2 x (0, x (— oo, — A -1 ). Row 4 gives the size of the £-disk 
which, under stereographic projection, corresponds to one hemisphere of S 2 . 



9.2. Two Physical Regimes. There are two natural physical regimes. Informally, both 
appear as limits 21 J, in the Regularized Picture, keeping the data rf fixed. They are 
distinguished by: 

• 2D (Scaling) Limit: a — 21 J, 0. 

• 4D (Scaling) Limit: a fixed, 21 J, 0. 

(The 4D Limits breaks Pole-Flip compatibility of rf . This will be discussed below.) 

Definition 9.1. The 2D Limit Assumptions are the hypothesis, with a = 21 and e = j, 
of Theorem \8. 1 1 on 21, DATA' 7 = rf , R and K up to and including condition d8.32| l. 

Remark 9.1. Explicitly, the 2D Limit Assumptions are: R > 4, K > 0, rf = when 
u < i, and 

0<|2l|<b, rf = Fli Pl ■ ?f a , max ||?f|| C B + 4 (D4(0)x(0i2)) <b 

cr6{-, + } 

Here, b £ (0, 1) depends only on R. The constant c £ (0, 1) in Theorem 18.11 depends 
only on R, K and max CT(E {_ +} |K||c«+ 2A "+e(D 4 (o)x(o,2))- 

The conditions on rf : M 2 x (0, oo) — ► C are independent of 21. Also the domain of 
definition Strip ( 1 , c) of the vacuum field in Theorem l8.1l is independent of 21. Therefore, 
Theorem l8.1l is consistent with the 2D Limit. The choice e = j is just for concreteness. 

Remark 9.2. The intuition behind the designations Regularized Picture and High Am- 
plitude Picture is immediately clear in the context of the 2D Limit. For the Regularized 
picture, see Remark 1631 In the High Amplitude Picture, the initial data at past null in- 
finity for the corresponding family of vacuum fields is DATA^ (£, u) = 2l~ 2 rf (£, u). It 
grows unboundedly as 21 J, 0. The Finite Mass Picture will be discussed momentarily. 
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Remark 9.3. From our perspective, MChrl investigates the 2D Limit (a = 21) in the 
Finite Mass Picture. Christodoulou's small parameter S > is to be identified with 
our 2l 4 . With this translation, the first equation in d9.1c| ) is precisely Christodoulou's 
"short pulse ansatz". For MChrL the "short pulse hierarchy" plays a central role (see 
equation (24) on page 20 in IChrl . and the following discussion). In our approach, 
this hierarchy plays no role at all. However, it can be recovered through the scaling 
transformations required to go from the Regularized Picture to the Finite Mass Picture, 
see J9.5cb and d9.6c| ) below. By contrast, our working picture, the Regularized Picture, 
merely contains a dichotomy: the ^P-even components display one behavior, the ?P-odd 
components another, see Remark 16.31 or ( 19.81 ). This dichotomy disappears in our 4D 
Limit. 



9.3. Trapped Spheres. 

Proposition 9.2. Make the 2D Limit Assumptions (Definitional} with K = 0. Set 

du\v a (Z,u)\ 2 



A(u-\ ) = min inf 

<ye{-,+} ee-D 4 (o) J 

Suppose Aiu^) > Ofor some ^ G (0, 1). Then ^ an d 76 are everywhere negative on 
the sphere 

Su.u ■ (u, = - \ A{u 1 ) 2T 2 ) 

whenever 21 G (0, b) is sufficiently small depending only on A(v^) and c. For instance, 
21 < Jq/c min{l, A(u])} will do. In other words, Su,u is a trapped sphere (see, Re- 
mark \23i . 

Remark 9.4. Clearly, there is an infinite dimensional family of pairs rf , 21 satisfying the 
assumptions of Proposition ^. 21 

Proof. By (18.34b . the components 7f , 7g of $ a = M%,% + u~ M & a satisfy: 



2„,2 



2l 2 u — u 

2 



U 



% 2 u — u 



< 



i 



< 



c\u\ 

1 

c \u\ 



for all (£, u, u) G Ai(0) x (0, 1) x (— oo, — c 1 ), where (see, Definition [5JJ 



wj(o)(e,«) 



<(0)(^M)=0 



To find a trapped sphere, let u = — A21 2 , where A > 0. Now, for all (£, u) G £>4(0) x 
(0, 1) and cA > 2l 2 : 



21 2 

If A- / ds|7 ? ' T (^s)| 2 + — < , then u, -A2T 2 ) < 0. 



If 



A 2 



cA 
21 



(9.2) 



2l 2 (2l 4 + A) cA 



— < , then 76(C,m,-A21" 2 ) < 0. 



Proposition l9.2l is a direct consequence of (E3 with A = i^(w!),if we also recall that 
Flip x does not change the sign of 73 and 7g , □ 
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9.4. The 2D Limit in the Finite Mass Picture for a Finite Duration Pulse. We make the 
2D Limit Assumptions (Definition ^. It . and the assumptions for a finite duration pulse: 

rf : R 2 x (0, oo) -> C has support contained in R 2 x (±, |). (9.3a) 

,3/4 

/ durf (£, u) = for all £ and cr. (9.3b) 

A/2 

Remark 9.5. Clearly, there is an infinite dimensional family of Pole-Flip compatible rf 
satisfying (|9.31 l. The conditions ( 19.31 ) are equivalent to the condition that rf {£,u) = 
^G a (^,u) for a function G a : R 2 x (0,oo) -> C with supp G ff C M 2 x (±, |). 

Remark 9.6. One can relax the conditions (19.3b by imposing suitable smallness condi- 
tions for rf on M 2 x (§,oo) and for du^K, u), depending on 21. In this case, 
contrary to our current convention, rf itself would have to be allowed to depend on 
21, say polynomially. A smallness condition depending on 21 is an "open" condition, as 
opposed to ( 19.3l l. 

We begin with an informal discussion of the 2D Limit a = 21 J, in the Finite Mass 
Picture. Fix u$ < 0. The shaded region in the figure (just below) is Strip(2l 4 , luol" 1 )- 
For |2l| small enough, it is contained in Strip(2l 4 , c2l~ 2 ), the domain, in the Finite Mass 
Picture, on which the solution to Theorem 18. 1 1 exists. Region I is Minkowski space 
A^i^. In the "pulse region" II, the solution grows unboundedly as 21 J, 0. Nevertheless, 
as we will show, ( 19. 3t implies that on each compact K C Strip(oo, |uo| _1 ) the far field 
formal solution converges as 21 [ 0. Note that as 21 J. 0, the compact K is eventually 
contained in the "after the pulse region" III = Illi U III2. One would like to have an 
analogous result for classical solutions. In this paper we take a step towards such a 
result, by controlling the solution in HIi, a strip moving and shrinking with 21. 




If ^ du |?7 cr ('C, u)\ 2 is positive and independent of £ and a, the limit 21 J, of the 
formal power series solution is the field corresponding to a Schwarzschild spacetime, 
whose future horizon is a level set of u, with u < uq, when |uo| > is sufficiently 
small. 
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We use the notation 

*Kcture ™d [ ££ cture ] = ^ ^Hcture(^) (£. 3*) 

for the field and for the formal power series, where Picture = R, H, F. Let RH, HR, 
FH, HF, FR, RF be the "transition matrices" between the different pictures, given by 

RH = — HR = diag (3, 3, 4, 4, 4, 2,4,3,3,3,2,2,4, 2, 3,4,5,4) (9.4a) 
FH = — HF = diag (4, 4, 8, 8, 8, 0,4,4,4,4,4,4,8,-4, 0,4,8,8) (9.4b) 
FR = -RF = diag (1, 1, 4, 4, 4, -2, 0, 1, 1, 1, 2, 2, 4, -6, -3, 0, 3, 4) (9.4c) 

Observe that FR = FH + HR. Then (use: the figure in Subsection |9.U equations ( 19.11 ), 
a = 21 and Definitions Il2l Il4l [33T > 

= 2l HR ^(£,w,2l 2 U ) (9.5a) 
^{^u,u) = 21™ !Pg (f, 2T 4 u, 2T 4 u) (9.5b) 
= 2l FR ^(£,2T 4 ;u,2n 2 u) (9.5c) 

The coefficient functions of the formal power series transform according to 

*S (*)&«) 

For all fc > 0, we have: 

(k){£,u) is a polynomial in 21 (9.7a) 
W (£i ^) i s a polynomial in 2l -2 without constant term (9.7b) 

Statement ( 19.71 ) is verified by induction over k. Just follow the construction of [3^ ] 
and [<?h ] in the proof of Lemma |6"T1 keeping in mind that the coefficient functions, in 
the Regularized Picture, of the Minkowski background [ A4a,a ] depend polynomially 
on 21, and that DATA^ is independent of 21. On the other hand, in the High Amplitude 
Picture, the Minkowski background [.Mi,]. ] is independent of 21, while DATA^ is pro- 
portional to 2l~ 2 . Incidentally, ( 19.7al ) has already been shown in Remark |6*31 There, it 
was also shown that 

*P-even (*}3-odd) components of \P£ (fc) (£, u) are even (odd) polynomials in 21 (9.8) 

The statement d9.8l ) also follows from d9.7a| ), d9.7b| ) and d9.6at , because the i-th com- 
ponent of i^(fc)(£, m) is ^5-even (^-odd) if (RH)ii is even (odd). Furthermore, 

the i-th component of u) has degree < (RH)^ - 2 + 2fc (9.9) 

as a polynomial in 21. 



= 2T R - 2fe ^(fc)(<e,u) (9.6a) 
= 2l FFl+4fc ^(fc)(£,2T 4 u) (9.6b) 
= 2l FR+2fe ^(fc)(£,2TV) (9.6c) 
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Lemma 9.1. The 18 components of each &p(k) (£, 2l 4 w), k > 0, are Laurent polynomi- 
als in 2l 2 . /fa component does not appear on the list 

M£(fc)(£,2l 4 u) fe = 0,l 
^ 2 )?(0)(£,2l 4 ^) 

(zi)£(fc)(£,2t 4 u) fc = 0,l,2,3 (9.10) 

(^ 2 )£(fc)(£,2l 4 w) fc = 0,l 

(z 3 )?(0)(e,2l 4 «) 

f/ie« /f is an actual polynomial in 2l 2 without constant term. 

Proof. By d9.4bt . ( |9.6bt . ( |9.7bt . all the components are Laurent polynomials in 2l 2 . By 
d9.4c| ). d9.6c| ). d9.7a| i. all those not in ( 19. 1 Ob are polynomials without constant term. □ 

Lemma 9.2. Suppose <|93]). Define U = R 2 x (|, co) andA a (£_) = f*J* du\if(^ u)\ 2 . 
For (£, m) G U, the functions in ( 19. 10b are polynomials in 2l 2 , and: 

( Wl )?(0)K,2l 4 1 i)=0 ( Zl )-( 2 )(^2l 4 2i )=2(e lil |-A lil )e lil |^ 

( Wl )?(l)(e,2t 4 u) = (zi)£(3)(£,2l 4 n) = G(2l 2 ) 

M£(0)(£,2l 4 u) = (z 2 )?(0)(^,2l 4 «) = (9.11) 

( Zl )£(0)(£,2l 4 «) = (z 2 )?(l)(^a 4 u) = -2e lfl |^ 

M(l)(£,2l 4 u) = (z 3 )?(0)(£,2l 4 u) = 

where 2^ - ^ + and e hl = |(1 + |£| 2 ) and Ai,i = -±( s fl + z£ 2 ). 

Proof. By d9.6ct . the equations d9.1 It can be translated from the Finite Mass to the 
Regularized Picture. In this proof, we work exclusively in the Regularized Picture. For 
convenience, we suppress the R and a indices as well as the argument (£, u), For ex- 
ample, u>i(k) means (wi)£(A;)(£, w). We use the shorthands e = e^.a and A = Aa,a 
(see, d4.3b "). Equivalently, e = 2lei i and A = S3tAi,i. For all the equations in d9.1 U 
concerning zeroth order coefficient functions, use Definition 15.11 and d9.6c| ). We only 
note that on IR 2 x (0, oo), 

z 3 (0) = -4(e^ + A)(e 1 | + 2A)a„ 1 r / -7 ? ^ 1 7/ + a„ 1 |77| 2 
which reduces to z 3 (0) = A" on U. For the rest of (19. lit , we use the equations 

N( Zl ) + ±D{z 2 ) = ^(Szt - 2Xz 2 - uazj) 

+ -±3 (2SXz 2 - Sluxz 3 + 6lu 3 z 2 + 4 : lu 8 z 1 - 4ZJ 4 z 2 ) 
N(z 2 ) + ±D(z 3 ) = + 25z 2 - 2cj 6 z 2 ) 

+ 4^ ( - 2cjiz 4 + 3w 3 z 3 + 2uj s z 2 - 3aJ4Z 3 ) 
L(uji) = -zx + \2% 2 lo 1 + 4j( ~2S% 2 uj 1 - 2uj 1 uj 2 ) 

The first and third appear (I5.7ab . For the second, we use ( I5.7ab and the constraint equa- 
tion 2/1 = 0. For the vector fields D, N, L, see d5.9t . We obtain, in succession, 
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on R 2 x (0, oo) 



onU 



HQ.) 
*a(l) 
wi(l) 




= 








Zl (2) 
Zl (3) 



= (not needed) 




(not needed) 



* 2 uzi(2)-§(e£-A)z 2 (2) 



For Zi(3), we have also used that 0)3(0), 0)4(0), /i(0), /2(0) all vanish on U (see, 
Definition l5.1b . We know from ( 19.8b that z\ (3) is a polynomial in 2l 2 . It has no constant 
term on U, because e, A and 2:2(2) are odd polynomials in 21. Now, use J9.6cb . □ 

Lemma 9.3. Suppose ( 19. 3b . Each component of each u), k > 0, is a Laurent 

polynomial in 2l 2 and a polynomial in u when u > |2l 4 . 

Proof. By d9.3ab and the construction of the formal power series in the proof of Lemma 
16.11 each ^(k)(^,u) is a polynomial in u on U = K 2 x (|, 00). Also recall (|9.7at and 
d9.8b . The lemma now follows from d9.6cb . □ 

Let (k)(£,u) be the polynomial extension in u of &j?(k)(£,u) from u > |2l 4 to 
u > 0. Then, [ W£ ] is a Pole-Flip compatible pair of formal solutions to (I5.7al ) on 
Strip ^ with Minkowski background [M^i ), and * ] =0. 

Lemma 9.4. Suppose d93T >. 77?en, (0)(£, w) is a polynomial in 2l 2 . More precisely, 



(0,0,-uA a ,0,0, 0,-A a ,0, 0,0, 0,0,uyT, 0, 0, yT, 2 m u ^ 0) + 0(2l 2 ) 



Proof. By direct calculation. □ 

Proposition 9.3. Suppose ( 19.31 ). £ac/i component of each iPf (k)(£,u), k > 0, is simul- 
taneously a polynomial in u and H 2 , for all € R 2 x (0,oo). 

Proof. They are polynomials in u by definition. The case k = is covered by Lemma 
19.41 The general case is shown by induction over k, using the fact that the equations 
(16.2b hold with Minkowski background [ M. 1,1 ] . In the present case, ( 16.2b are equations 
for polynomials in u. The generic terms Vk on the right hand sides in ( 16.2b are, by 
the inductive hypothesis, polynomials in 2l 2 . When using d6.2ab through d6.2r| ) in this 
order to determine the components of 'Pp(k), only polynomials in 2t 2 are generated. 
If ^ appears on the left hand side, then the non-constant terms as a u-polynomial of 

the corresponding component of Wp (k) are determined uniquely by the right hand side. 
The constant term of integration is determined by the restriction of fy F to u = 2l 4 , that 
is $ F (k)(^, 2l 4 ), which is itself a polynomial in 2l 2 , by Lemmas |9~T1 and l9~2l □ 

Proposition 9.4. Suppose d9.3b . For each k > 0, let &F<&=o(k) be the constant term of 
^(k) as a polynomial in 2l 2 . Then [tf^ a=0 ] is the unique formal solution to d5.7ab 



*£(o)(e,a) 



(9.12) 
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with Minkowski background [Mi.i] and characteristic initial data 

*F,a=o(0)(£, Si) = 

(0,0,-^,0,0, 0,-^,0,0,0,0,0,^, 0,0, ^^uei.i 4vL CT ,0) (9.13a) 



[#?,«=„ ](£,0,u) = 



1 






© 



-A' 




V o / 



V 








J 



(9.13b) 



Its coefficient functions are polynomials in u. Moreover, [^ a=0 ] = ^ - 

Particularly, if A a (t;) = A is independent of £ and a, then it represents Schwarzschild 

spacetime, with mass m = 2~ 3 / 2 A 



Proof. The first part follows from Lemma |9TT| 19. 2119. 3119. 4| and Proposition ^. 31 In par- 
ticular, d9TT3ab follows from (|9~T2l , and d9.13bt follows from d97TTb . 

If A a {£) = A is inde pendent of £ and a, then ^ a=0 ( )(C> u) and [#g a ](g, 0, u) 
in d9.13at and (19. 1 3bb correspond to spherically symmetric initial data for J5.7al > with 
Minkowski background [A^i.i]. Therefore, [f F " a=0 ] is spherically symmetric, and 
therefore, by a formal Birkhoff Theorem, a formal expansion of a Schwarzschild vac- 
uum field. The component 72 of [A^i,i] + u~ M [{'p^J vanishes on the sphere 
(u,u) — {—A, 0), and therefore this sphere is a section of the Schwarzschild horizon. 
Its area is equal to 2-kA 2 , which gives the formula for the mass to. □ 



The discussion of the formal solution f \P£ 



tions. The bound ( 18.341 implies that for all u < 

K 



is finished. We now turn to classical solu- 

-1 



sup 

|a|<_R-3 



k=Q 



C0(D 4 (0)x(|,l)) 



< 



c \u 



K+l 



This bound, in turn, implies the Finite Mass Picture bound (use d9.5c| ), d9.6c| >) 



sup 

QGN* 
\a\<R-3 



K 

^f(-,«)-E 

k=0 



c°(n 4 (o)x(|2i 4 ,a 4 )) 



< 



ty2K-iR+& 



\K+\ 



when u < — 2l 2 c 1 . The power of 21 on the right hand side arises as 2K — 4R + 8 = 
2(K + 1) - 4(R - 3) - 6. Given R > 4, we choose K = 2R - 3. (Then, the constant 
c depends only on R and max ff£ {_ + } ll 7 ? cr ||c 5fi (n4(o)x(o,2))-) Altogether, we obtain: 

Proposition 9.5. Suppose i9.3i . For each uq < and each R > 4, the limit as 21 J, of 



sup \u 

u<ua 



2R-2 



sup 

\a\<R-3 



()' 



2 |^ 3 %(fc)(.) 

fc=0 



co(r> 4 (o)x(ja 4 ,a 4 )) 



is zero. Here, the solution the functions ^(k), and the u-interval (|2l 4 ,2l 4 ) de- 



pend on 21. Under appropriate conditions ( see, Proposition ^ A\ , the Schwarzschild vac- 
uum field can be approximated arbitrarily closely on the strip IIIi. 
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Remark 9. 7. To obtain the last result, we had to explicitly calculate the first four orders 
of the far field expansion, in particular for the component z\. 

So far, we have provided complete, detailed arguments for each of our statements. At 
this point of the paper, the character of our discussion changes. For the rest of Section|9j 
we sketch additional applications of our overall hybrid method and give informal argu- 
ments to support our informal assertions. We will give rigorous discussions in another 
place. 



9.5. 2D Limit in the Regularized Picture: beyond the far field regime. Theorem 18. II 
produces vacuum fields on Strip(l, c), with c > independent of 21 (see, Remark l9~Tl ). 
In this subsection, we informally answer the question: How can one control these vac- 
uum fields when u > — - and get closer to the (expected) singularity? The strategy, as 
before, is to generate an appropriate formal solution to the initial value problem (see, 
Section 0, and then to construct a classical solution by estimating its deviation from a 
truncation of the formal solution. The far field expansion has run its course. We need 
an expansion in 2l 2 . 

In this subsection, we make the 2D Limit Assumptions (see, Definition 19. It , work 
exclusively in the Regularized Picture, and suppress the R and a indices. For example, 
<P means <£> a R . 

In analogy with (I5.5ab . (I5.5cb . ( 15.6al i. ( I5.6H . set 

P = diag(l, 1, 0, 2, 2) © diag(0, 0, 1, 1, 1, 0, 0, 0) © diag(0, 1, 0, 1, 0) (9.14a) 
= diag(2,2,l,l,l)©diag(l,0,0,0,0,0,l,l,l)ediag(l,0,l) (9.14b) 
<P = 2l p (2T p 0) (9.14c) 

<P* = 2l p3 (%- p) ¥) (9.14d) 

The third line indicates that we wish to write <P as 21 p times a new field. In order not 
to introduce yet another name, we write the new field as 2l _p <£. Similar for 2l -P "^>". 
We make formal expansions of 2l~ p ^ in powers of 2l 2 . The properties of the far field 
expansion, for instance ( I9.7ab and ( 19.81 . suggest that the ansatz ( 19.141 ) is consistent. 
However, this must be checked. 

A formal power series { / } in 2l 2 on an open subset U C Strip^ with values in a 
vector space X is a formal sum 

oo 

{/} = E (« 2 )' /{*}(*) (9 - 15) 

e=o 

For each I > 0, the coefficient f{£} :U — * X is smooth and independent of 21. 

Let {2t~ p .Ma,a } be the formal expansion of 2l -P .Ma,» in powers of 2l 2 (see, 



Definition |4JJ|. It is defined on Strip^ and takes values in 1Z. 

Our ansatz is to write the field 2l~ p <£ as a formal series { 2t _p ^ } on some open set 
U C Strip^ with values in 7Z. In this context, the far field ansatz (Section [5} becomes 

{2T P 0} = { 2T p .Ma,a } + U - M { 2r p !f } (9.16) 

To define the associated formal constraint field, see Definition 12.41 and (I5.6bb . we fix 
the weight functions Ai,A2,A3,A4 by (I5.6cb . as before. Then, by direct inspection, 
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{ 2T P <2>« } or, equivalently, { 21" p ^ } are also formal power series of the form 
( 19. 1 5b . They are defined on U and take values in 1Z. For each i > 0, the coefficient 
(%- pi &){l} is determined by (2T p 0){m}, < m < I. Similar for (%r p *&$){£}. 

We want to formally solve the same characteristic initial value problem as before: 

• (EZD with # and Ma,si replaced by 2l p { 21 } and 2l p { 2T p .M a , a }, 
. {2T P V«}=0, 

• formal asymptotic initial conditions (I5.3ab . d5.3bb with, say, u = i. 

It is understood that DATA = 77 (see, Definition l5.ll ) is fixed and independent of 21. 

Remark 9.8. The meaning of the formal asymptotic initial condition (1 5 . 3 ab is 

00 

^(2l 2 )' lim (2r p ^)m(^,«) = 2T P ^(0)(£,«). 

£=0 

Equality is in the sense of TvL-valued formal power series in 2l 2 . The right hand side is 
actually a polynomial in 2l 2 of degree 1. The function ^(0) is given in Definition 15. II 

Constructing a solution { 2l~ p S r } to this formal initial value problem requires solv- 
ing an infinite family of differential equations. All but a finite number of them are linear. 
It is possible to arrange these equations so that, when they are solved step by step, the 
"angular derivatives" are only applied to functions that have already been 

constructed. An essential ingredient is 

D = 0{%) N=£+0(V1 2 ) L=(l + ^f 3 )£ (9.17) 

as 21 -> 0. Here, C(2l fc ) stands for 0(21*0^ + C(2t fe )^, when k = 1,2. In this 
sense, one only has to solve 2-dimensional problems in the (u, u) plane. 

Observe that property ( 19.9b of the - expansion implies that the coefficient function 
(2l _p tf r ){€}, £ > 1 is of the order 0(|it|~ £+1 ) as u — ► —00. For this reason, we expect 
that all the arguments in Section[8]can be applied, with minor modifications, when the 
function 'Pk in (S2) is replaced by the truncation 

K+2 

2i p J2 (a 2 ) £ (2l- p ^){£} 

e=o 

of 2l p { 2l _p if - }. One should be able to conclude, in analogy with Theorem 18. II that 
both a classical solution ^ and a formal power series solution 2l p { 2l _p !f - } exist on 
Strip(l, c), and that 

A'+l 

2T p <f - (^Y (2T P #)U}, (9-18) 

and all its partial derivatives up to some finite order, are estimated, in absolute value, 
by < c _1 2l 2K+4 |it| _A_1 on Strip(l, c). Here, smallness conditions similar to those in 
Theorem l8.1l must be made. In particular, c > has to be sufficiently small. 

The fact that the difference J9. 1 8b goes to zero as 21 J. 0, uniformly on Strip(l, c), 
means that the formal 2l 2 expansion "has not yet been exhausted". To better understand 
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what happens, let us examine the formal expansion in just a little more detail. We only 
discuss the zeroth coefficient, (2l _i V){0} or, equivalently, (2l~ p ^){0}, see (19.16b . 
Set a = 72{0}/e3{0} and b = 76{0}. The constraint equations U2{0} = it3{0} = 
u 6 {0} = and the equation ^e 3 {0} = 2e 3 {0} 3?7s{0} derived from d5.7ab yield the 
system S-a = —2ab and S-b = —2ab. The initial conditions (see, Remark |9~8l i 





e 3 {0} 


uj 2 {0} = u 2 72 {0} 


^6{0} = U 2 ( 76 {0}-i) 




1 








u — > — oo 


1 


-V|DATA| 2 






select the unique solution 



= m& h 7e{0} - 27T & h(Z,u,u)=u 2 -{<p(t,u)) 2 (9.19) 



where (p > and ip 2 — 2 d u 1 d u 1 |DATA| 2 . The solution ( 19.19b is defined on Uo, where 



Ue = {(£,«,«) e Stri Poo | u < -e- <p(£,u)} 



for every e > 0. The solution (19.19b and therefore the formal solution { 2t _p ^ } break 
down at u = —<p{i,u) < 0, for example because 76{0} diverges. Conversely, it can 
be shown that the whole formal solution { 2l -P <£ } to the initial value problem exists 
on Uo, that is, no "earlier" breakdown occurs. At the formal level, the scalar curvature 
invariants are in general unbounded as u | —<p(£,, u) < 0. We refer to u = — u) < 
as the formal (naive) singularity. Observe that 72(0} < and 76{0} < on Uo. 

It follows from the structure of the matrix P, in particular its nonzero entries, that the 
components e3{0}, 7i{0}, 7 2 {0}, 76{0}, 77{0}, 7s{0}, u>i{0}, W3{0}, ws{0} of the 
coefficient function (2t~ p <£){0} satisfy the quasilinear symmetric hyperbolic system 
and the constraints in Proposition 12.41 This system has been investigated in situations 
with higher symmetry, for example in [Sze |. In our present context, however, the fields 
depend on all four coordinates. The collapse of the frame as 21 — * 0, see (19.17b . is 
responsible for reducing the four-dimensional system to a family of two-dimensional 
systems, one for each £. It is possible to quasi-explicitly solve these two-dimensional 
systems near the formal singularity. That is, there is a formal solution given by an appro- 
priate expansion in the "distance" from the formal singularity which is an asymptotic 
expansion to the true classical solution (of the two-dimensional system). The behavior 
of the solution to this two-dimensional system leads us to speculate that the 2l 2 expan- 
sion exhibits an instability close to the formal singularity. This instability appears to 
drive the full four dimensional system into a new regime in which the classical vacuum 
solution may display features of the BKL scenario. See, [BKL| and references therein. 



We conclude this subsection with a discussion of the figure 
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O 



tp(£, «) formal (naive) singularity 



* // 



• Christodoulou HChrl constructs strongly focused gravitational wave solutions on I 
(see, Remark l9~3l l. Recall that u ~ — ^ is the place where trapped spheres first 
form, see Proposition l9.2l 

• In this paper, the far field expansion has been used to construct vacuum fields on the 
larger I U II = Strip(l, c), where c > is sufficiently small. See, Theorem l8.ll 

• The 2l 2 expansion outlined in this subsection allows one, using appropriate energy 
estimates, to construct classical vacuum fields on at least I U II U III = Strip ( 1 , e'~ 1 ) n 
U e . Here, e, e' > are arbitrary constants (in the figure, < e < e'), and |2l| is 
sufficiently small, depending on e, e'. Moreover, the formal 21 2 power series solution 
is an asymptotic expansion to the true classical solution as 21 — > 0, uniformly on 
I U II U III. In other words, the 2l 2 expansion allows one to construct and control the 
solution up to any "finite distance" from the formal singularity. 

The justification of the last statement relies on the fact that a suitable truncation of 
the 2l 2 expansion is an approximate vacuum field, with error terms going to zero 
uniformly on I U II U III as 21 — > (by a compactness argument). Furthermore, 
these error terms decay quickly enough as u — > — oo to be "integrable". See, the 
discussion of ( 19.181 ). 

Remark 9.9. Motivated by IBKI and IIAnReL we expect that coupling the gravitational 
field to a massless scalar field will make it possible to construct, under suitable generic 
conditions, strongly focused solutions from past null infinity all the way into a piece of 
the singularity. 

Remark 9.10. Observe that: 

• One only has to solve linear equations to inductively construct the far field expan- 
sion [ if - ] in Lemma |6TI 

• By contrast, the construction of the leading term of the 2l 2 expansion { 2l~ p !? , } 
requires the solution of a nonlinear (effectively two-dimensional) system. In other 
words, one is expanding around a non-trivial "background". 



9.6. 4D Limit. Recall that in the 2D Limit, the ^-odd components of the coefficients 
of the far field formal solution [<1>r] go to zero as 21 — ► (in the Regularized Picture). 
On the other hand, the ^-even components do not in general go to zero. See, Remark 
16.31 or (19.81 1. This asymmetric feature of the 2D Limit disappears in the more general 
4D Limit (see, Subsection |9.2t , as one can see by looking at the far field expansion. It 



76 



is for the purpose of taking the 4D Limit that the two scaling parameters, a and 21, have 
been carried along through the whole paper. 

To compare the 4D Limit with the 2D Limit, it is useful to formulate the second 
smallness condition in (18.32b in a picture in which the Minkowski background and the 
stereographic coordinates £ are fixed (independent of a and 21), for example the High 
Amplitude Picture. The smallness condition becomes 



^max^ ||^DATA^|| co(Dj(0)x(01)) < ^ 



OLI+OL2 



for all a = (ai, a 2 , a 3 ) € Ng with \a\ < R + 4. Notice that in the 4D Limit, |^| 
large factor, and there is one factor for each "angular derivative". 



A. Index of Notation 

This is a partial list of symbols used in this paper. It refers to their main / typical usage. 
Warning: These symbols can have different meanings. However, these meanings will 
always be made clear in each particular context. (For example, the entry for the symbol 
S in the index below refers to equation (14.41 ), which corresponds to the main /typical 
meaning of the symbol S in many sections of this paper. Nevertheless, the symbol S 
stands for a field transformation in Section[3] and it stands for a set in the local context 
of Proposition l7. II ) Symbols which only appear in the Appendices are not listed. In the 
third column, a selected reference is given. 



Symbol 


Typical Meaning 


See 


(*), (**) 


frame nondegeneracy conditions 


Definition |2JJ 


[A, B]= AB- BA 


commutator of operators 




<p 


parameter p dependent bound 


Conventionl7.1l 


21, a 


scaling parameters 


Definitions EUgT] 


si, a, 3, £ 


field transformations 


Section[3] 


A(<2>), f(<P) 


constituents of (SHS) 


Definition^. 31 


A(&,x), f(&, x) 


constituents of (I5.7ab 


Propositionl5.1l 


A(#), t{$,d x $) 


constituents of (SHS) 


Proposition^. 31 




constituents of d5.7bb 


Propositionl5.1l 


BH = QE + Src 




(S4) in Section[8] 


B S « - Q ~» 




(S2) in Section[8] 


b, c 


small constants 


Theoreml8.1l 


C 


complex conjugation operator 


Conventionl8.3l 







Conventionl7.3l 


DATA 


data at past null infinity 


Definitional] 


D,D, N,L 


complex frame vector fields 


Definition^. 11 


Dr, B r 


open disk and open ball 


Convention^. 11 


e, A, p 




Definitionl4~T1 


E k Af}(t) 


energy 


Definition O 


<t> = (e, 7, w) 


1Z valued field 


Section[2] 


& = (t,u,v) 


1Z valued constraint field 


Definition[2~4l 




Pole-Flip transformation 


Definitionl3.6l 


I, J 


intervals 




J, J, H, (?*,... 


generic symbols 


Definitions [8Tll8~2l 
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Symbol 


Typical Meaning 


See 


j" 


energy current vector field 


Section|7] 


K 


truncation index, see also '^k' 


Theorem|8.1| 


Ai, A2, A3, A4 


weight functions 


Definition^. 31 


M 


far field ansatz matrix 


Section[5] 




doubly scaled Minkowski field 


Definitions OE21 


M(q,O)9 = h(q,0) 


general symm. hyp. system 


Section|7] 


M(q,E)E = h(q,Z) 


a particular symm. hyp. system 


(S10) in Section[8] 


N 


set of integers > 




No 


set of integers > 




u \°)> <-Aso, 0, t) 


families of subsets of R 3 


{tAt), (KtLU) in oCC. | / | 


7T, TT 


certain permutation matrices 


(S3), (SI) in Sec. [8] 


V 


parity held transformation 


KemarkU.yl 


V, V*, Vh 


generic symbols 


Definitions 15.31 lo.ll 


Qa 


multi-derivative 


Definition 17. 11 


a— 1 
O 


an integration operator 


(15.2b 


9 X , d q , 8q 


gradient operator w.r.t. x, q, q 




Q 
O 


~ ^(afi" — 'si*) 


Definition[5Ul 


q = (q Q ,q\q 2 ,q 3 ) 


coordinates 


Convention^. 11 


q= {q\q 2 ,q 3 ) 






spatial components of q 


Convention^. 11 


Q.K. 


general subsets of R 


(EB4) in Section|7] 


y, al 


particular subsets of R 3 


(S9) in Section[8] 


r> j j 771 

H, H, r 


three pictures 


i - , • 1771 

Section|9] 


r£ 


differentiability index 


Theorem 18. 11 


ft, ft 


real vector spaces 


dm>. dm> 


3?, 3 


real / imaginary part operators 




R, C 


the real and complex numbers 




b 




d4.41i 


Sup^{/}(i) 


supremum norm 


Definition ITT1 


(SHS) , (SHS) , (subSHS) 


symmetric hyperbolic systems 


Section|2] 


Strip(/x, A) 


family of open subsets of R 4 


<i4~n> 


<7 


stereographic chart superscript 


Proposition^. 41 


f 


time coordinate, equal to u + u 


Convention^. 11 




general open subset of R 4 


Section|2] 


E, E» 


-K~ lr R. and Tx~ x Ti, valued fields 


(S3), (SI) in Sectiong] 




coordinates 


Section|2] 




coordinates 


Section|2] 




either (£\£ 2 ) or £ + i£ 2 




z = 5Rz — i 


complex conjugation 




#=(/,w,«) 


ft valued field 


<ED 




formal power series, coefficients 


dO 




truncated formal power series 


(S2) in Section[8] 


= {s,p,y) 


ft valued constraint field 


Propositionl5.1l 


&t(k) 


formal power series, coefficients 


Section|6] 
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B. Generalized Vacuum Equations 

Our main reference for this appendix is |f-Y|. 

The vacuum Einstein equations, written in local coordinates x 1 , x 2 , x 3 , x 4 on a 
connected open set U in M 4 , are a nonlinear system of partial differential equations for 
the ten metric tensor fields g M „. Namely, 

W aiP = (B.l) 

where, 

r>5 d pS d pS , pfJ- pS pp, pS / p[i p^ \ pS 

are the components of the Riemann curvature tensor for the Levi-Civita connection 

p _ 1 I d „ _i d_„ d_ n \ 

J- HvX — 2 V dx^ J vX ' dx" •"/ jA dx^ Jw) 

associated to the metric g^ v . 

There are patent mathematical advantages to introducing more fields and equations 
that, in the presence of appropriate constraints, collapse to the vacuum Einstein equa- 
tions. The purpose of this appendix is to introduce a particular generalized system of 
vacuum equations and explain how it will be used. In this appendix and in AppendixICl 
we work with real quantities. In Appendix iDl we employ a complex tetrad formalism, 
see INPI . as in the main body of this paper. 

Definition B.l. A generalized spacetime is an open subset U q/M 4 together with 

• 16 frame fields E a ^ and the associated vector fields 

It is assumed that E\, E% E3, E4, are frame vector fields. That is, they are linearly 
independent tangent vectors at every point oflA. 

• A constant, symmetric matrix g a b with three positive and one negative eigenvalues. 
For our purposes, the matrix g a b and its inverse g ab are 

/l \ /l \ 

,„ \ 0100 , M 010 

U«0 - 00 -1 ( 9 1 - 00 -1 (B ' 2) 

\00-1 0/ \00-1 0/ 

• 24 connection fields r a b c that are antisymmetric in the indices b and c. 

• 10 Weyl fields W a bk( characterized by 

W abke = -W bak i W ajke + W a£]k + W ak tj = 

W abke = -W aMk g ak W abM = 

Wabki = Wklab 

Convention B.l. Small Latin indices a, b, c ... are frame indices and always run from 
one to four. Small Greek indices A, p, v ... are coordinate indices and also always run 
from one to four. Frame indices are raised and lowered with the constant tensor g ab . 
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We associate to every generalized spacetime 

• A Lorentzian metric g determined by g(E a , E b ) = g a b- 

• A connection V specified by 

g(V Ea E b , E c ) = r abc or, equivalently, V Ea E b = E ab E c 

where r ab c = g cd r ab d- The antisymmetry of r abc in the last two indices, expresses 
the property that the connection V is compatible with the metric. 

• 24 connection torsion fields 

T ab » = r ab c E c » - r ba c E/ - E a {Ef) + E b (E a ») 

They measure the deviation of V from the Levi-Civita connection for the metric g. 
That is, T ab ^ vanishes if and only if 

r abc = \ ( - g(E a , [E b , E c ]) + g(E c , [E a ,E b ]) + g(E b , [E c , £?„])) 

• 36 curvature torsion fields 

Ukiab = E a (r b ek) — E b (r a £k) 

+ r b e m r am k — r a e m r bm .k — (r ab m — r ba m )r m ek — Wuab 

The curvature tensor R k i ab for the connection V is given by 

R k lab = E a (Tw ) — E b (r a £ k ) 

, m t~i k T^mT-i k ( m m rp m\ k 

T J- M 1 am — J- at 1 bra ~ (/ ab ~ 1 ba ~ + ab ) i ml 

= T ab " l r m £ k + U k £ ab + W k £ab (B.3) 

where T ab m E^ n = T ab ^ L . In the event that T ab ^ vanishes, the curvature torsion fields 
Ukiab measure the deviation of the Riemann curvature for the Levi - Civita connec- 
tion from the Weyl tensor Wkiab- 

The curvature tensor Rkiab for the connection V has the symmetries 

Rkiab — —Rlkab Rkiab = ~ Rkiba 

Warning. The customary pair exchange symmetry and cyclic identity do not neces- 
sarily hold when the torsion T ab ^ is not zero. 

• 16 Bianchi fields 

V abi jk = ViWabjk+^kWabij +VjW abki 

or, equivalently, contracted Bianchi fields 

Vbjk = g^Vabijk = V W a bjk 

For the definition of ViW ab jk, see (lB.4t below. The fields V ab ijk vanish if and only 
if the Weyl tensor W a bki satisfies the Bianchi identities with respect to the connec- 
tion V. Similarly, the fields Vbjk vanish if and only if the Weyl tensor satisfies the 
contracted Bianchi identities with respect to the connection V. 
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Remark B.l. The contracted Bianchi fields are equivalent to the Bianchi fields. This 
fact is an immediate consequence of the following algebraic identity. Suppose, A a ujk 
is antisymmetric in the first two indices and totally antisymmetric in the last three. Set 

Aaij = g tk Aabijk A^ ake = \tklijAa^ 

Then A a ujk = § £ijfc £ (^a«, — A Ua + A\ ab ), 

Convention B.l. We define E, r, W, T, U, V to be covariant tensors (vector field 
valued in the case of E and T) on U whose components with respect to the fixed frame 
E a are given by 

E a = E a ^-^ T(E a , Eb) = Tab^-^ 

r(E a , Eb, E c ) = r a bc U (E a , Eb, E c , Ed) = U a b c d 

W{E a ,E bl E c , E d ) = W a bcd V{E a ,E b , E c ) = V abc 

From this perspective: 

• If Xi, X 2 , X 3 are vector fields on U, then T{Xi ,X 2 ,X 3 ) = Xf X* X§ r abc . 
Here X, = X?E a for i = 1, 2, 3. In general, r(X 1 ,X 2 , X 3 ) + g(X7 Xl X 2 , X 3 ). 

• Covariant derivatives of all these tensors are well defined. For example, 

Vj W a bcd = Ei(W a bcd) — r ia m W m bcd ~ r^Wamcd ~ r ic m W a bmd ~ Ei^Wabcm 

(B.4) 

Definition B.2. We refer to the system 

(T, U,V) = (B.5) 
as the generalized vacuum field equations. 

Proposition B.l. A generalized spacetime E a ^, r a bc, W a bk£ on an open subset hi of 
K 4 is a solution to the generalized vacuum field equations ( |B.5t if and only ifW is the 
Levi-Civita connection for the metric g, and the associated Riemann curvature tensor 
coincides with the Weyl tensor W. In this event, g is a solution to the vacuum Einstein 
equations (IB. Il l on U. 

Proposition B.2. The tensors T, U, V have the algebraic symmetries: 

(B.6a) 
(B.6b) 
(B.6c) 

They satisfy the generalized Bianchi equations 

(1, 11,23) =0 (B.7) 

where, by definition, 

Sa M = e a ijk ( ViTj-fc" - U m uj E m " - T^-^Ef ) (B.8a) 

ticab — ' ( ViUabjk — U m ijk r ma b + ^Vabijk — ftrTT Aab ) (B.8b) 

W jk = VbV b jk + U a mab W mb jk - \U m3ab W abm k (B.8c) 

+ \U m kab W a j — \T a b^ Q^JlW a jk 



± ab — 




Vkab 


— —Vkba 


Uabki = 


— Ubakl 


v\ b 


= 


Uabkl = 


— Uabik 


V a bc + Vbca + V ca b 


= 
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Here e a bcd is totally antisymmetric and £1234 = — 1. Furthermore, V, is the tensor 
derivation that acts on frame indices as Vj and ignores coordinate indices. Explicitly, 

v i-L jk ^iy^jk ) 1 ij mk 1 ik ^jm ■ 

Remark B.2. The generalized Bianchi equations ( IB. 7b are identities: they hold for all 
generalized spacetimes. Both ( IB. 51 ) and ( IB. 7b are quadratically nonlinear. Each, has 
exactly one linear term. Respectively, —W^tab and —V a bijk m the equations U = and 
il = 0. The only coordinate index appears in the T = and 1=0 equations. Observe 
that, for fixed EJ 1 , r a f, c , W a bjk, the equations ( IB. 71 ) are linear and homogeneous in 

?a6 , Uabjk and V a bijk- 

Our goal is to construct physically interesting vacuum spacetimes. In this appendix 
we have traded in the 10 traditional metric tensor fields for 50 frame, connection and 
Weyl fields and an additional 76 connection torsion, curvature torsion and Bianchi 
fields. How can this formalism be of any practical use? Not only are there 126 fields, 
but both the generalized vacuum and Bianchi equations are overdetermined, since the 
tensor V vanishes whenever T and U both vanish. 

Here is a rough outline of our strategy. Regard the frame E a ^ and general connection 
2~jjfc as vector fields with values in K. 16 and M 24 respectively. We conceptualize abstract 
gauge conditions as fixed affine linear subspaces £ C K. 16 and Q C K 24 . The frame and 
connection are gauge fixed when E a ^(p) G £ and (p) G Q for all points p EU. No 
conditions are imposed on the Weyl tensor. There are 

dimf + dim£ + 10 < 50 

independent gauge fixed frame, gauge fixed connection and Weyl fields. An abstract 
gauge fixed, generalized spacetime is summarized by a field <P on U taking values in 

£®g®r 10 . 

Abstract gauge conditions should have three properties: 

• Property\B\l. They are always "locally realizable". That is, near each point in every 
spacetime there is a coordinate system and a frame such that the components of the 
frame, with respect to the coordinate vector fields, and the components of the Levi - 
Civita connection, with respect to the frame, lie in the gauge subspaces £ and Q. In 
other words, abstract gauge conditions should not exclude, a priori, any spacetimes. 

• Property\B\2. They are "symmetric hyperbolic". In the present context, a symmetric 
hyperbolic system of partial differential equations for the column vector v is a system 

A a E a {v) = f 

where A a is a symmetric matrix and A 3 + A 4 is strictly positive definite. Now, pick 
bases for £ and Q, and rewrite the 60 equations T = and U = explicitly in terms 
of the components of the field <P. It is required that one can select exactly dim£ 
linear combinations of the 24 connection torsion equations (T = 0) and exactly 
dimCf linear combinations of the 36 curvature torsion equations (U = 0) equations 
and 10 linear combinations of the 16 Bianchi equations (V = 0) which taken together 
comprise a (quasilinear) symmetric hyperbolic system for <P which we will refer to 
as (SHS) . 
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Property iBl 2 is not just wishful thinking, as it may first appear. Only the principal 
parts of dB.51 l, 

E a (E b n - E b (E a n 

Ej(rkba) — Ek{rj a t) 

Ei(W abjk ) + E k (Wabij) + E 3 {W abkl ) 

have to be considered in the quest for symmetric hyperbolic equations. Furthermore, 
only the frame and connection fields in the principal parts have to be written out in 
terms of <P. It is unnecessary to open up and look at the occurrences of the frame fields 
inside the first order differential operators E a . At this level, it is required that there 
are, in turn, linear combinations of the principal parts that are symmetric hyperbolic. 
In principle, the field ^, that contains all information about the generalized spacetime, 
is now uniquely determined, given appropriate data, by (SHS) . However, there is an 
important catch. (SHS) and the abstract gauge conditions imply that some part of the 
tensors T, U and V vanish, but not all. The remaining components are summarized 
in the constraint field If (SHS) is satisfied and <P$ vanishes, then E^, r abc , W ab jk 
is a solution to the generalized vacuum field equations. 
• Property \E\3. They are "dual symmetric hyperbolic". If (SHS) is satisfied, it is re- 
quired that, in an entirely similar way, judicious linear combinations of the general- 
ized Bianchi equations (lB.7t can be brought into the form of a linear, homogeneous 
symmetric hyperbolic system for <P$ which we refer to as (SHS) . In particular, if the 
data for any well posed problem for the system (SHS) vanishes, then the constraint 
field & vanishes everywhere. 

It is much simpler to carry out this general methodology in practice than to formulate it 
in broad conceptual terms. Different problems require different gauges and symmetric 
hyperbolic systems. In Appendix ICl we introduce the wavefront gauge for Lorentzian 
manifolds. In Appendix iDl we fix the abstract wavefront gauge and select symmetric 
hyperbolic subsystems from the generalized vacuum and Bianchi equations that are 
particularly suited to the problem we are solving in this paper. 

Proof (of the Generalized Bianchi Equations ( 1B.7I >). 

• T = 0: Write T(X, Y) = V X Y - V r l - [X, Y]. Repeatedly exploiting the total 
antisymmetry of e a ^ k and then using the Jacobi identity, 

e aijk (ViVjE k - V.-Vi-Efc - V [EtiEj] E k ) 

= e a ^ k {v i (T(E j ,E k ))+T(E i , [E jt E k ])} 

It follows that e ai ^R b ijk E^ = e ai i k {V + T jk b r bl c E c ^ - T jk u -^E^) . Sub- 
stituting ( IB. 3b for the curvature R b ij k , the identity T = follows. 

• il = 0: Apply the operator t a% ^ k Vj to the identity ( IB. 3b in the form 

Rbcjk = Ubcjk + Tj k m r, ncb + W bc j k 

and use the standard Bianchi identity for the curvature tensor corresponding to a 
connection with torsion, e a ^ k : (V \ Rbcjk + Tij e R bc i k ) = 0. 
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• 2J = 0: The divergence \7 b V b jk = - 5 (V a V h - ^b^a)W ab 3k . Express the com- 
mutators in terms of the curvature tensor: 

(v c v d - y d v c )w abJk 

= —Tc/ViWabjk ~ R acdWlbjk ~ R l bcdW aljk — R i jcdWablk — R l kcdW abjl 

Substitute ( |B.3t for the curvature, contract indices, write out the covariant derivatives, 
rearrange, collect terms and cancel to obtain 23 = 0. □ 

C. The Wavefront Gauge for Lorentzian Manifolds 

Here, we introduce the geometric wavefront gauge in the language of Lorentzian ge- 
ometry. The abstract wavefront gauge, in the language of generalized spacetimes, is 
introduced in AppendixiDl 

Proposition C.l (Geometric wavefront gauge). Every point on any Lorentzian 
A-manifold (M, g) has an open neighborhood on which there are coordinates 

and an oriented frame (Ei, E2, -E3, E4) such that g(E a , Eb) = g a b, see dB.2b , such 
that E3 and E4 are both future directed vector fields, and such that 

(a) the coordinate functions u and u are solutions to the eikonal equation, that is, 
g ab E a (u)E b (u) = 0andg ab E a (u)E b {u) = 0. 

(b) the vector field E± is minus the gradient ofu, that is g{Ei, ■ ) = —du. 

(c) the coordinates £ , £ 2 are constant along the integral curves of E4. 

(d) the function = E^u) is strictly positive and the vector field e% E3 is minus the 
gradient of 'u, that is g(esE 37 ■ ) = —du. 

(e) E4 and e^E^ are null geodesic vector fields. 

(f) the frame vector fields E\ andEi satisfy 

g(\7 E4 E u E 2 ) = 

where V is the Levi-Civita connection. 

Proof (Informal). To start with, suppose M is Minkowski space. Let X = (X°, X) = 
(X°, X 1 , X 2 , X 3 ) £ K x R 3 be standard Cartesian coordinates. For every fixed e > 0, 
define u,u : {X e M 4 : X° > — e} — > R by 

u(jp) = sup X 1 , u(jp) = — inf X . 
(-e,x)e/-(p) (-«>x)eJ-(p) 

Here, I~{p) (resp., I + (p)) is the chronological past (future) of p, that is, all points 
that can be reached from p by traveling along past (future) directed, piecewise smooth, 
timelike curves. Note that 

• by construction, u(p) < u(q) and u(p) < u(q) for all q <E I + (p), 

• there is a sufficiently small e > and an open neighborhood U C M of the origin 
X = 0, on which u, u are smooth and du, du are linearly independent. 
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It follows that, du and du are timelike or null on U. 

Suppose du is timelike at a point p 6 U. Then, the level set £ of u passing through 
p would be (locally) a smooth spacelike hypersurface. If q 6 I + {p) is sufficiently close 
to p, then every past directed timelike curve from q intersects S, and 

u(q) = sup X 1 = sup sup X 1 = sup u(p) = u(p). 

(-e,x)e/-( g ) p'esm-(q) (-e,X)e/-(p') p'eznr-(g) 

because I~(q) flfl = U P 'ei:ni-( ? ) n with # = {-^ € M 4 : X° = -e}. 

This contradicts the assumption that du(p) is timelike. Therefore, du is null. Similarly, 
du is null. 

Now, fix a point po on an y Lorentzian manifold M, and let (X° , X) be smooth local 
coordinates that vanish at p. with — dX° a future directed 1-form. Precisely the same 
construction for u and u works on a suitably small neighborhood U C M of p$. 

Define vector fields L and L by g(L, ■ ) = — du and g(L, ■ ) = — du. They are future 
null. Define = L, e 3 = L{u) > and E3 = e^L- In particular, <?(£ , 3, E4) = —1. 

Condition (e) is equivalent to V^L = and V l_L = 0. These are consequences of 
the general fact that for any function w, the acceleration V^l 7 ^ of its gradient field W 
is the gradient field of the function ^g(W, W). 

Let K\ and K2 be spacelike, orthonormal vector fields, perpendicular to £3 and E4. 
Define ( %\ ) = ( _ c °f n Q Q *™% )( ^ ) where a satisfies the differential equation £4 (a) = 
— g(V E 4 Ki, K2) along the integral curves of E4. □ 

Remark C.l. The spacelike vector fields E\ and £2 and the null geodesic vector field 
E\ are tangent to each level set of u. Each level set of u is a union of null geodesies, 
the lines of constant £\ £ 2 and u. 

Similarly, E\, E2 and the null geodesic vector field E3 are tangent to each level set 
of u. Each level set of u is a union of null geodesies. They are, in general, not given by 
the lines of constant £\ £ 2 and u. 

We refer to the intersections of level sets of u and u as wavefronts. They are spacelike 
and their tangent space is spanned by E\ and E^- 

Proposition C.2. Fix a coordinate system and frame as in the geometric wavefront 
gauge of Proposition \C~l\ Let E a ^ be the components of the frame field with respect to 
the coordinate system and let r a b c be the components of the Levi-Civita connection V, 

E a = E a ^ g!^, r a b c = g(V E a Eb, E c ). 

Then 

/* * 0\ 

(*-)-::§: 

\0 e 3 0/ 

where * is a generic symbol. Moreover, 

f* Kn K12 Xl1 Xl2 * \ 




* X2I *22 X21 X22 * 

* * * * 



\0 * * / 
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where the matrix indices (be) run over the ordered sequence 

(12), (31), (32), (41), (42), (34). 

Here, Xab an< ^X AB are me second fundamental forms of the wavefronts in the normal 
null directions E 3 and E4. As such, 

XAB=XBA KaB^KbA' 

Moreover, the entries filled with the generic symbol -j^ satisfy 

-T434 = r 3A i A = 1,2. 

Proof. We first verify the O's and l's in these matrices. The entries of (E a ^) follow 
directly from Proposition IC.ll for example, E 3 (u) = du(E 3 ) = —g{Ei, E 3 ) = 1 by 
(b). The zeros in (-T a bc) are accounted for by (f), by V.E4-E4 = 0, see (e), and by the 
fact that V e 3 E 3 is proportional to E 3 , see (e). Next, A23 - ^213 = (x 12 = X 21 ) 
and r"i24 — I214 = (X12 = X21) follow, by (d) and (b), from [Ei,E 2 ](u) = and 
[Ei , E 2 ] (u) = 0, respectively. Finally, 1^34— /3A4 = follows from [Ea, E 3 ] (u) = 0. 
□ 

D. The Abstract Wavefront Gauge 

In this Appendix, we leave the realm of Lorentz 4-manifolds, and speak exclusively in 
the language of generalized spacetimes, as in AppendixlBl 

We now define the abstract wavefront gauge and show that it has the Properties |B] 1 
through|B]3. 

It is convenient to introduce the complex frame 

(F 1 ,F 2 ,F 3 ,F 4 ) = (D,D,N,L), F a =F a ^ 

where 

D = 2~^(E 1 +iE 2 ), D = 2~i(E 1 -iE 2 ), N = E 3 , L = E 4 . 
These fields are sections of the complexified tangent bundle. 

Convention D.l. Let t be the constant matrix for which F a M = t a E^, and t^ 1 its 
inverse. For every tensor field S in Appendix IB] with components S ai ...o„ " bm with 
respect to the real frame (E a ), the corresponding components with respect to the com- 
plex frame (F a ) are distinguished by boldface letters and defined by 

Sa 1 ...a n 6l -"^=ta 1 il -"*a n ^5i 1 ...i n J * 1 -" iw (r\ 6l ...(t- 1 ) im 6 » (D.l) 

If the tensor field S also carries coordinate indices, they are unaffected by ( ID. 11 1. An 
equivalent statement to dD.lt is: if S is the tensor field for which 

S(E ai ,...,E an ) = S ai ... an b ^ bm E bl ®m---®RE bm 
and S is extended complex linearly in its n arguments, then 

S(F ai ,...,F an ) = S ai ... an bl - b ™F bl ® c ---®cF bm 
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The transformation ( ID. Il l commutes with contraction of indices. Accordingly, the in- 
dices of boldface fields have to be raised and lowered with 



g 



ab 



Sab — tb^ Qij 



see ( 1B.21 i and ( I2.3l l. The complex components S Ql ... ar 



are only introduced for 



notational convenience. The corresponding tensor field 5 will, however, always be real 
in the sense that if all its n arguments are real, then the result is real. 

For the particular covariant tensors E, r, W, T, U, V (see, Convention 
transformation dD.lt becomes 



the 



1 a g x n 



F a 

r a bc = r(F a ,Fb, f c ) 

Wabcd = W(F a ,F b ,F c ,F d ) 
Definition D.l. Let 

<P = (e, 7, w) : U -> K C C 5 



r-p H d 



T(F a ,F b ) 
Uahcd = U (F a ,F b , F c , Fd) 
V abc = V(F a ,F b ,F c ) 



C 5 



(ED 



be a sufficiently differentiable field satisfying conditions (*) and (**) in Definition \2.1\ 
The Abstract Wavefront Gauge Spacetime 

M* = (f ", r ahc , W afccd ) 

« defined just as in Definition \2.1\ 

Remark D.l. Definition IP. 1 1 implicitly fixes abstract gauge conditions in the sense of 
Appendix [B] The affine spaces £ and Q have real dimensions, respectively, 7 and 14. 
That is, the field <P has 3 1 real components. 



Remark D.l. Observe that Proposition IC.2l is the statement that the abstract wavefront 
gauge in Pefinition lP.ll is locally realizable in the sense of AppendixlBlCPropertvlBl 1). 

Proposition D.l. Let M$ be an abstract wavefront gauge spacetime. Let Ai, A2, A3, A4 
be strictly positive weight functions on Li. Then, there are unique fields 



(t,u,o): U 
(t, u, v) : hi 
and such that 



n c c 5 

a c c 5 



see d2TTb 
see d2.6t 



(iti 


it 2 


0\ 


U 


h 





u 


ts 





-tl 


-t 2 


is 


-tl 


-t 2 


*3 


V u 


ts 


-ta 0/ 



u 



(ah)(jfe) 



^1*2 +U2 U 7 - U S U 8 - U7 

u 9 -u 7 -u 6 

-Ug -U 6 -U 7 
— Ul U 4 — M3 

Ul — U3 U 4 

\u 2 - "2 u 7 + U 8 Us + U7 



-U 3 - u 4 u 4 + u 3 u 8 - u 8 \ 

u 5 -u 6 -u 5 

-«6 U 5 -U 5 

-Ui -U2 -U3 + U4 

-U 2 -Ui U 4 - U 3 

-u 3 + u 4 u 4 - U3 u 8 + u 8 / 
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/^(-t>2+«l 

^(»2-«l) + ^«3 
^■(»3 - 03 - "2 + «2 
V ^(-«2+U 2 ) 
1 



T7«2 



) £( 



-04 
1 



«2j 



A: 



-^1 



^(t> a -«i) 



Ai 



if 



"3 



D 4 + 
1 



V2) 
V3) 



' ! 3 



— f- 

A 2 ^ 



2 
2 



"0 + XT U 3 



w l) ^("3 + »3 - «2 - V 2 ) 



The matrix indices (ab), (jk) run over the ordered sequence 

(12) (31) (32) (41) (42) (34) 

Proof. In general, the tensors T, U and V lie (pointwise) in spaces of real dimension 
24, 36 and 16, respectively (see, equations (IB.6I 1). By direct inspection, the equations 



TV = 
U3441 = U4134 



-31 



= 



9U 3 132 =0 



Tab' 1 = 
3U4142 =0 



hold for every field <P, and, consequently, the associated tensors T, U and V lie in 
subspaces of real dimension 24 — 9 = 15, 36 — 4 = 32 and 16 — = 16, respectively. 
By construction, the matrices on the right hand sides of the equations above lie in these 
subspaces. The linear map from (t, u, t>) © (t, u, v) to these matrices has maximal rank 
63, which is 15+32+16. Therefore, the fields (t, u, t>) and (t, u, v) exist and are unique. 
□ 



Proposition ID. II defines a unique "splitting" of the nonzero components of (T, U, V) 
into two sets, (t, u, 0) and (t, u, v). This splitting and the role of the weight functions 
Aj is clarified by the following proposition and the subsequent remarks. 

Proposition D.2. For every choice of strictly positive weight functions Ax, A2, A3, A4 

on IA, the system of equations 

(t,u,o) = 

is equivalent to the system (SHS) , A.(<P)<P = f (<P), in Definition \2.3\ which is a (quasi- 
linear) symmetric hyperbolic system for the field <P = (e,7, w) provided e% > 0. In 
particular, the abstract wavefront gauge of Definition \D.l\ has Propertv\E\2. 

Proof. The equivalence of (t, u, t>) = with the symmetric hyperbolic system (SHS) in 
Definition ^. 3 l is by direct (machine) verification. □ 

Remark D.3. The term "symmetric" is a slight misnomer, in the sense that the matrices 
(<P) determining the principal part (<P) are complex Hermitian rather than 
real symmetric. 

Remark DA. (SHS) is of a form which is particularly suited to constructing solutions. 
The reason is that the first two blocks Ai($), A2(^) of the principal part, correspond- 
ing to the principal parts of t = and u = 0, are diagonal and only L or N appear. 
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Remark D. 5. Note that 



Ai x 



Principal Part Operator 



•02 
W 



N 


D 


D 


L ' 



_ A 2 x 



jV 


D 


A3X 


D 


L~ 


N 


D 




D 


L 



A4X 



V 



N 


D 


D 


L 



W4 



The dotted lines in the schematic diagram for the 5x5 matrix on the right hand side 
indicate that the overlapping entries are the sums \L + Ai+iiV. Each ( ^ ^ ) block is 
symmetric hyperbolic, and consequently, so is the 5x5 matrix for any choice of the 
strictly positive weight functions. 

Remark D.6. The weights have a natural interpretation in terms of energies. The en- 
ergy current naturally associated to (SHS) , A(<P)<P = f(^), is the vector field = 
^A^(<P) <P. Estimates are obtained by applying the divergence theorem to 

dpf = #t(0 M A")* + 20?(<2> t f (<£>)). 

The energies are integrals over the spacelike components of the boundary. The functions 
Ai, A2, A3, A4 appear in the boundary integrals and play the role of weights for the 
components wi , w 2 , W3 , w 4 , W5 . 

Proposition D.3. Let A\, A2, A3, A4 be strictly positive weight functions on U. The 
components of (t, u, v) are given by ( 12.71 l in Definition \2.4\ The field = (t, u, v) : 
hi — ► 1Z is called the constraint field associated to <P = (e, 7, w). 

Proof. By direct (machine) calculation. □ 



The generalized vacuum field equations ( IB. 51 reduce, in the abstract wavefront gauge, 
to (SHS) and = 0, see Proposition lD.il How can we ensure that a solution to (SHS) 
also satisfies = 0? The answer is given in Proposition |D3] 

Proposition D.4. Assume that <P = (e, 7, w) satisfies > and solves (SHS) or, equiv- 
alently, (t, u, t>) = 0. Let the Latin indices in X a ^, iX a bc, %3ab denote components of 
the fields dB.81 > with respect to the complex frame field F a (see, Convention \D. 1 h The 
subsystem of the generalized Bianchi equations (IB.7I > given by 



( iUu \ 

2(^-412 + it-434) 

il 2 14 
iliu 
il-223 

2 (-U112 + i&134) 
i(il 212 + il 234 ) 
-^332 



X 41 
i(2J 12 +2J 3 4) 
»23 



= 



(D.2) 



is equivalent to the system (SHS) , A(<P)^ = f d x <P)<I>K in Proposition \2.3\ which 
is a linear homogeneous symmetric hyperbolic system for = (t, u, v). In particular, 
the abstract wavefront gauge of Definition \D.l\ has Propertv\B~\3. 



89 



Proof. By assumption, (t, u, t>) = 0. The equivalence of ( ID. 21 ) with the linear, homoge- 
neous symmetric hyperbolic system (SHS) is by direct (machine) verification. □ 

E. Symmetries: Proofs 

In this section, we prove that the field transformations introduced in Section [3] are field 
symmetries (see, Definition 13. lb . 

Recall from Section [3] the definition of a field transformation S. In this Appendix we 
take a slightly different perspective and regard i£M and x' 6 W as two sets of (global) 
coordinates on the same 4-dimensional manifold. Similarly, we regard <P, A, F a , V, W 
and their primed counterparts as objects on this 4-manifold. Here, V is the connection 
associated with r, and W is viewed as a 4-covariant tensor (see, Appendix [B] and 
Convention ID. 11 1. Field transformations are defined, in this Appendix, through their 
action on the coordinates x, the complex frame vector fields F a , the connection V, the 
4-covariant tensor W, and the strictly positive weight functions A. The definitions of 3, 
21 given below are equivalent to the corresponding definitions in Section [3] Those of £, 
3 are slight generalizations, because <£ , C 2 and £ are allowed to depend on (x 1 , x 2 , x 4 ). 

Convention E.l. For the rest of this appendix, it is implicitly assumed that x' = S ■ x 
is a local diffeomorphism on K 4 . With this understanding, it is unnecessary to specify 
the ranges of the x and x', because the discussion is purely algebraic. A dot, ■ , always 
denotes a group action. 

Angular coordinate transformation £ Let x 2 , x ), <L (x , x 2 , x 4 ) be functions. 



x 



{€ 1 {x 1 ,x 2 ,x A ), € 2 (x\x 2 ,x 4 ), x 3 , x 4 ) 



€■ V = V 




£• W = W 



€-A 



A 



U(l) transformation 3- Let ( = ((x 1 , x 2 , x A ) £ U(l). 
x = 3 • x = X 



3- v = v 
3 • W = w 




3- A 



A 



Global Isotropic Scaling 3- Let Z > be a constant. 



X - -^j ■ x ; X j -\JX j \§X J 



a- v = v 
3 • w = z 2 w 




Z-A 



A 





21- V = V 



diag(l,2l 2 ,2l 4 ,2t 6 )yl 
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Remark E.l. The action on the frame induces a global conformal transformation of the 
associated metric: € ■ g = g, 3 ■ g = g, 3 • g = Z 2 g, 21 • g = 2l 2 g. 

Remark E.2. €, 3, 3, 21 preserve the wavefront gauge and, consequently, induce an ac- 
tion on <P = (e, 7, w). We illustrate this important fact by three examples. First, 

( a ' F ) 3 M 3WY = *~* F ^ ^ = %~* afr + ^ & + 

compatible with the wavefront gauge. Necessarily, (21 • e)j = 2l -3 for i = 4, 5. 
Second, abbreviating F' a = 2t • F a , we have 

(21 • «;) 6 = (21 ■ W)(Fl, F' 2 , F%, F$ = 2T 6 (21 • W)(F 3 , F 2 , F 3 , F 2 ) 

= 2T 6 2t 2 W(F 3 , F 2 , F 3 , F 2 ) = 2T 4 w 5 . 

Third, abbreviating F' a = 3 ■ F a 

(3 • r){F[, F{, F 2 ) = (3 • <?)((3 • V) n Fi,Fi) =(g(V Fl F 1 ,F 2 )+F 1 (0 
(3 • r)(F>, F3, F4) = (3 • <?)((3 • V)^F^,Fi) = Cg(V Fl F 3 ,F A ) 

which is consistent with the wave front gauge if and only if (3 • 7)3 = C73 + \ -Pi(C) 
and (3 • 7)4 = C _1 74 + \ -^(C -1 )- By direct calculation, the present definitions of £, 
3, 3, 21 are seen to be equivalent to those in Section [3] (generalizations in the case of £ 
and 3). 

Proposition E.l. Let S be one of€, 3, 3, 21. Then, separately: 

• (*) and (**) are preserved (see, Deftnition \2.1j . 

• (t,u,t>) = if and only if {S ■ S ■ u, 5 • t>) = 0. 

• (t, m, u) = if and only if ' [S ■ t, S ■ u, S ■ v) = 0. 

/n particular, S is a field symmetry in f/ie sense of Definition \3.1\ 

Proof. Let Ricm be the Riemann curvature tensor associated to g, considered as a 4- 
covariant tensor. For each S, there are complex functions m, n 2 , K 3 , K4 and a constant 
ft > such that KiK 2 f2 2 = K3K4J? 2 = 1 and 

(S-F«)"g$jjr = K.Fa /1 5 iir a = 1,2,3,4 
S • V = V S • ,g = I2 2 g S ■ Ricm = I2 2 Ricm S ■ W = fi 2 W 

We abbreviate F' a = S ■ F a = n a F a . By the definition of T, U, V in Appendix|B] 

(5 • TYiF'^FDggy = K aKb T^(F a7 F b )^ (E.l) 
(S ■ U){F' a , Fl,F' c , F' d ) = Q 2 n a K b K c K d (lI(F a , F b , F c , F d ) (E.l) 

+ g {F a ,F b )T»{F c ,F d )±^) 
(S ■ V){F' a ,Fl,F' c ) = K a K b n c V(F a , F b , F c ) (E.3) 
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By the definition of (t, u, t>) and (t, u, v) in Proposition lD.il 





€ 


3 


3 


21 


t = 


<= 


S 


■ t = 


and t 


= 4= 


=> S 


t = 







true 


true 


true 


u = 


o <s= 


=> s 


u = 


and t 




=> S 


t = 





true 




true 


true 


= 


<= 


=> s 


= 


and v 


= <= 


=> S 


■ v = 


= 


true 


true 


true 


true 



where 

• in the case of t and t we use ( IE. II ). observing that d ^.,y is proportional to for 
H = 1,2, 3, 4 if S is one of 3, 5, 21, 

■\ tl-ka (iiirQ /"\ T ii 'in/I ni >ira hpq /II^j /n ci ari t-i r\ rr f V» o t 

dx f 



in the case of u and it we use ( IE.2I ), observing that = for a,fx= 1, 2, 3, 4 if 



S is one of £, 3, 21, 
» in the case of o and w we use (IE. 3b and the transformation law for A. 



The remaining cases are discussed separately. Case 1: If S = £, we have 

ae 1 a , ac 2 a /„ _ i o\ a _ a 



0x" 



ax" a{x'Y "t ax" a(Fp 



such that t = <^=> € ■ i = and t = < 
Case 2: If 5 = 3, note that for all a, 6, c, d, 

1 8k 



d(x'y 

<t-t = follow from Proposition lD.il 



g(F a ,F b )T»(F c ,F d )±%% = g(F a ,F b ) £. =12 T*(F C , F d )£§^ 



by the structure of the torsion-matrix (last column vanishes) and 



Oh.,, 
i).r"' 



= 0. Because 



of K3 = tii = 1 and the term g(F a ,F b ), the expression vanishes unless (a, b) 
(1, 2), (2, 1). At this point, one verifies directly that 



(t,u) = 
(t,u) = 



3-u = o, 
3-« = o, 



(3-t,3-u) = o =>■ u = o, 
(3 • t, 3 ■ «) = o u = o 



This concludes the proof. □ 



F. An Estimate for Pole-Flip 

Lemma F.l. Let B c R 2 fee o/?en. For a/Z < ri < T2 /ef 

A(n,ra) = {(&«.«) G K 4 : (a,u) eB, n < |£| < r 2 ). 
For all \a\ > 1, a/Z integers R > 0, and all C -fields <f> : A(|a|ri, | or | r"2 ) ~~ * 7£> 

|| Fn Pa ' ^|| c n (j 4(M, M)) ~(H,ri,r 2 ) ||^||c R (A(|a|ri, |a|r 2 )) ( F 'l) 

For FlipQ,, see Deiinition \3.6\ The same estimate holds for the C R norms on the image 
of the sets A('-^, ^j-) and A(\a\ri, \a\r2) under the change of coordinates from x = 
(£, u, u) to q = (t, £,u), see Convention ]? '.71 

Remark F.l. The point here is that the constant in ( IF. lb is independent of |a| > 1 and 
B. Lemma IeTI is used in Step 9 of the proof of Theorem 18. II If Theorem 18.11 would 
be stated with the additional condition a = 21, then Lemma IeTI would not appear in its 
proof. 
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Proof. Without loss of generality, we can assume a > 1, because Flip a = Flip_ a . 
Let € a be the angular coordinate transformation (Definition 13. 21 ) given by <£(£) = a£. 
Let $ be the angular coordinate transformation with (£(£) = Let 3 be the f/(l) 
transformation (Definition [3~3T l with C(£ ) = — £/£. We have Flip a = SoC^o^o^i/a = 
£ a o3°5° £i/a- Decompose <P = (e,j, w) = <S>\ © ^ 2 where <Pi = (e 1; e 2 , e 4 ,e 5 ) 
and = ( e 3i 7ii ■ • ■ , 78j w 1j ■ ■ ■ > w s)- Introduce the notation e(l) = 1 and e(2) = 0. 
For all (3 = (pi,fo,0, 0) £ N 4 , with \f3\ < R, all < a < b, all A > 0, all < r < R, 
all N = 1,2, and all fields <£, 

||^(CA • #)jv||c"(A(Aa,A6)) = A^Hfl || || C o (A(a>6)) 

\\(£\-$) N \\cr(A(\a,Xb)) < A e(Ar) - r |l^||c"(A(a, b) ) for 0<A<1 

|| (ff- ^)jv||c*(A(6- 1 ,o- 1 )) $(H,a,6) ||*iv||c"-(A(a,6)) 
11(3 • *)l||cf(A(o,6)) $(ii,o,6) ll*l||c»-(A(a,6)) 

||(3 ■ ^)2||c(A( a ,6)) <(H,a,6) Pi || C"*(A(a,6)) + II ^2 || C r (A(a,b)) 

Set X = (R, n, 7*2)- The above estimates imply 
||^(FHp a -$) 2 || C o (A(M)) 

= H^((CaO3offoCi)-#) 2 || C (A( ^ )) 

= a- |/3| ||^((3°^oe:i).<?) 2 || c o (A(7 L j j_ )) 

<X OT^ ||(£i • $)i\\cMWt u t*)) + a- lPl \\(€i. • *) a ||ol/»l(A(r 1 ,r,)) 
<X OT \\$l\\c\0l(A(ari,ar 2 )) + 11*2 1| c^l (A(ari ,ar 2 )) 
<X ||*||olPI(A(ari,Q!ra)) 

Similarly, ||9^(FUp Q • <P)i\\ c »(A(^,^)) <x Pllc!0i(A(an,m- 2 ))' Together, 

||3' 3 (Flip Q •<?)|| c o (A( ^^ )) < A - ||*|| |/»|(A(ar ll ar !1 )) 

The last estimate, and the fact that and both commute with Flip Q , imply 
Lemma IrTI □ 



G. Supplement to Proposition §□] 

Here we make explicit all the polynomials J, £T, Ti, T-i (see, Definition |8.U that appear 
in Proposition 18. II In particular, it will be clear, by inspection, that these polynomials 
are independent of K, as required. We use the C 5 © C 9 © C 4 block-notation, as in (S6) 
and (S7), and the complex conjugation operator C. 

• Recall that, for each fj,, the square matrices (q, 0), -B 2 (q, 0) and B^ (q, 0) are 5x5, 
9x9 and 4x4, respectively. See (S6). 
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In equation dOat . we have B?(q,Q) - B? = u~ 2 Qk (that is, H = 0) for all 
i = 1,2,3 and \i = 0, 1, 2, 3, except in the cases (i, /a) = (1, A) with A = 1,2, when 



51(9,0)-^ = 



and similar for B\ (q, 0) 
For (IQbll . 
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where the operators £ + tp and £_ tp are defined by ip 
tively, where ip is any complex valued function. 
For(l83dli 



- ip C and i(ip — ip C), respec- 




Q33 — Q33 

1 +C 

-C 1 





2 
10 

0, 



For ( I8.7bb , ( I8.7cl ), ( I8.7db . recall from Remark 1531 that f(q, Of) is a quadratic poly- 
nomial in >P, & without constant term. Let f (q, \P) = fm (q, &) + f(2) (q, be its 
decomposition into homogeneous (over R) parts. By definition ( I8.3cb . 



-ill 

ds I 



__ Q(q,s~)II 



dsi ds 2 I (si ,s 2 )=0 



(2) 



(q, 7r(si^ + s 2 II)) 
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It follows from direct inspection of ( I5.8bb that f( 2 ) (q, W) is a polynomial in \P, \P 
whose coefficients are Laurent polynomials in -, with complex coefficients. Now 
one reads off from ( 15.101 ) that the Laurent polynomials have the structure recorded in 
diHbl), d877cl i and (fOdl . 

For ( I8.7ab . recall from Remark 1531 that A(q, \P) is affine linear (over R) in <P. Let 



A M (q, ty) = A^ (q, ty) + (q, W) be its decomposition into homogeneous parts. 
By d8.3bb . we have | s=0 B^(g, sE) = 7r _1 A^ (g, 7r S)7r. Writing out the result 
in the notation (h, a, £) — ir (Si, S2, £3) of ( 18.21 ). we obtain for /i = and fj, = 3, 



d.s 



s=0 B^(g, S H) =diag(0, £/i 3 , 
and for /x = A = 1,2, 



c&s, TjUs) © (7^3 1 9 ) ©0 4 



dsls=0 



(q,sS) 







1 h A 
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_+3 







r/lA+3 











© 9X 9 © {^h A+3 1 4 



References 

AnRe. Anderson, L. and Rendall, A.D., Commun. Math. Phys. 218, (2001) 479-5 1 1 

BK. Belinskii, V.A., Khalatnikov, I.M., Sov. Phys. JETP36, (1973) 591-597 

BKL. Belinskii, V.A., Khalatnikov, I.M., Lifshitz, E.M., Adv. in Phys. 31, (1982) 639-667 

BBM. Bondi, H.. van der Burg, M.G.J, and Metzner, A.W.K., Proc. Roy. Soc. Lond. A 269, (1962) 21-52 

Cha. Chandrasekhar, S., The Mathematical Theory of Black Holes (Oxford U., 1983) 

Chr. Christodoulou, D., The Formation of Black Holes in General Relativity (EMS, 2009) 

ChrKl. Christodoulou, D. and Klainerman, S., 

The Global Nonlinear Stability of the Minkowski Space, (Princeton U., 1993) 
Fr. Friedrich, H„ Proc. R. Soc. Lond. A 375, (1981) 169-184 
John. John, F., Partial Differential Equations (Springer, 1978) 139 
NP. Newman, E. and Penrose, R., J. Math. Phys. 3, (1962) 566-578 
Pen. Penrose, R., Phys. Rev. Lett. 14, (1965) 57-59 
Sze. Szekeres, P., J. Math. Phys. 13, (1972) 286-294 

Tay. Taylor, M. E„ Partial Differential Equations 111 (Springer, 1997) 360-370 
Tay2. Taylor, M. E., Partial Differential Equations I (Springer, 1996) 434-439 



